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COMPUTER PROGRAMS FOR CALCULATING THE STATIC 
LONGITUDINAL AERODYNAMIC CHARACTERISTICS 
OF WING-BODY-TAIL CONFIGURATIONS 

by Michael R. Mendenhall, Frederick K. Goodwin, 

Marnix F. E. Dillenius, and David M. Kline 
Nielsen Engineering & Research, Inc. 

SUMMARY 

This document is a user 1 s manual for four computer programs developed 
to calculate the longitudinal aerodynamic characteristics of wing-body and 
wing-body-tail combinations. The R1307 program is based on a linear method 
and is limited to the small range of angles of attack for which the lift 
and moment characteristics of wings and bodies are linear with angle of 
attack. The CRSFLW program is based on a crossflow method of predicting 
the forces and moments on bodies alone or wing-body combinations over a 
large angle-of-attack range. The method states that the normal-force 
distribution on a body is made up of a potential term given by slender- 
body theory and a viscous crossflow term modified by Newtonian theory. The 
SUBSON program predicts the longitudinal aerodynamic characteristics of 
wing-body-tail combinations at subsonic speeds and at angles of attack for 
which symmetrical pairs of vortices are shed from the body nose and the 
leading and side edges of the lifting surfaces. Program SUPSON predicts 
the longitudinal aerodynamic characteristics of wing-body-tail combinations 
at supersonic speeds in the same angle-of-attack range. 

This program manual contains a description of the use of each program, 
instructions for preparation of input, a description of the output, program 
listings, and sample cases for each program. 

INTRODUCTION 


An engineering prediction method for determining the longitudinal 
aerodynamic characteristics of wing-body-tail configurations including 
nonlinear aerodynamics of components and interference between them is 
presented in reference 1. Specifically, particular attention is paid to 
the nonlinearities associated with symmetrical vortex shedding from the 
nose of the fuselage and with leading-edge and side-edge separation 


vortices from the lifting surfaces. Four computer programs were developed 
to calculate the longitudinal aerodynamic characteristics of wing-body and 
wing-body-tail combinations. This document is a user's manual for these 
four computer programs. Principal reliance is made herein to reference 1 
for a description of the methods and calculation procedures. Reference 1 
also contains calculated results and comparisons with data for various 
types of configurations. 

The first program, called R1307, is based on the method of reference 2. 
Use of this program is limited to the small range of angles of attack for 
which the lift and moment characteristics of wings and bodies are linear 
with angle of attack. In most cases, the upper angle limit of usefulness 
of R1307 is approximately 10°. 

The second program, called CRSFLW, is a crossflow method of predicting 
the forces and moments on bodies-alone or wing-body combinations over a 
large angle-of-attack range. This program is based on the method described 
in references 3 and 4 which states that the normal-force distribution on a 
body is made up of a potential term gi'ren by slender-body theory and a 
viscous crossflow term modified by Newtonian theory. Nonlinear effects due 
to vortex shedding from the nose of the fuselage or the leading edge of 
the lifting surface are not included in this method. 

The third program, called SUBSON, predicts the longitudinal aerodynamic 
characteristics of wing-body-tail combinations at subsonic speeds. It is 
based on the extension of the method of reference 2 to angles of attack for 
which symmetrical vortices are shed from the body nose and the leading and 
side edges of the lifting surfaces. The lifting surfaces are described by 
a vortex-lattice scheme. 

The fourth program, called SUPSON, predicts the longitudinal aerodynamic 
characteristics of wing-body-tail combinations at supersonic speeds. It 
uses the same procedure as program SUBSON with the exception that a constant- 
pressure-panel lifting-surface theory describes the wing and tail surfaces. 

Each of the above programs is presented in a separate self-contained 
section of this manual which contains a description of the use of the 
program, instructions for preparation of input, a description of the output, 
a complete program listing, and sample cases. A common list of references 
follows the fourth section. 
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PART I - R1307 COMPUTER PROGRAM 


Introduction 

This computer program automates the method presented in NACA Report 
1307, reference 2, for determining the lift, pitching moment, and center 
of pressure of wing-body-tail combinations. The method is restricted to 
bodies of circular cross section with wings and tails which do not have 
swept-forward leading edges or swept-back trailing edges. It is further 
restricted to small angles of attack and small angles of. wing and tail 
incidence in which the forces are linear with angle. For a complete 
description of the method the user of this program should consult 
reference 2. 

The following sections of this write-up will present a description 
of the program, a description of the input, a description of the output, 
a program listing, and sample cases. The notation used is that of refe- 
rence 2. The list of symbols from reference 2 is included herein for 
reference purposes. 

The program is written in Fortran IV for the IBM 360 series machines. 
No tapes, drums, or disks other than the standard input/output units are 
required. The running time for a typical case on the IBM 360/67 is two 
to three seconds. To run the program on other machines such as the 
CDC 6600 minor changes are required since the inverse sine and cosine 
routines are used. On the IBM 360 these are called by ARSIN and ARCOS 
while on the CDC 6600 they are ASIN and ACOS . One or both of these 
routines are used in the following subprograms: 


EQ19 

EQ24 

EQ24L 

EQ30 

EQ30L 

VOLOG 


List of Symbols 
tail-alone aspect ratio 




wing-alone aspect ratio 


c 


r 


chord at wing-body juncture or tail-body juncture 
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tip chord of wing or tail 

lift coefficient based on wing-alone area except tail-alone lift 
coefficient based on tail-alone area 

lift-curve slope for angle of attack, per radian 

lift-curve slope for wing or tail incidence, per radian 

pitching-moment coefficient based on wing-alone area 
pitching-moment-curve slope for angle of attack 

pitching-moment-curve slope for wing-incidence angle 


body diameter 
wing vortex semi span 

ratio of lift component to lift of wing alone or tail alone for 
variable wing or tail incidence 

ratio of lift component to lift of wing alone or tail alone for 
variable angle of attack 

ratio of lift of body nose to lift of wing alone 


length of wing-body-tail combination 

distance from most forward point of body to intersection of wing 
leading edge and body 

distance from most forward point of body to center of moments 


reference length 


distance from most forward point of body to shoulder of body 
nose 


distance from most forward point of body to intersection of 
tail leading edge and body 


distance from most forward point of body to center of pressure 
position 


lift force 

cotangent of leading-edge sweep angle 





N 


U W 


N 


"W 


x,y,z 

X 


*h 


e 

SA 

6 

X 


N 

T 


free-stream Mach number 
body radius 

body radius at shoulder of nose 
body radius at wing 
body radius at tail 

maximum semispan of wing or tail in combination with body 
cross-sectional area of nose at maximum section 

reference area 

tail-alone area 

wing-alone area 

volume of body nose up to shoulder 

streamwise, spanwise, and vertical coordinates, respectively 

distance to center of pressure measured from intersection of 
wing leading edge and body for wing quantities and from inter- 
section of tail leading edge and body for tail quantities 

distance from intersection of wing leading edge and body to 
wing hinge line 

angle of attack of body centerline 

VK - ~\ 

wing-alone or tail-alone effective aspect ratio 
wing-or tail-incidence angle, degrees 

taper ratio, 

sweep angle of leading edge, degrees 

Subscripts 

body nose 
tail 
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w 


wing 


B (T) body in presence of tail 
B (W) body in presence of wing 
T(B) tail in presence of body 
W (B) wing in presence of body 


Description of Program 

The R1307 computer program consists of a main program, thirteen 
function subprograms, and four subroutine subprograms. The main program 
performs the calculations as shown on the calculating form presented in 
Table I of reference 2. All input and output takes place in the main 
program. The function and subroutine subprograms provide the quantities 
which are obtained from the charts and equations of reference 2 when 
following the procedure outlined in Table I of that reference. Unless 
otherwise noted, equations, figures, charts, and appendices referred to 
in this section are those in reference 2 . 


Function APENB calculates the tail interference factor using the 
equations of Appendix B. These values are plotted in Chart 7. 

Function APENC calculates the tail interference factor using the 
equations of Appendix C. This function is used at supersonic speeds for 
rectangular tails when the wing vortex is inboard of the tip of the tail. 

Function CHRT 8 calculates the wing or tail lift-curve slope at super- 
sonic speeds from the curves of Chart 8 . Values of 3 c L a from the curves 
are tabulated in CHRT 8 and linear interpolation in 3A, X, and zero sweep 
location is performed. 


Function EQ14 
tion (14) . Values 

Function EQ19 
tion (19) . Values 

Function EQ21 
tion (21) . Values 

Function EQ24 
using equation (24) 


calculates the value of ^( 3 ) or k t(b) us i n 9 equa- 
obtained from this equation are plotted in Chart 1. 


calculates the value of ^( 3 ) or k i>(B) us i ni 3 equa- 
obtained from this equation are also plotted in Chart 1. 

calculates the value of Kg ^ or ^ (t) us i n 9 equa- 
obtained from this equation are plotted in Chart 1. 

calculates the value of Kg (w) ( 3 C La ) w or Kg (T) ( 0 C La ) T 
. This function is used for the high-aspect-ratio range 
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at supersonic speeds when there is an afterbody behind the wing or tail 
and 01^ or 01^ is greater than one. Values are plotted in Chart 4(a). 

Function EQ24L calculates the same quantity as EQ24 for the case 
where there is no leading-edge sweep, 01^ = « or 01^ = oo. The limiting 
form of equation (24) for this case is 



Function EQ26 also calculates the same quantity of EQ24 but for 
values of 011^ < 1.0 or 0m T < 1.0. Equation (26) is used for this calcu- 
lation and values are plotted in Chart 4(a). 

Functions EQ30, EQ30L, and EQ31 are analogous EQ24, EQ24L, and EQ26, 
respectively, for the case where there is no afterbody behind the wing or 
tail. Equations (30) and (31) are used for the calculations and values 
are plotted in Chart 4(b). The limiting form of equation (30) when 
01^ = oo or 0m T = oo is 



Function V0L0G calculates the volume of an ogive nose using the equa- 
tions in the appendix of reference 3. In the notation of reference 2, 
the volume is 
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where 



Subroutine CH1011 calculates the wing-or tail-alone center of pressure 
at supersonic or subsonic speed. Curves for determining this quantity are 
given in Chart 10 and Chart 11 of reference 2. Values obtained from 
these curves are tabulated in the subroutine and linear interpolation in 
PA, X, and zero sweep location is performed. 

Subroutine CH1416 determines the center of pressure of the lift 
transferred to the body by the wing or tail. This quantity is presented 
in Charts 14, 15, and 16 of reference 2. For large aspect ratios at 
supersonic speeds. Chart 14 is used. For all other supersonic cases, 

Chart 15 is used. For all subsonic cases, Chart 16 is used. Values have 
been read from these charts and are tabulated in the subroutine. The sub- 
routine selects the chart to be used on the basis of the Mach number and 
aspect-ratio parameter and then determines the value of (x/c^).,-.. or 
(x/c r ) B ^ T j by linear interpolation in the tables. 

Subroutine CH56 calculates the lateral position of the wing vortex 
(f - r) w / (s - r) at the wing location using Chart 5 for subsonic speeds 
and Chart 6 for supersonic speeds, values have been read from these 
charts and are tabulated in the subroutine. Linear interpolation in PA, 

X, and zero sweep location is performed. 

Subroutine KFACT calls the appropriate function subprograms used to 
calculate K^g^, K 0 ^ , and or Kg^, and k T ( B j. After 

these are determined, equation (33) of reference 2 is used to calculate 

(W) ° r k B (T) ' 


Description of Input 

Variable definitions .- The format of the input cards for the R1307 
program is shown in figure 1. In this figure the program variable name 
is shown as well as the card columns in which the value is punched and 


8 


the format in which it is punched. The remainder of this section consists 
of a table listing these program input variables along with the algebraic 
symbol used in reference 2, if applicable, and the input variable definition. 
The algebraic notation used in defining the configuration is shown in 
figure 2. A discussion of the preparation of the input is presented in 
the section following the table. All input length and area quantities 
are dimensional and should have consistent dimensions. 


PROGRAM ALGEBRAIC 

NOTATION NOTATION 


DEFINITION 


Item 1 

NHEAD Number of cards of information which 

identify run. 


Item 2 
HEAD 
Item 3 
NTAIL 

NOSE 

NAFTBW 

NAFTBT 


Identifying information. 


Is a tail present? 

NTAIL = 0; no 
NTAIL = 1; yes 

Ogive nose? 

NOSE ■ 0; no 
NOSE » 1; yes 

Afterbody behind wing trailing edge? 
NAFTBW - 0; no 
NAFTBW ■ 1; yes 

Afterbody behind tail trailing edge? 
NAFTBT « 0; no 
NAFTBT = 1; yes 


Item 4 

FMACH 

SLM 


NL 


M 


Eree-stream Mach number. 

Distance from most forward point of body 
to center of moments. 


REFS 


Reference area to be used in calculated 
lift and moment coefficients. 


REFL 


Reference length to be used in calculated 
moment coefficients. 
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PROGRAM ALGEBRAIC 


NOTATION 

NOTATION 

DEFINITION 

Item 5 



RW 

r w 

Average body radius at wing location. 

SW 

s w 

Maximum semispan of wing in combination 
with body. 

CTW 

C t 

Tip chord of wing. 

CRW 

c r 

Chord at wing-body juncture. 

WLESWP 

^LE 

Sweep angle of wing leading edge, degrees. 

SLW 


Distance from most forward point of body 
to intersection of wing leading edge 
and body 

XHW 

x h 

Distance from intersection of wing leading 
edge and body to wing-hinge line 

Item 6 



RT 

r T 

Average body radius at tail location. 

ST 

S T 

Maximum semispan of tail in combination 
with body. 

CTT 

C t 

Tip chord of tail. 

CRT 

c r 

Chord at tail-body juncture. 

TLESWP 

^E 

Sweep angle of tail leading edge, degrees. 

SLT 

^T 

Distance from most forward point of body 
to intersection of tail leading edge and 
body. 

Item 7 



RN 

r N 

Body radius at shoulder of nose. 

SLS 


Distance from most forward point of body 
to shoulder of body nose. 

VS 

v s 

Volume of body nose up to shoulder. 
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PROGRAM ALGEBRAIC 


NOTATION 

NOTATION 

DEFINITION 

Item 8 



CLAW 

< C L a )w 

Wing lift-curve slope based on exposed 
wing area, per degree. 

CLAT 

(CL a* T 

Tail lift-curve slope based on exposed 
tail area, per degree. 

CLAN 

^° L a^N 

Nose lift-curve slope based on body 
cross-sectional area at shoulder of 
body nose, per degree. 

Item 9 



ALFI 

a i 

Initial body angle of attack for which 
calculation is to be performed, degrees. 

ALFF 

a f 

Final body angle of attack for which 
calculation is to be performed, degrees. 

DALF 

Aa 

Angle-of-attack increment to be used 
between and a^, degrees. 

Item 10 



DWI 

V 

X 

Initial wing incidence angle for which 
calculation is to be performed, degrees. 

DWF 


Final wing incidence angle for which 
calculation is to be performed, degrees. 

DDW 


Wing incidence angle increment to be used 
between 6^ and 6 Wf , degrees. 

Item 11 



DTI 

6 T. 

1 

Initial tail incidence angle for which 
calculation is to be performed, degrees. 

DTF 


Final tail incidence angle for which 
calculation is to be performed, degrees. 

DDT 

A6 t 

Tail incidence angle increment to be used 


between 6«p.j_ and 6 Tf , degrees. 

Input preparation.- A discussion of the input variables will be 
presented in this section as an aid in the preparation of the data deck. 
Before beginning this discussion a few words need to be said as to what 
the computer program treats as the wing and what it treats as the tail. 
If a configuration has one set of lifting surfaces, this is the wing 
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regardless of its axial location on the body, and data describing the set 
are input as wing data. If there are two sets of lifting surfaces, the 
set closest to the nose is the wing and the aft set the tail. For example, 
if the configuration has a set of canards near the nose and a wing further 
aft on the body, the canard data are input as wing data and the wing data 
input as tail data. 

Item number 1 of the input data is an index NHEAD which indicates 
how many cards of information, item number 2, are to follow to identify 
the run. The value of NHEAD must be one or greater. Item number 2 is a 
set of NHEAD cards containing Hollerith information which the user wishes 
to use to identify the run. This information can be punched anywhere in 
the cards and is reproduced in the output just as it is read in. 

Item number 3 contains four indices. The first, NTAIL , specifies 
whether a tail is (NTAIL =1) or is not (NTAIL = 0) present. The second 
index, NOSE, specifies whether the nose is (NOSE =1) or is not (NOSE = 0) 
an ogive. The purpose of this index is to provide a computation within 
the program of the nose volume for an ogive. For non-ogive noses, the 
volume must be input in Item 7. The third index, NAFTBW, specifies whether 
the body extends (NAFTBW = 1) or does not extend (NAFTBW = 0) behind the 
wing trailing edge. The last index, NAFTBT, specifies the same thing with 
respect to the tail. If NTAIL = 0 then NAFTBT should be input as zero. 
These last two indices are only used at supersonic Mach numbers. Thus, 
their values are immaterial at subsonic speeds. 

Item numbers 4 through 7 are self explanatory if the table in the 
preceeding section and figure 2 are referred to. Item number 6, which 
contains the tail data, is omitted from the input deck if there is no 
tail, NTAIL =0, in item number 3. If the nose is an ogive, NOSE =1 in 
item number 3, then the nose volume, VS of item number 7, is input as 
zero and the program calculates the volume. 

Item number 8 contains the lift-curve slopes, per degree, of the 
wing alone, CLAW, tail alone, CLAT, and nose, CLAN. The first two are 
determined by joining the exposed panels together. If experimental values 
of these quantities are known, they should be used. At subsonic Mach 
numbers CLAW and CLAT must be input. They can be obtained from, for 
example, reference 5. At supersonic speeds they can be input as zero and 
the program will determine them using Chart 8 of reference 2. If there 
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is no tail, NTAIL = 0, CLAT is input as zero. If the nose lift-curve 
slope is not known, it is input as zero and the program uses the slender- 
body value of 0.0349 per degree. 

The last three items of input, item numbers 9, 10, and 11, specify 
the ranges of angle of attack, wing incidence angle, and tail incidence 
angle for which calculations are to be performed. The first number on 
each card is the initial value of the angle, the second number is the 
final value, and the third number is the increment to be used in going 
from the initial to the final value. Calculations are performed for all 
combinations of these angles. If there is no tail, NTAIL = 0, zeros 
should be input for all three numbers in item 11. 

Sample cases.- Listings of the input data decks for two sample cases 
are presented in figure 3 and sketches of the two configurations are shown 
in figures 4 and 5. Sample case 1 is the example used in the computing 
form presented as Table I in reference 2. Sample case 2 is the configu- 
ration used in the tests of reference 6. 

Sample case 1 is a wing-body-tail combination. The nose is not an 
ogive, NOSE =0, so that the nose volume, VS, is input. The Mach number, 
FMACH, is 1.99. Thus, the lift-curve slopes, CLAW, CLAT, and CLAN, are 
input as zero and the program is allowed to calculate them. 

Sample case 2 is a canard-body-wing combination. The nose is an 
ogive, NOSE =1, so that the nose volume, VS, is input as zero and the 
program calculates this quantity. The Mach number for this case is sub- 
sonic, FMACH = 0.13, so that the wing (canard) and tail (wing) lift-curve 
slope values, CLAW and CLAT, are input. In this example the nose lift- 
curve slope value, CLAN, is also input. 

Description of Output 

The output produced by the R1307 computer program for sample case 1 
is shown in figure 6. The first output produced by the program, fig- 
ure 6(a), is a tabulation of most of the input data, Items 1, 2, and 4 
through 7 of figure 1. The next output, figure 6(b), lists quantities 
which are calculated by the program for the wing, tail, and nose and the 
lift-curve slopes which were either read in or calculated. The wing 
quantities which are tabulated are on the following page. 
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OUTPUT 

ALGEBRAIC 

NOTATION 

NOTATION 

EXPOSED AREA 

S W 

BETA*AR 


BETA*D/CR 

0d/c r 

AR PARAM 

0A W (1 + l/in0) (1 

TAPER RATIO 

X 

R/S 

(r / s) w 

SM*BETA 

m3 

CKWB 

^(B) 

CKBW 

^(W) 

SKBW 

*W(B) 

SKBW 

k B(W) 

(XBAR/CR) W (B) 

(x/ c r ) W (B) 

(XBAR/CR) B (W) 

* X//c r*B(W) 

(F - R)W/(S - R)W 

(f - r) w /(s - r), 


Except for the last quantity in the above list, the same quantities are 
tabulated for the tail. The listed nose quantities are: 


OUTPUT 

NOTATION 

BASE AREA 
CRN 


ALGEBRAIC 

NOTATION 





NOSE CENTER J 

OF PRESSURE ^N 


The definitions of the three lift-curve slopes are: 
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(CLA)W = 

(Q ) - 

da ^ L ; W 

_d_/ 

f L w \ 

da ' 


(CLA)T = 

(c ) = 

da '*'Ii T 

_d_ 1 

^ T N 

da 1 


(CLA)N = 

da ^ I/N 



da \ 

^ q oo S N/ 


The first output on figure 6(c) is a series of lift and moment curve 
slopes for the complete configuration without including wing-tail inter- 
ference. These six quantities appear in boxes 88 through 93 of Table I 
of reference 2. All of these slopes are evaluated at a = 6 W = 6,p = 0 
(A = DW = DT = 0 in the output notation) and the coefficients are formed 
using the input reference area S R and reference length £ R . For example 


and 


D (CL) /D (DW) = 


_d 

d6. 


W 


(C Lc } 


D (CM) /D (DW) = jf- (C mc ) 



The next quantities tabulated are lift and moment curve slopes for 
configuration components and the complete configuration including inter- 
ference of the wing vortices on the tail. These are also evaluated at 
a = 6^ = 6^ = 0. The column identified "BODY" is the nose component and 
that identified "WING-BODY" is the wing-body combination including the 
nose and wing-body interference. The column identified "TAIL-BODY MINUS 
NOSE" includes only the tail lift and that produced by tail-body inter- 
ference. The last column pertains to the complete configuration including 
wing-tail interference. The quantities tabulated in this block of output 
appear in boxes 124 through 127 of Table I of reference 2. 

The last output listed in figure 6(c) gives the lift and pitching- 
moment coefficients and the center of pressure of the complete configura- 
tion, both with and without wing-tail interference, as a function of angle 
of attack, A, wing incidence angle, DW, and tail incidence angle, DT. The 
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quantity denoted as LM - L/LR is the center of pressure measured from 
the center of moments and made dimensionless by the reference length. It 
is calculated as CMC/CLC. The ranges of these three angles for which 
calculations were to be performed were read in as items 9, 10, and 11 of 
the input data, see figure 1. The nine quantities in this table appear 
in boxes 94, 95, 96, 100, 104, 105, 120, 122, and 123 of Table I of 
reference 2. 

The last page of output, figure 6(d), tabulates for the same angle 
ranges the lift and pitching-moment coefficient components including those 
due to wing interference on the tail. The columns headed "NOSE" are nose- 
alone quantities due to angle of attack. The columns headed "WING A+DW" 
are wing quantities, including wing-body interference, due to angle of 
attack and v/ing incidence. The columns headed "TAIL A+DT" are tail quan- 
tities, including tail-body interference, due to angle of attack and tail 
incidence. The remaining two columns headed "TAIL WING INT" are the coeffi- 
cients produced by the tail due to wing vortex interference. The quantities 
tabulated on this page are not included in Table I of reference 2. 

Figure 7 contains the output for sample case 2. The format is iden- 
tical to that which has just been described. Since the configuration for 
sample case 2 is a canard-body-wing combination, "WING" in the output 
refers to the canard and "TAIL" refers to the wing. 


Program Listing 

The R1307 computer program consists of the main program, thirteen 
function subprograms, and four subroutine subprograms . Each source deck 
is identified in columns 73 through 80 by a four-character identification 
and a three-digit number sequencing the cards within that deck. The 
program listing is given on the following pages. The table below will 
act as a table of contents for the listing. 


PROGRAM 

MAIN 

Functions : 

APENB 

APENC 


IDENTIFICATION 

LN01 

LN02 

LN03 


PAGE NO. 
18 


21 

21 
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PROGRAM 

IDENTIFICATION 

PAGE NO 

CHRT8 

LN04 

21 

EQ14 

LN05 

22 

EQ19 

LN06 

22 

EQ21 

LN07 

22 

EQ24 

LN08 

22 

EQ24L 

LN09 

22 

EQ26 

LN10 

22 

EQ30 

LN11 

23 

EQ30L 

LN12 

23 

EQ31 

LN13 

23 

VOLOG 

LN14 

23 

Subroutines: 



CH1011 

LN15 

24 

CH1416 

LN16 

24 

CH56 

LN17 

26 

KFACT 

LN18 

26 
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ITEM 3 FORMAT (415) ; 1 card 
I NTAII4 NOSE INAFTBWINA 



Figure 1.- Input format for R1307 program. 
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SAMPLE CASE 

1 WING 

• BUOY * 

TAIL CONFIGURATION 

SAMPLE CALCULATION i 

OF TABLE I 

OF REPORT 1307 

this is no 
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TABLE III 

OF THAT 

REPORT 

1 o 

1 0 
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1.25 

45,0 

0,562 

3.19 

1.56 



0,0 

0,0 

0,0 



0,0 

15,0 

5,0 



0,0 

4,9 
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9,16 


l, 37b 


(a) Sample case 1. 
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•AMPLE CASE 2 canard 
REF, NASA TM X*643 
1110 
*13 39,3 406, 

1,03 5.03 0, 

2. 15,55 0, 

.«« 7,5 0, 

,0305 ,0418 ,0341 

0, 24, 4, 

15, 5, 

9 * 0 , 0 , 


CONFIGURATION 
BRADY, PAGE, KOENIG 

19,11 

B, 05 63,43 7,5 

25, 59, 29, 


• BODY • WING 

JAN, 1962 


3,5 


(b) Sample case 2. 


Figure 3.- Sample input decks for R1307 program. 
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Figure 4.- Sample case 1 configuration. 
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LIFT AND MOMENT CURVE SLOPES OF COMPLETE CONFIGURATION 
NO WING-TAIL INTERFERENCE 
PER DEGREE 

D ( CL ) /D ( A ) DICLI/DIDWI D(CL)/DIDT) D(CM)/D(A) DfCMI/OIDW) DtCM)/D(DT> 

0.080T7 0.04714 0.01581 -0.00606 -0.00053 -0.00714 
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Figure 6.- Concluded 
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Figure 7.- Output from R1307 program for sample case 
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Figure 7.- Continued 
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Figure 7.- Concluded. 


PART II - CRSFLW COMPUTER PROGRAM 


Introduction 

This computer program is based on the method presented in NASA TN 
D-6996 and TN D-7228 (refs. 3 and 4) to calculate the normal force, axial 
force, and pitching moment of a body alone or wing-body combination. The 
bodies must be slender but may have circular or noncircular cross sections, 
and the method is not restricted to small angles of attack. The method of 
calculation of normal-force and pitching-moment coefficients is based on 
the concept from reference 7 that the normal-force distribution over a 
body is made up of a potential term given by slender-body theory and a 
viscous crossflow term modified by Newtonian theory. Empirical information 
on crossflow drag coefficients as a function of Mach number and Reynolds 
number is incorporated into the program. This allows the procedure to be 
applied over a wide range of angles of attack, Mach numbers, and Reynolds 
numbers . 

The following sections present a description of the program, a 
description of the input, a description of the output, a program listing, 
and sample cases. A list of symbols from references 3 and 4 is included 
for reference. 

The program is written in FORTRAN IV for the IBM 360 series machines. 
No tapes, drums, or disks other than the standard input/ output units are 
required. The running time for a typical case on the IBM 360/67 is three 
seconds. The program will run on other machines such as the CDC 6600 with 
no modifications. 


A 

P 

A 

r 

a,b 


List of Symbols 
body base area (at x = £) 

planform area 

reference area 

semimajor and semiminor axes of elliptic cross section 

F 


axial-force coefficient, 
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c 


n 


C D 

C L 

C m 


"N 


n 


d 

D 

k 


i 

L 


^co 

r 

Re 

s 

V 

V oo 

W 

X 

X 


crossflow drag coefficient of circular cylinder section, 
F 


n 


% (A V d cy 


drag coefficient, 

< 3oo A r 


lift coefficient. 


lift 

^» A r 


pitching-moment coefficient about station x m from nose, 
pitching moment 

<2oo A r x 


normal-force coefficient, 


n 


^ A r 

local normal-force coefficient per unit length 


body cross-section diameter 
drag 

ratio of corner radius to body width for bodies of square cross 
section 

body length 

lift 

free-stream Mach number 


1 2 

free-stream dynamic pressure, pV TO 
body cross-section radius 

P V 00 d 

free-stream Reynolds number, — - — 

wing semispan 
body volume 
free-stream velocity 

body width 

reference length 

axial distance from body nose 
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X 


ac 


x 


X 


m 


a 


7 

e 


*1 

M* 

P 


cy 

LE 

Newt 

SB 

SF 

TE 

W 


distance from nose to aerodynamic force center 
distance from nose to centroid of body plariform area 
distance from nose to pitching-moment reference center 
angle of attack 

ratio of specific heats (taken as 1.4 for air) 
wing planform semiapex angle 

crossflow drag proportionality factor 
viscosity coefficient of air 
density of air 


Subscripts 

cylinder 
leading edge 
Newtonian theory 
s lender-body theory 
skin friction 
trailing edge 
wave or pressure 


Description of Program 

The CRSFLW computer program consists of a main program, eight func- 
tion subprograms, and two subroutine subprograms. All input and output 
takes place in the main program. The function and subroutine subprograms 
provide quantities from the curves and equations of references 3 and 4 and 
other specific services to the main program. 

Function subprogram CNSB computes the ratio (C n /C no ) for winged 
circular cross-section bodies and winged elliptic cross-section bodies 
with the major axis parallel to the crossflow velocity or normal to the 
crossflow using equations (13) , (14) , and (15) of reference 4, respectively. 

This same routine is used for winged-bodies with varying cross sections in 
which the above ratio changes with x-distance. 
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Function CNNT computes the ratio ( c n/ c n 0 ) Newt for the same wing-body 
configurations considered above using equations (16) , (17) , and (18) of 

reference 4. 

Function CNRSB computes (C n /C no ) gB for bodies alone with similar 
cross sections over their length. The result is identically 1.0 for 
circular bodies, and results for elliptic bodies with their major axis 
perpendicular to the crossflow velocity, or normal to it, are computed 
using equations (21) and (22) of reference 3, respectively. For bodies 
with square cross sections with rounded corners, (C n /C nQ ) gB = 1.19 at 
k = 0 (no corner radius) and (C n /C no ) gB = at ^ = (circular cross 

section) ; therefore, linear interpolation between these two end points 
is used for intermediate cases. 

Function CNRNT computes ( c n/ c no^Newt ^ or bodies a l° ne with similar 
cross section over their length. The result is identically 1.0 for cir- 
cular bodies, and results for elliptic bodies with their major axis per- 
pendicular to the crossflow, or normal to it, are computed using equa- 
tions (21) and (22) of reference 3, respectively. Bodies with square 
cross sections with rounded corners are considered using equation (23) 
of the same reference. 

Function CAW computes the wave or pressure contribution to the axial 
force, Ca w , for various nose shapes and body combinations at M > 1. 
Forward facing conical -nosed bodies are considered using equation (10) 
of reference 3. C^ w for tangent ogive noses and for Newtonian 
minimum drag noses are obtained from correlation curves in figure 6 of 
the above reference. For circular bodies with flat noses, it is assumed 
that C£ w is equal to the stagnation pressure coefficient and figure 7 
of reference 3 is used. for conical-nosed bodies of elliptic cross 

section is computed using equation (28) of the same reference. 

Function CDN computes the crossflow drag coefficient, C^, of a cir- 
cular body as a function of crossflow Mach number and crossflow Reynolds 
number from figures 1, 2, and 3 of reference 3. 

Function ETA computes the ratio of the crossflow drag coefficient 
for a finite length cylinder to that for an infinite length cylinder using 
figure 4 in reference 3. 
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Function FUN calculates the integrand in equations (7) and (8) of 
reference 4 for C N and C m for the cases with variable cross-sectional 
shapes over the length of the body or wing-body combinations. 

Subroutine SIMP is a Simpson's Rule integration package used to 
evaluate the integrals for C N and C m when the cross-sectional shape is 
variable over the length of the body. 

Subroutine CEL 2 computes the complete elliptic integral of the second 
kind, E. This is used in equation (10) of reference 4 to compute a modi- 
fication factor, X, defined as the ratio of the lift of the triangular 
wing alone by linearized theory to that by slender-body theory. 

Description of Input 

Variable definitions.- The format of the input cards for the CRSFLW 
program is shown in figure 1. The variable names are shown as well as 
the card columns in which the value is punched and the format in which it 
is punched. The remainder of this section consists of a table listing 
the program input variable along with the appropriate algebraic symbol, 
and the variable definition. A discussion of the preparation of the input 
is presented following the table. The algebraic notation used to define 
the configuration and described in the input table is shown in figure 2. 
All input length and area quantities are dimensional and should have 
consistent units. Also shown in figure 2 are sketches of a body alone 
with the positive sense of the forces and moments illustrated. 


PROGRAM 

ALGEBRAIC 


NOTATION 

NOTATION 

DEFINITION 

Item 1 

LHEAD 
Item 2 


80 columns of alphanumeric information 
used for identification. 

RL 

i 

Length of the body. 

XM 

x m 

Position about which the pitching moments 
are to be taken. 

AB 

A b 

Base area. 

AR 

A r 

Reference area. 
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PROGRAM 

NOTATION 

ALGEBRAIC 

NOTATION 

DEFINITION 

AP 

A P 

Planform area. 

V 

V 

Volume of body. 

XC 

x c 

x-coordinate of body centroid 

XX 

X 

Reference length. 

Item 3 
NSHP 


Body shape index: 


NSHP = 1 Circular body. 

= 2 Elliptical body with major 
axis horizontal. 

■ 3 Elliptical body with major 
axis vertical. 

= 4 Square body with rounded 
corners . 

Body shape similar over entire 
length. 

Body shape varies with length. 
(This option used only when 
NSHP = 2 or 3 . ) 

NWING NWING - 0 No wing present. 

/ 0 Wing present. 

NOSE NOSE ■ 1 Conical nose. 

= 2 Ogive nose. 

« 3 Newtonian minimum drag nose. 

= 4 Flat nose. 

NOTE : If NSHP - 4, then NOSE - 4 and NWING =0 is the only acceptable 

combination of variables. 

Item 4 

NXS Number of x-stations along the body where 

body shape information is input. 

2 < NXS < 50 

Item 5 

X(I) ,1=1, NXS x Values of x at each station at which 

body information is to be input. (If 
there are more than 8, continue on 
following cards.) X(l) = 0.0 and 
X (NXS) = RL. 



45 


PROGRAM ALGEBRAIC 


NOTATION 

NOTATION 

DEFINITION 

Item 6 


Item 6 consists of parts (a) , (b) , 

or (c) depending on body shape. 

(a) 

Circular Body 
(NSHP = 1) 




R(I) ,1=1, NXS 

r(x) 

The radius of the body at each x- 
station. 



or 


(b) 

Elliptical Body 
(NSHP = 2 or 3) 




(1) A ( I) , 1=1 , 
NXS 

a (x) 

The value of the semimajor axis at 
x-station; use as many cards as 
necessary, 8 values per card. 


(2) B ( I) ,1=1, 
NXS 

b (x) 

The value of the semiminor axis at 
each x-station „ 



or 


(c) 

Square Body 
(NSHP = 4) 




(1) RK 

k 

Ratio of the corner radius to the body 
width. Only one value allowed for 
entire body. (0 < k 0.5) 


(2) SSQ(I), 
1=1, NXS 

w 

Length of square side at each x- 
station. 

Item 7 


Item 7 consists of one card plus 
parts (a) or (b) depending on type 
of body. Omit Item 7 for blunt-nosed 
body. (NOSE « 4) 

SNOSE 

*n 

Length of nose. 

(a) 

Circular Body 
(NSHP = 1) 




RNOSE 

r 

Radius of body at base of nose. 



or 


(b) 

Elliptical Body 
(NSHP = 2 or 3) 

* 



ANOSE 

a 

Semimajor axis at base of nose. 


BNOSE 

b 

Semiminor axis at base of nose. 
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PROGRAM ALGEBRAIC 


NOTATION 

NOTATION 

DEFINITION 

Item 8 


Omit Item 8 if no wing is present. 
(NWING = 0) 

XLE 

*LE 

x-coordinate of the wing leading edge 
at the wing-body juncture. 

XLET 

*“tip 

x-coordinate of leading edge of wing 
tip ' * LB ti P ^ XLE 

XTET 


x-coordinate of trailing edge of wing 
tip- x > 

TE tip ^ E tip 

XTE 

X TE 

x-coordinate of wing trailing edge at 
the wing-body juncture, x ^ x 

TE TE tip 

SSPAN 

s 

Wing semispan. 

EPS 

€ 

Wing planform semiapex angle, radians 

Item 9 


Item 9 consists of (a) or (b) 
depending on body shape. Omit if 
no wing is present. (NWING = 0) 

(a) Circular Body 
(NSHP = 1) 



RLE 

r 

Radius of the body at the juncture of 


the wing leading edge and the body. 


or 


(b) Elliptical Body 
(NSHP = 2 or 3) 


ALE 

a 

Semimajor axis at juncture of wing 
leading edge and body. 

BLE 

b 

Semiminor axis at juncture of wing 
leading edge and body. 

Item 10 



ALFI 

a i 

Initial value of angle of attack, 
degrees. 

ALFF 

a f 

Final value of angle of attack, 
degrees. 

DALF 

Aa 

Increment in angle of attack. 
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PROGRAM 

ALGEBRAIC 


NOTATION 

NOTATION 

DEFINITION 

CASF 

C A 

SF 

Axial-force coefficient due to skin 
friction. 

GAM 

y 

Ratio of specific heats, typically 
1.4 for air. 

RE 

Re 

Reynolds number based on body diameter 
and free-stream properties. 

FMACH 

M 

00 

Mach number. 

. . . . . _ ^ o 


NOTE: The angle-of-attack range for bodies alone is 0 <_ a <. 180°. 


The angle-of-attack range for wing-body combinations is 0 < a 90°. 

Input preparation.- A discussion of the input variables is presented 
in this section as an aid in the preparation of the data deck. 

Item number 1 is a single card containing any information which the 
user wishes to use to identify the run. This information is reproduced in 
the output exactly as it is punched on this card. 

Item 2 is geometry and reference information and all the variables are 
explained in the previous table or in figure 2. 

Item 3 contains four indices controlling the type of body to be 
considered. The first, NSHP, specifies the body cross-sectional shape to 
be circular (NSHP = 1) , elliptical (NSHP = 2 or 3) , or square with rounded 
corners (NSHP = 4) . The second index, NVAR, specifies whether the body 
shape is similar (NVAR = 1) or not (NVAR = 2) over the length. This index 
is concerned only with elliptical cross sections (NSHP = 2 or 3) . The 
third index, NWING, specifies whether a wing is (NWING > 0) or is not 
(NWING = 0) present. The last index, NOSE, specifies the nose shape to 
be a cone (NOSE = 1) , an ogive (NOSE = 2) , a Newtonian minimum drag shape 
(NOSE = 3) , or a flat nose (NOSE = 4) . 

Item 4 is a single index, NXS, specifying the number of stations along 
the body at which shape information is to be input. This number must be 
equal to or less than 50. 

Item 5 contains the x-stations at which the body information is to be 
input. The first value should be the tip of the nose and the last should 
be the base of the body. There are eight values per card in increasing 
order of x distance up to the total of NXS stations. 
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Item 6 contains the variables defining the body shape at the x-stations 
in Item 5. For a circular cross section (NSHP = 1), Item 6 is the radius 
of each station. For an elliptic station, (NSHP = 2 or 3) , Item 6 is made 
up of two cards. Part (1) contains lengths of the semimajor axis at the 
x-stations, and part (2) contains the lengths of the semiminor axis at the 
same stations. If the body is square with rounded corners, the first part 
of Item 6 is a single card containing the ratio of the corner radius to 
the body width. Part (2) of Item 6 is the length of the side at each 
x-station. 

Item 7 is made up of two cards, the first of which contains the length 
of the nose, SNOSE. If the body cross section is circular, (NSHP = 1) , 
the second card of Item 7 contains the radius of the nose at the base, 
(RNOSE) . If the body cross section is elliptic (NSHP = 2 or 3) , the second 
card of Item 7 contains the length of the semimajor axis, ANOSE, and the 
length of the semiminor axis, BNOSE . Item 7 is omitted if NOSE = 4. 

Item 8 contains the variables describing the geometry of the wing. If 
no wing is present (NWING = 0) , Item 8 is omitted. The six variables in 
Item 8 are self explanatory. 

The variables in Item 9 are a function of the body shape in the 
vicinity of the wing. If no wing is present, Item 9 is omitted. For a 
circular body (NSHP = 1) , Item 9 contains the radius of the body (RLE) 
at the juncture of the wing leading edge and the body. For an elliptic 
body (NSHP = 2 or 3) , Item 9 contains the lengths of the semimajor and 
semiminor axes, ALE and BLE, respectively, at the juncture of the wing 
leading edge and the body. 

Item 10 is the last card making up a particular run. The first three 
variables are the initial angle of attack, ALFI, the final angle of attack, 
ALFF, and the increment in angle of attack DALF. The next quantity is the 
axial-force coefficient due to skin friction, CASF. GAM is the ratio of 
specific heats. RE is the free-stream Reynolds number based on body 
diameter and FMACH is the free-stream Mach number. The angle-of-attack 
range for bodies alone is 0° to 180° and that for wing-body combinations 
is 0° to 90°. 

Input decks may be stacked for multiple runs. A second case starting 
with Item 1 can be placed directly after Item 10. 
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Sample cases. - Listings of the input decks for five sample cases are 
presented in figure 3 and sketches of the configurations are shown in fig- 
ure 4. Sample cases 1 and 2 are body-alone configurations taken from 
figure 9 of reference 3. 

Sample case 1 is body number 2 of that reference and is a flat-nosed 
cylinder. Sample case 2 is body number 5 of the same reference and is the 
same cylindrical body with a conical nose attached. Sample cases 3, 4, 
and 5 are wing-body configurations taken from figure 6 of reference 4. 
These latter configurations all have the same body length, base area, and 
aspect ratio. 


Description of Output 

The output produced by the CRSFLW computer program for sample case 4 
is shown in figure 5. The first page of output is a summary of input 
quantities. Various notes are printed describing the specified components 
of the configuration. The last items printed on this first page are the 
flow conditions: Re, y and M^. 

The next page is headed by the identification information on the first 
card of the input deck. Following this are the calculated results printed 
on two lines. As indicated by the heading, the first line contains a, 

C N , c A , C m , x ac , C L , C D , and L/D. The second line contains components of 
Cjj, C A , and C m . CN1 and CMl are the potential portions of the normal 
force and pitching moments, and CN2 and CM2 are the viscous crossflow 
portions of the lift and pitching moments. The axial-force components 
are defined in the following table. The last variable is the crossflow 
drag coefficient, C^. The output variables are defined as follows. 


OUTPUT 

NOTATION 

ALBEGRAIC 

NOTATION 

ALPHA 

a, 

degrees 

CN 

C N 

*1$ 

II 

CA 

C A 

II 

CM 

C m 

M 

qA r X 
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OUTPUT 


ALGEBRAIC 


NOTATION 

NOTATION 

XAC 

/ X 
( m 
X ac \ X 

CL 

C T = C. T cos 
L N 

CD 

C D = °N Sin 

L/D 

V C D 

CDC 

C <3 n 


X 

** c 

~A 

V 


I. For nonvarying cross-sectional shape and no wing present. 


OUTPUT 

ALGEBRAIC 

NOTATION 

NOTATION 

CNl 

A, 

sin 2a cos j | 

CN2 

’' c d n si >> 2 “ (J?-') 

CMl 

rp - 

CM2 

r a_ /x - x 

[^n ( x 

CAl 

C A W 

CA2 

c Asf 

CA3 

c *b 


O ' SB 


in 2a cos 


sV- 


] si 

-) 




2 JV C ”o) 


SB 
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IX. For 3 varying cross-sectional shape or a body with lifting 
surfaces. 


OUTPUT 

NOTATION 


CNl 


ALGEBRAIC 

NOTATION 


(ib 

\A r 


sin 2a cos 




dx 


SB 


CN2 


CMl 


2r)C ( j n sin z a 



r dx 



CM2 


2r|Cd n sin 2 a 

A r X 



- x) dx 


Program Listing 

The CRSFLW computer program consists of the main program, eight 
function subprograms and two subroutine subprograms. Each source deck 
is identified in columns 73 through 80 by a four-character identification 
and a three-digit number sequencing the cards within that deck. The 
program listing is given on the following pages. The table below will 
act as a table of contents for the listing. 


PROGRAM 

IDENTIFICATION 

PAGE 

MAIN 

CF01 

54 

CNSB 

CF02 

56 

CNNT 

CF03 

57 

CNRSB 

CF04 

58 

CNRNT 

CF05 

58 

CAW 

CF06 

58 

CDN 

CF07 

59 

ETA 

CF08 

60 
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PROGRAM 

IDENTIFICATION 

PAGE NO. 

FUN 

CF09 

60 

SIMP 

CF10 

60 

CEL2 

CF11 

61 
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Figure 1.- Input format for CRSFLW program 












ITEM 6c FORMAT (8F10.5) 
(NSHP = 4) ,1 
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ITEM 10 FORMAT (8F10.5) 

I 1 ALFI | u ALFF I 21 DALF I 31 CASF 1 1,1 GAM 1 51 RE 1 61 FMACH 








Sample 

CASE 1 

NASA T N 

0-6996 

80DY NO, 2 

FLAT-NOSE CYLINDER 


30,46 

1.5,24 

11.40 

11,40 

116,13 

347,32 

15,24 

3,81 

1 

2 

1 0 

4 






0,0 

30,46 







1,905 

1,905 







0,0 

180,0 

10,0 

0,0 

1.4 

125000, 

2,66 


sample 

CASE 2 

NASA TN 

0-6996 

600 Y NO, 5 

CONE- 

CYLINDER 


41.91 

20,96 

11,40 

11,40 

137,90 

390,74 

23,66 

3,81 

i 

3 

1 0 

1 






O 

* 

o 

11,43 

41.91 






0,0 

1.905 

1,905 






11,43 








1,905 








0,0 

180,0 

10,0 

0,0 

1.4 

125000, 

2,66 



SAMPLE CASE 3 NASA TN 0-7228 MING-CONICAL BODY 1.5 ASPECT RATIO 


25,16 


25,16 

36,96 

36,96 

237,95 

310,22 

18,68 

25,16 

1 

2 

1 

1 

1 






0,0 


25,16 







0,0 


3,43 







25,16 









3,43 









0,0 


25,18 

25,18 

25,16 

9,45 

,359 



0,0 









0,0 


90,0 

5,0 

0,0 

1.4 

2180000, 

1.97 


SAMPLE 

CASE 4 

NASA TN 0-7228 WING-ELLIPTICAL 

BODY , MAJOR 

AXIS 

HORIZONTAL 

25,18 


25,16 

36,96 

36,96 

237,95 

310,22 

16,66 

25,16 

2 

1 

1 

1 






2 









0,0 


25,18 







0.0 


5,94 







0,0 


1,96 







25,18 









5,94 


1,98 







0,0 


25,16 

25,16 

25,18 

9,45 

.359 



0,0 


0,0 







0 i 


90, 

5, 

0,0 

1.4 

2180000, 1 

.97 


sample i 

CASE 5 

NASA TN 0*7228 

> MING-ELLIPTICAL 

BODY , MAJOR 

AXIS 

vertical 

25,18 


25,16 

36,96 

36,96 

237.95 

310,22 

16,86 

25,18 

3 

1 

1 

l 






2 









0,0 


25,16 







0,0 


5,94 







0,0 


1,98 







25,18 









5,94 


1,96 







o.o 


25,18 

25,18 

25,16 

9,45 

,359 



0,0 


0,0 







0,0 


90,0 

5,0 

0,0 

1.4 

2180000, 

1.97 



Figure 3.- Sample input decks for CRSFLW program. 
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Sample Case 2 


!<> 


l - 41.91 
d - 3.81 



25.18 cm 
5.94 
1.98 
9.47 
20.61° 

£ m 25.18 cm 
a - 5,94 
b - 1.98 
s - 9.47 
€ - 20.61° 

Figure 4.- Configurations used in sample cases 
for CRSFEW program. 
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ELLIPTICAL 8U0V ( N$MP ■ 2) BODY SHAPE S|*ILA« OVER ItNUTh 
TABULATION UP A AND 0 
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Figure 5.- Output from CRSFLW program. 
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Figure 5.- Concluded 


PART III - SUBSON COMPUTER PROGRAM 


Introduction 

This computer program predicts the static longitudinal aerodynamic 
characteristics of wing-body-tail combinations at subsonic speeds. It is 
an extension of the method of reference 2 to angles of attack for which 
symmetrical body vortices are shed from the nose of the configuration and 
leading-edge and side-edge separation vortices are shed from the wing and 
tail. The body is limited to circular cross-section shapes. 

The program is written in FORTRAN IV for the IBM 360 series machines. 
No tapes, drums, or disks other than the standard input/output units are 
required. Minor changes are required to run the program on other machines 
such as the CDC 6600. Typical running time on the IBM 360/67 for a wing- 
body— tail configuration is the order of one to two minutes. Actual time 
is dependent on the type of vortex lattice used to represent the lifting 
surfaces and whether or not a trim condition is calculated. Some specific 
running times are noted in the discussion of the sample cases. 

The following sections present a description of the program, a 
description of the input, a description of the output, a program listing, 
and sample cases. The algebraic notation used in this section is the same 
as that used in reference 1. A list of symbols from reference 1 is 
included for reference. 


List of Symbols 

JR aspect ratio 

a local body radius 

b semispan 

C, axial-force coefficient 

A 

C dc crossflow drag coefficient 

C D ^ induced drag coefficient 

C L lift coefficient, 

. . M 

C pitching-moment coefficient, — t-t 

m qt> i 
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N 

qS 


normal-force coefficient, 
local chord 

section-lift coefficient 
section leading-edge suction coefficient 

t 

vortex-lift ratio, figure 4 

lift force 
reference length 

pitching moment about center of moments, or free-stream Mach 
number 

normal force 

free-stream dynamic pressure 

body radius 

radius of base of nose 

reference area 

perturbation velocities along x,y,z directions, respectively 
free-stream velocity 

configuration coordinates with origin at body nose, figure 2 
x location of center of moments 

x position for onset of separation from body nose 

center-of-pressure location 
z location of center of moments 

body angle of attack 

Vl - M 2 

right body-vortex strength, positive counterclockwise when 
viewed from rear of configuration 

n'th separation-vortex strength on right wing panel 
trail ing-vortex strength on right wing panel 



6 

0 N 

A 

P 

a 

0 

A 

avg 

B 

B (T) 
B(W) 
c 

CP 

e 

HL 

LE 

N 

P 

root 

SE 

T(B) 

TE 

t 

tip 

v 

W(B) 

w 


lifting-surface deflection angle, positive trailing edge down 
nose angle, degrees 

sweep angle 
density 

complex vortex position, y + iz 
dihedral angle 

Subscripts 

afterbody 

average 

body 

body in presence of tail 
body in presence of wing 
canard 

center of pressure 
tail or empennage 
hinge line 
leading edge 
nose 

potential 
root chord 
side edge 

tail in presence of body 
trailing edge 
trailing vortex 
tip chord 
vortex 

wing in presence of body 
wing 
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Description of Method 

A brief description of the method is presented herein. The user 
should consult reference 1 for a complete description and details of the 
theoretical approach. 

An axisymmetric nose at some moderate angle of attack sheds a symme- 
tric pair of body vortices. These shed body vortices, whose strength and 
position are determined from data correlations, are tracked downstream 
past the wing using slender-body techniques in the crossflow plane. The 
vortex-induced velocities are computed at the wing control points and 
combined with the Beskin upwash induced by the body to obtain the total 
upwash induced on each wing panel. This, added to the free-stream contri- 
bution, results in a total local incidence angle distribution over the 
wing. 

The lifting surfaces (wing and horizontal tail or canard) are modeled 
by a vortex-lattice scheme (ref. 8) which has the capability to include 
velocity fields from external sources. The total upwash at the control 
points must be cancelled by the wing-circulation-induced velocity to 
satisfy the tangency boundary condition of the vortex-lattice method. 

The wing loading and trailing- vortex strength and position are obtained 
from this vortex-lattice calculation. The distribution of leading-edge 
suction and side-edge suction (if present) and their associated vortex 
positions and strengths are also obtained from the vortex-lattice calcu- 
lations. The leading-edge separation vortex lift is obtained from the 
suction distribution with the help of the Polhamus vortex-lift analogy 
(ref. 9) and correlation curves. 

The trajectories of the body vortices, the wing trailing vortex, and 
the wing leading-edge separation vortices are computed downstream past the 
afterbody and horizontal tail. These trajectories are computed in the 
crossflow plane considering mutual interference between the vortices and 
interference from their images in the body. The induced velocity field 
at the tail is computed, and the tail loading is obtained in a manner 
similar to that just described for the wing. The forces on the body due 
to the presence of the wing and tail are computed by the method of ref- 
erence 2. The free vortex-induced forces on the body are computed using 
the method of Sacks (ref. 10) . 
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The forces and moments on the entire configuration are obtained by 
summing the contributions of the various components. These forces are 
resolved into normal and axial force (excluding frictional drag) , and 
lift and induced drag. 

The subsonic prediction method includes an option to compute the trim 
conditions of a wing-body-tail configuration at some specified angle of 
attack. This is carried out by an iterative process in which the incidence 
of the tail or wing (canard) is varied until a zero pitching moment is 
achieved. 


Description of Program 

The SUBSON computer program consists of a main program and fifteen 
subroutines. The main program (SB01) accepts all the input, prints most 
of the output, and generally directs the flow of the calculation. The 
subroutines or groups of subroutines provide specific services to the main 
program during the calculation procedure. The following is a list of the 
subroutines and their general purpose. 

Subroutine LATTUS sets up the horseshoe vortex-lattice arrangement 
for the lifting surfaces. It locates the coordinates of the control 
points, calculates the influence coefficient matrix, and computes any 
geometry-related parameters connected with the lifting surfaces. 

Subroutine SHAPE does a table look-up for the body radius and slope 
and local lifting-surface semispan at any prescribed axial station. 

Subroutine BDYVTX uses tables derived from data correlations to look 
up the strength and position of the pair of symmetric vortices shed from 
the body nose. 

Subroutine CNVNZ computes the nose vortex-induced normal force and 
pitching moment on the nose of the body using the method of reference 10. 

Subroutine TRJTKY computes the trajectories of the free vortices past 
the configuration using the subroutines FCT, OUTP, HPCG, and SHAPE. 

Subroutine FCT computes the derivatives in the equations of motion 
for each free vortex. 

Subroutine HPCG is a predictor-corrector integration package which 
uses a Runge-Kutta starting procedure. 
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Subroutine OUTP stores the vortex positions in a table at specified 
intervals in x. When necessary, some diagnostic information on the vortex 
trajectories is available as optional output. 

Subroutine EXTV EL computes the vortex-induced velocity at wing or 
tail control points. 

Subroutine VTXLAT computes the strengths of the bound vortices on the 
lifting surface. It also computes the leading-edge suction distribution, 
the strength and spanwise position of the associated separation vortex, 
and the strength and position of the trailing vortex. This subroutine 
calls subroutines LQADl and INVERS. 

Subroutine LOADl uses the circulation distribution from the previous 
subroutine to compute the span-loading distribution and the forces and 
center of pressure on the lifting surface. 

Subroutine INVERS solves a system of linear simultaneous equations 
for the circulation strengths. 

Subroutine INFWW computes the influence function for a horseshoe 
vortex . 

Subroutine ZVTX determines the vertical position of the leading-edge 
separation vortex using a table look-up of correlated data for delta wings. 

Subroutine CNVTX computes the vortex- induced force and center of 
pressure on the afterbody using the method of reference 10. 

Description of Input 

Variable definitions.- The format of the input cards for the SUBSON 
program is shown in figure 1. In this figure the program variable name is 
shown as well as the card columns in which the value is punched and the 
format in which it is punched. The following is a table of the input 
variables along with the algebraic symbol where applicable. The input 
length and area quantities are dimensional and should have consistent 
dimensions. The variable is defined and its limits shown where necessary. 
The algebraic notation used in defining the configuration is shown in fig- 
ure 2. A discussion of the preparation of the input is presented in the 
section following the table. 
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PROGRAM 

NOTATION 

Item 1 

NHEAD 

NTBL 

NPRINT 


MPRINT 


NOSEV 

NTRIM 

NCWW 

MSWW 

NOTE : (NCWW 

NCAMW 

NPSIW 


ALGEBRAIC 

NOTATION DEFINITION 


MSWW) £ 100. 


Number of heading cards. 


Number of entries in table of body- 
coordinates . 

5 < NTBL < (96 - MSWW - MSWT) 

Output option: 

NPRINT < 0 Minimum output, final aero- 
dynamic characteristics 
only. 

= 0 Standard output. 

= 10 Optional additional output. 


Output option in trajectory calculation: 
MPRINT = 0 No additional output. 

> 0 Output vortex trajectories 

as calculated. (This option 
should be used as diagnostic 
only if program has prema- 
turely terminated execution 
during a previous trajectory 
calculation. ) 


NOSEV - 0 No nose separation vortex 
pair. 

« 1 Nose separation vortex pair 
included. 


NTRIM = 0 

> 0 

= 1 
= 2 


No trim calculation. 

Trim condition calculation, 
as follows. 

Wing incidence variable. 
Tail incidence variable. 


Number of chordwise vortices on wing. 
2 < NCWW < 10 

Number of spanwise vortices on wing. 

2 £ MSWW < 25 


NCAMW « 0 No wing camberline slopes to 
be input. 

■■ 1 Camberline slope at each wing 
control points to be input. 


NPSIW * 0 Unbroken wing leading edge 
and trailing edge. 

> 0 Input leading- and trail ing- 
edge sweep angles at each 
specified spanwise station. 
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PROGRAM 

ALGEBRAIC 


NOTATION 

NOTATION 

DEFINITION 

NSEFW 


Number of leading-edge separation vortices 
shed from wing. 1 NSEPW 2 

NCWT 


Number of chordwise vortices on tail. 
2 < NCWT <. 10 

MSWT 


Number of spanwise vortices on tail. 
2 < MSWT <_ 25 

NOTE: (NCWT X MSWT) < 100. 


NCAMT 


NCAMT - 0 NO tail cairiberline slopes to 
be input. 

> 0 Caniberline slopes at tail 

control points to be input. 

NPSIT 


NPSIT - 0 Straight tail leading edge 
and trailing edge. 

> 0 Input leading- and trailing- 
edge sweep angles at each 
specified tail spanwise 
station. 

NSEPT 


Number of leading-edge separation vortices 
shed from tail. 1 £ NSEPT £ 2 

NBODY 


Index controlling body upwash. 

NBODY < 0 No body upwash included in 
wing and tail interference 
calculation. 

— 0 Body upwash included. 

Item 2 

TITLE 


*Any alphabetic or numeric identification 
information. Number of cards equal to 
NHEAD. 

Item 3 

XBDY ( J) 


x-station at which body coordinates are 
defined* 

RBDY (J) 


Body radius at above stations. 

Item 4 

XM 

x m 

x-coordinate of moment center. 

ZM 

z m 

z-coordinate of moment center. 

EMACH 

M 

Mach number. 
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PROGRAM ALGEBRAIC 


NOTATION 

NOTATION 

DEFINITION 

REFS 

S 

Reference area. 

REFL 

i 

Reference length. 

THETAN 

0 N 

Nose angle, degrees. 


DXOUT 


x-increment in output table of free vortex 
trajectories (DXOUT < RAVGW) , typically. 

DXI 


Maximum integration interval for vortex 
trajectory calculation. (DXI DXOUT) 

Item 5 



XLEW 


x-coordinate of wing leading-edge inter- 
section with body. 

XTEW 


x-coordinate of wing trailing-edge inter- 
section with body. 

XHLW 

X. 

x-coordinate of wing hinge line at wing- 
body juncture. 

ZHLW 

*“V 

z-coordinate of wing hinge line at wing- 
body juncture. 

XCPW 


x-coordinate of alternate wing center of 
pressure location. 

BS2W 

b/2 

Wing semispan. 

RAVGW 


Average body radius at wing. 

YSEPW 


Spanwise location of 2 leading-edge 

separation vortex. If NSEPW = 1, 
YSEPW = b/2. 

Item 6 



PHIW 

4> 

Dihedral angle of wing, positive tip up, 
degrees. 

PSILEW(l) 


Leading-edge sweep angle at first wing 
station adjacent to body, degrees. 
Sweepback is positive. 

PSITEW(l) 


Trailing-edge sweep angle at first wing 
station adjacent to body, degrees. 


Sweepforward is negative. 
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PROGRAM 

NOTATION 

ALGEBRAIC 

NOTATION 

DEFINITION 

CMTEST 


Tolerance for C m in trim calculation if 
NTRIM / 0, CMTEST > 0.0. 

Item 7 



YW 

y 

y-coordinate of outboard side of lifting- 
surface panels on wing. Last value must 
equal b/2. MSWW values must be input. 

Item 8 


Delete if NPSIW = 0. 

PSILEW 

Vr 

Leading-edge sweep angle corresponding to 
values of y in Item 7, degrees. MSWW 
values must be input. 

Item 9 


Delete if NPSIW = 0. 

PSITEW 

***» 

Trailing-edge sweep angle corresponding 
to values of y in Item 7, degrees. 
MSWW values must be input. 

Item 10 



ALPHLW ( J) 

a 

w 

Tangent of local angle on wing due to 
camber and twist. (Values are input 
from leading edge to trailing edge, from 
root to tip.) MSWW x NCWW values are 
required. 

Item 11 



XLET 

** 

m 

CD 

x-coordinate of tail leading-edge inter- 
section with body. 

XTET 

xte . 

x-coordinate of tail trailing-edge inter- 
section with body. 

XHLT 


x-coordinate of tail hinge line. 

ZHLT 

ZHL e 

z-coordinate of tail hinge line. 

XCPT 

Xcp e 

x-coordinate of alternate tail center of 
pressure location. 

BS2T 

b/2 

Semi span of tail. 

RAVGT 

rav g e 

Average body radius at tail. 


YSEPT 
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Spanwise location of second leading-edge 
separation vortex, if NSEPT - 1, 

YSEPT - b/2. 


PROGRAM 

NOTATION 


ALGEBRAIC 

NOTATION 


DEFINITION 


Item 12 



PHIT 

0 

Tail dihedral angle, degrees. 

PSILET(l) 

^e 

Leading-edge sweep angle of tail at first 
station adjacent to body, degrees. 

PSITET (1) 

A TE e 

Trailing-edge sweep angle of tail at first 
station adjacent to body, degrees. 

Item 13 



YT 

y 

y-coordinate of outboard side of lifting- 
surface panels on tail. Last value must 
equal b/2 . 

Item 14 


Delete if MPSIT - 0. 

PSILET(J) 


Leading-edge sweep angle corresponding to 
values of y in Item 13, degrees. 

Item 15 


Delete if MPSIT = 0. 

PSITET (J) 


Trailing-edge sweep angle corresponding to 
values of y in Item 13, degrees. 

Item 16 



ALPHT(J) 

a e 

Tangent of local angle on tail due to 
camber and twist. (Values are input 
from leading edge to trailing edge, from 
root to tip.) MSWT X NCWT values are 
required. 

Item 17 



NDEX 


Index controlling next case of input. 

NDEX ■ 1 Execute program using variables 
on this card. 

* 0 Ignore this card and return to 
beginning for new case. 

ALPHAD 

a 

Angle of attack of configuration, degrees. 

ALPIW 

6 w 

Incidence angle of wing relative to body 
axis, degrees. 

AKVLWI 

K v 

LE w 

Fraction of leading-edge suction converted 
to lift in inboard wing region. 

(0 < KJus £. 1-0) 
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PROGRAM 

NOTATION 

ALGEBRAIC 

NOTATION 

DEFINITION 

AKVLW2 

K v 

LE w 

Fraction of leading-edge suction converted 
to lift in outboard wing region. 

1-0) 

AKVSW 

K v 

SE 

w 

Fraction of side-edge suction converted to 
lift on wing. (0 < K* SE < 1.0) 

ALPIT 

6 e 

Incidence angle of tail relative to body 
axis, degrees. 

AKVLTl 

0) 

w 

* 

Fraction of leading-edge suction converted 
to lift in inboard tail region. 

(0 < k£ le < 1.0) 

AKVLT2 

K* 

V LE e 

Fraction of leading-edge suction converted 
to lift in outboard tail region. 

(o < k Vle 1 1.0) 

AKVST 

SE e 

Fraction of side-edge suction converted to 
lift on tail. (0 < k£ se £1.0) 

WLIMIT 

W / V |max 

Limit on vortex- induced velocities at wing 
and tail control points. 

(0 < WLIMIT < 1.0) 


Input preparation.- A discussion of the input variables is presented 
in this section as an aid in the preparation of the input data deck. If 
a configuration has one set of lifting surfaces, this is denoted the wing 
regardless of its axial location on the body and data describing this 
lifting surface are input as wing data. If there are two lifting surfaces, 
the set nearest the nose is the wing and the aft set is the tail. For 
example, in a canard-body-wing configuration, the canard data are input 
as wing data and the wing data are input as tail data. In the following 
discussion, necessary geometric relations are illustrated in figure 2. 

Item number 1 of the input data deck (fig. 1) is a card containing 
indices specifying particular program options. NHEAD indicates the number 
of identification cards following in Item 2. NTBL is the number of entries 
in the table describing the body shape. Note again that only circular 
cross-sectional bodies are permitted. NPRINT is an index which determines 
the quantity of output obtained from the program. Typically, this number 
is zero. For diagnostic purposes, a provision for detailed output 
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information on lattice slopes, induced velocities, and circulation strengths 
is provided with NPRINT =10. An abbreviated output summary is obtained by 
setting NPRINT = -1. Examples of the output variatons are described further 
in the sample case output. MPRINT is a special index controlling the quan- 
tity of output obtained during vortex trajectory calculations. Typically, 
MPRINT = 0, but if it is greater than zero, the trajectory information is 
printed as it is computed. The additional trajectory information is useful 
only if the program fails to compute a trajectory for some particular case, 
and the additional output may give some clue as to why the calculation 
failed. This option should only be used after a computational problem has 
been discovered. 

The index NOSEV specifies whether or not there is (NOSEV =1) or is 
not (NOSEV = 0) a symmetrical pair of vortices shed from the nose of the 
configuration. If the nose angle (THETAN) is less than or equal to four 
degrees, NOSEV is automatically set equal to zero. The next index, NTRIM, 
determines whether or not a trim calculation is made. No trim calculation 
is made when NTRIM = 0, and when NTRIM > 0, a trim condition is calculated. 
If NTRIM = 1, the wing incidence is varied to achieve trim and if NTRIM = 2, 
the tail incidence is varied. 

The following five indices are associated with the wing. NCWW is 
the number of chordwise rows in the wing lattice. The only quantitative 
restrictions is that 2 < NCWW £ 10. MSWW is the number of spanwise 
columns in the wing lattice, and it must fall within the range 2 <_ MSWW 
25 • T* 1 * 2 total number of wing panels (NCWW x MSWW) must be less than or 

equal to 100. Some guidance in choosing a proper lattice arrangement for 
various shape wings is provided in reference 8. NCAMW specifies if wing 
camberline slopes are nonzero and must be input (NCAMW = 1) or are zero 
and need not be input (NCAMW = 0) . The index NPSIW identifies straight 
leading edges and trailing edges (NPSIW = 0) which requires input of only 
one value for leading— edge and trail ing-edge sweep angles, or broken 
leading and trailing edges (NPSIW = 1) which require input of the sweep 
angles at each spanwise station on the wing. NSEPW specifies the number 
°f leading-edge separation vortices shed from the wing, and it must be 
either 1 or 2. If NSEPW = 2, a special value of YSEPW is required in 
Item 5. 

The provision for multiple separation vortices is included to handle 
wings with breaks in leading-edge sweep. It has been observed that 
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leading-edge separation vortices are shed from the wing regions inboard 
and outboard of the break. A maximum of two leading-edge separation 
vortices are allowed. 

The next five indices in Item 1, NCWT, MSWT, NCAMT, NPSIT, and NSEPT 
are tail indices analogous to the previous five wing indices and subject 
to exactly the same restrictions. If no tail is present, all five must 
be zero. 

The last index in Item 1, NBODY, determines whether or not the upwash 
field around the body is included (NBODY > 0) in the wing and tail inter- 
ference calculation or not included (NBODY < 0) . This index is used to 
aid in determining the magnitude of the body- interference effect and 
generally should be set equal to zero. 

Item 2 is a group of NHEAD cards containing identification information 
which is printed at the beginning of the output. 

Item 3 is a group of NTBL cards describing the body shape. Each card 
contains an x-station, XBDY, and the corresponding body radius, RBDY. The 
cards should be in ascending order in x and there should be less than 
100 cards in this item. The program internally sets up its own table of 
coordinates which is stored in the XBDY and RBDY arrays and is limited to 
100 entries. A good rule of thumb to follow in inputting Item 3 is the 
following. 

NTBL < (96 - MSWW - MSWT) 

Some care is required when describing the body shape via XBDY and RBDY. 
Linear interpolation is used throughout; therefore, where the body shape 
is changing rapidly, more points are required. There should be a minimum 
of five entries in the nose region ahead of the wing and there must be 
entries at x-stations identically equal to XWLE and XWTE, and XTLE and 
XTTE if a tail is present. The last entry in the table must be greater 
than XWTE or XTTE, whichever is greater, by an amount not less than DXI. 

If the body is made up of a nose section followed by a cylindrical after- 
body, there should be two points on the cylinder very close together near 
the beginning of the cylinder. Points on a cylinder can be spaced large 
distances apart, but if the cylinder is followed by a section with 
changing radius, the last two points on the cylinder should be close 
together . 
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Item 4 consists of a single card containing XM and ZM, the coordinates 
of the center of moments, EMACH, the free-stream Mach number (0 <_ EMACH 
< 1.0) , and the reference area and reference length, REFS and REFL, 
respectively. THETAN is the nose semiapex angle in degrees (see fig. 2) . 

The final two variables are associated with the free vortex trajectory 
calculations. DXOUT is the approximate increment in x at which trajectory 
coordinates are stored for use in induced velocity calculations. A lower 
limit for this variable is about 0.5 percent of the overall length of the 
body because of storage limitations. Typically, a reasonable value for 
DXOUT is about one half the maximum radius of the body. DXI is the initial 
integration interval for the trajectory calculations. The integration 
package will cut the interval in half if necessary for reasonable accuracy, 
and this halving process can occur ten times before the program automati- 
cally terminates execution with an appropriate message. If DXI is made 
too large, the program will stop because of unacceptable accuracy, and if 
DXI is made too small, the running time will become large. Experience 
has indicated that a value of DXI between 2 and 5 percent of the body 
length will work for most cases. Under rare circumstances when two 
vortices get very close together or when a vortex gets very near the wing 
or body, a smaller value of DXI may be required. 

Item 5 contains geometric information for the wing. XLEW is the 
distance from the nose of the body to the intersection of the wing leading 
edge with the body. XTEW is the location of the trailing-edge intersection 
with the body. The wing hinge line at the wing-body juncture is located 
by the next two variables, XHLW and ZHLW. If an experimental center of 
pressure location is to be used for moment calculations, XCFW must contain 
the appropriate value. Otherwise, the program computes a center of 
pressure and XCFW must be identically zero. BS2W is the wing semispan 
measured from the centerline of the body. RAVGW is the average body radius 
in the vicinity of the wing. YSEFW is the y-station at which the wing is 
assumed broken for purposes of having two leading-edge separation vortices. 

If NSEFW = 1, then YSEFW must be equal to BS2W. If NSEPW = 2, YSEPW must 
be given some value greater than RAVGW and less than BS2W and the chosen 
value should coincide with one of the breaks in the lattice layout. That 
is, YSEPW will be equal to one of the values of YW to be described in 
Item 7. it is advised that there be at least three values of YW on either 
side of YSEFW to achieve reasonable accuracy in the separation vortex 
strength and position calculation. 
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Wing parameters are contined in Item 6. The first variable PHIW, is 
the dihedral angle in degrees for the entire wing. No breaks in dihedral 
are permitted. The second variable, PSIWLE(l), is the sweep angle, in 
degrees, of the leading edge at the wing-body juncture. If the leading 
edge has no breaks in sweep (NPSIW = 0) , this value is the only sweep 
angle associated with the leading edge. PSITEW(l) is the sweep angle, in 
degrees, of the trailing edge at the wing-body juncture. If there are no 
breaks in sweep, it must be the trailing-edge sweep angle. Remember that 
a swept forward trailing edge has a negative sweep angle. The last vari- 
able in this item, CMTEST, is the convergence tolerance on pitching moment 
for a trim calculation (NTRIM > 0) . A typical value for this quantity is 
about 1 percent of the magnitude of the untrimmed pitching moment. If it 
is made too small, the computer time required to converge to a trimmed 
solution can be very large. 

Item 7 is a list of the spanwise locations, YW, of the outboard side 
of each column of vortices. These quantities are dimensional spanwise 
distances measured from the body centerline. There are MSWW values input. 
The last value must be equal to the wing semispan, BS2W. These spanwise 
columns forming the wing lattice need not be equally spaced, but for 
convenience in preparing input, it is quite acceptable to use an equal 
spacing along the wing. 

Items 8 and 9 are optional and are included in the input deck only if 
there are breaks in sweep of the wing leading and trailing edges (NPSIW > 0) 
Item 8 includes the leading-edge sweep angle, in degrees, of each column of 
vortices from the wing-body juncture to the tip. If the wing sweep angle 
is continuously changing as in an ogee wing, the sweep angle at the center 
of each column of the lattice should be used. 

Item 9 contains the wing trailing-edge sweep angle corresponding to 
the leading-edge angles in Item 8. If the trailing edge is unbroken and 
has constant sweep, the values must still be input even though they are 
all the same. 

Item 10 is also optional and is included only if the wing is cambered 
or twisted (NCAMW > 0) . If such is the case, ALFHLW, the tangent of the 
local camber angle of each element of the lattice is input. There 

are MSWW cards corresponding to the number of spanwise columns forming 
the lattice, one card for each column. The camber angles on each card 
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run from leading edge to trailing edge with NCWW values per card. If there 
are more than eight chordwise rows, the ninth and tenth values follow on 
the next card. The ALPHLW value for the most forward area element in each 
column must start on a new card. 

Values of ALPHLW are obtained as follows. Consider the sketch in 
figure 3 which shows the cambered and twisted section of the lifting 
surface at some spanwise station. At point P, corresponding to a control 
point on the wing mean surface, a tangent to the wing mean surface is 
constructed, which makes an angle with the wing root chord. The 

positive sense of is shown in this figure. The input value required 

is ALPHLW = tan a^. Near the leading edge of the section shown in fig- 
ure 3, a ^ is negative. Item 10 completes the input description of the 
wing. 

If a tail or aft lifting surface is not present on the configuration 
(MSWT *= 0) , the next portion of input is Item 17 which specifies angle of 
attack and other nongeometric-related parameters. If a tail is present 
(MSWT > 1), Items 11 through 16 are required input. These items specifying 
the tail geometry are analogous to the equivalent wing parameters in 
Items 5 through 10 and the rules and restrictions regarding preparation 
tail input are the same as those described above for the wing. 

Item 17 is a group of cards, one card for each run, which specifies 
the variables which are considered changeable for a given geometric con- 
figuration. The first entry on the card is the index, NDEX, which is 
simply used to control the stacking of additional cases. NDEX = 1 on 
each card represents a new angle of attack or incidence angle condition. 

If NDEX = 0, the card is ignored and the program returns to read in a 
new case beginning with Item 1. Thus, a blank card is used to separate 
( 3 ; *-ff eren t cases. When NDEX ji 0, the next value on the card is the 
configuration angle of attack in degrees, ALPHAD, taken as the angle 
between the axis of the body and the free-stream velocity. The second 
quantity is the incidence angle of the wing root chord in degrees, ALPIW. 

Its sense is such that a positive incidence is a leading edge up condition. 
The next three variables are the K* factors which relate the actual 
realized vortex lift from the leading and side edges to that which is 
theoretically available. AKVLWl is the fraction of leading-edge separation 
vortex lift which is obtained on the inboard portion of the wing if 
NSEPW =2 or on the entire wing if NSEPW = 1. AKVLWl is a nuntoer between 
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zero and one and is generally geometry dependent. Its value can be ob- 
tained for sharp-edged delta wings from the correlation curves in figure 4. 
The source of figure 4 is described in detail in reference 1. The correla- 
tion curve can be used to get AKVLWl for any swept wing, but since figure 4 
was obtained for sharp-edged delta wings specifically, some judgement is 
necessary when other wings are considered. Instead of using the actual 
wing aspect ratio for nondelta wings, it is possible that some equivalent 
aspect ratio given by the delta wing expression 


would give a more reliable value. The factor is included as an input vari- 
able so that its effect can easily be examined by making a series of runs 
with AKVLWl varied between zero and one. 

If a wing is broken into two leading-edge vortex regions (NSEPW = 2) , 
then AKVLW2 is the K* LE factor which applies to the outboard portion of 
the wing. It is acceptable for AKVLWl and AKVLW2 to be equal. In the 
case of wings with breaks in sweep, the appropriate values can be obtained 
from figure 4 using an effective aspect ratio in each region as calculated 
by equation (1) . At the present stage in the development of program SUBSON 
there is no reliable method of choosing the correct factor. Many 

more data comparisons for double-delta wings or variable sweep wings 
should be made for this purpose. 

For unswept leading edges, the vortex lift from the leading edge is 
usually small. In this case, the full amount of vortex lift should be 
retained and both AKVLWl and AKVLW2 should be unity. For a wing with 
nonzero tip chord, the side-edge suction lift is generally very small 
compared to the potential lift except for very low-aspect-ratio wings. 
Comparisons with rectangular wing data indicate that the side-edge factor, 
AKVSW, should be unity at all times. When the tip chord is zero, such as 
on a delta wing, AKVSW should be identically zero. 

The next four variables, ALPIT, AKVLTl, AKVLT2 , and AKVST are the 
corresponding tail parameters. They fall under the same rules and guide- 
lines set up for the respective wing parameters. If no tail is present, 
ail four values should be set equal to zero. 
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The final quantity on this card is WLIMIT, the maximum allowable 
vortex- induced velocity nondimensionalized by free-stream velocity. The 
purpose of this variable is to limit the magnitude of the vortex-induced 
velocities on the wing or tail. In the course of program development, a 
canard-wing-body configuration developed a lift curve which appeared to 
exhibit a discontinuity around a = 8°. Close investigation showed that 
this occurred when the canard trail ing-vortex trajectory abruptly changed 
from passing beneath the wing to passing over the wing. At this point, 
the vortex-induced velocities changed character rapidly and, because of 
the close proximity of the vortex to the wing, the velocities were large. 
This created very large local angles of attack on the wing, and the vortex- 
lattice scheme predicted large changes in wing loading. This, of course, 
is an unrealistic situation because a true viscous vortex does not behave 
as a potential vortex and induce infinite velocities at its center. For 
this reason a limit was introduced which arbitrarily sets any vortex- 
induced velocity greater than WLIMIT, equal to WLIMIT. 

Generally, WLIMIT should be set equal to 1.0. if, in the process of 
running the program, unusual variations in the lift or pitching moment 
with angle of attack occur which can be attributed to unrealistic vortex- 
induced interference, WLIMIT can be used to limit the magnitude of the 
large induced velocities causing the problem. A value of WLIMIT = 0.1 
has been used in some specific examples to reduce the apparent discontinuity 
in the predicted lift and moment curves and resulted in good agreement with 
experiment. 

This discussion is not meant to suggest that an arbitrary velocity 
limit will cure the problems with the near flow fields of potential 
vortices. it is simply included to note that a simple, approximate fix 
is available. If WLIMIT = 0.0, the effect of the free vortices on the 
lifting surfaces is completely eliminated. 

The above discussion includes all the input required for a typical 
run. The sample cases in the following section cover the options available 
in the program. 

Sample cases.- Some sample cases are now presented to illustrate the 
preparation of input decks for various types of configurations . The airplane 
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configurations chosen for these examples are the canard-wing-body combina- 
tions of reference 6 shown in figure 5 and the wing-body-tail configuration 
from reference 11 shown in figure 6. 

in figure 7 (a) , the complete input deck for a canard-wing-body con- 
figuration is shown for sample case 1. The geometry corresponds to the 
sketch in figure 5 and this sample case considers the presence of both ^ 
lifting surfaces. This series runs is for four angles of attack (a = 4 , 
i2°, and 16°) with no canard deflection (6 W ■ 0°) and five angles of 
attack (a - 0°, 4°, 8°, 12°, and 16°) with positive canard deflection 
(6 - 10°) . Sample case 1 requires approximately 720 seconds on the 

IBM 360/67 computer; however, this time is much less than would be required 
if each run were made individually. A single run of this type requires 
approximately 100 seconds. 

The input deck for the second sample case is shown in figure 7 (b) . 

This deck is for the same configuration examined in the case above, but 
with the canard removed. In this case, the leading- and trailing-edge 
sweep angles are input at each spanwise station to illustrate the procedure 
for a wing with breaks in sweep. 

Sample case 3 is the wing-body-tail combination shown in figure 6. 

This input deck specifies a minimum amount of output and a trim calculation 
with the tail incidence variable. Only one angle of attack is specified 
because of the uncertainty in the amount of time required to converge on a 
trimmed solution. This particular run requires approximately 100 seconds 
on the IBM 360/67. Note that the vortex- induced velocity is limited to 
0.1 by the variable WLIMIT on the last card. This was necessary because 
of the large effect the wing shed vorticity had on the tail loading. The 
relative position between the trailing vortices and the tail was such that 
small changes in tail angle resulted in large nonlinear changes in tail 
loading which prevented convergence on a trim condition. 


Description of Output 

The output produced by the SUBSON computer program for sample case 1 
is shown in figure 8. The first page of output from the program, figure 8 (a 
is a tabulation of the input data in Items 1 and 2 of figure 1. The next 
page of output, figure 8(b), is a summary of the geometry of the configu- 
ration by component. The first quantities at the top of the page are the 
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first angle of attack to be considered, the Mach number, and beta. The 
next two items are the reference area and the reference length from Item 4 
of figure 1. If the wing and /or tail have alternate center-of-pressure 
locations input in Items 5 and 11 of figure 1, these are printed here. 

If a trim calculation is requested (NTRIM > 0) , the value of CMTEST speci- 
fied in Item 6 of figure 1 is printed here. These are followed by the 
geometry of the wing. 

The wing quantities which are tabulated are: 


OUTPUT 

NOTATION 

ALGEBRAIC 

NOTATION 

XLE 

*“root 

XTE 

"“root 

B/2 

b/2 

XHL 

*HL 

ZHL 

St 

C (ROOT) 

°root 

C (TIP) 

C tip 

Y (SEP) 

YSEPW (Item 5) 


The following block of data contains the wing dihedral angle in degrees 
and the leading-edge and trailing-edge sweep angles at the various y 
stations. The first entry denoted Y (RT) represents the wing-body juncture 
and the corresponding initial sweep angles of the wing leading and trailing 
edges. The following entries under the heading Y(WING) represent the y 
stations defining the spanwise lattice layout on the wing. The sweep 
angles at these stations are noted in the next two columns. 

The same quantities are tabulated for the tail surface if one is 
present. The following quantities are listed for the body. 
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OUTPUT 

NOTATION 

ALGEBRAIC 
■ NOTATION 

THETA 

9 (Item 4) 
N 

FINENESS 

*LE/ r N 

R(BASE) 

r N 

AVERAGE RADIUS WING 

rav 9 w 

AVERAGE RADIUS TAIL 

r avg e 

CENTER OF MOMENTS X 

x m 

CENTER OF MOMENTS Z 

z m 

DXI 

Ax *\ 

\ Item 4 

DXOUT 

- J 

X 

X 

R 

r 

S 

s or 
w e 

DR/DX 

dr/dx 


This concludes the general geometric description of the configuration. 
This information is output once at the beginning of each case. The follow- 
ing output is dependent on the information input on each card of Item 17; 
that is, the angle of attack, incidence angles, K* factors, and induced 

velocity limit. 

Figure 8(c) is the first page of output for each run within the series 
of runs making up sample case 1. The first line summarizes the information 
input in Item 17 as follows; 

OUTPUT ALGEBRAIC 

NOTATION JSOTATIOl^ 

ALPHA a 
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OUTPUT 

NOTATION 

ALGEBRAIC 

NOTATION 

INCIDENCE WING 

6 w 

INCIDENCE TAIL 

*e 

WING KVLE* 

K 

LE w 

WING KVSE* 

K 

SE 

w 

TAIL KVLE* 

LE e 

TAIL KVSE* 

K 

SE e 

W/V LIMIT 

(w / V) max 

The next block of output on this page is a summary of the strength 
and position of the right-hand vortex (if present) of the symmetrical pair 
of vortices shed from the nose of the body. The last entry, XS/RB, is the 
body-vortex separation location. 


OUTPUT 

ALGEBRAIC 

NOTOVTION 

jtotatioi^ 


r 0 

GAM/2 *PI*V*RB 

£ 

27rVr N 

Y/RB 

V r H 

Z/RB 

V r N 

XS/RB 

V r H 


The following block of output is the induced Beskin upwash at the 
wing control points (x,y,z) due to the presence of the body. The induced 
velocities, V/V(INF) and W/V(INF), expressed as a fraction of the free- 
stream velocity, are positive in the positive y- and z-directions, respeo 
tively. A summary of the nose-vortex position at the wing leading edge 
follows. The next block of data are the velocities induced at the same 
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wing control points by the vortex pair shed from the nose of the body. 

These velocities have the same positive sense as the body-induced upwash 
above. The final block of data in figure 8(c) is the total induced 
velocity at each control point. 

The next page of output, figure 8(d), contains the results from the 
lifting-surface calculations for the wing in a wing-alone coordinate system. 
Under FLOW CONDITIONS, the angle of attack is the incompressible angle of 
attack of the wing, including incidence. The next printed information is 
the lattice layout followed by the heading REFERENCE QUANTITIES. Under 
this heading, the actual exposed planform area of the wing and the average 
chord are listed. The aerodynamic coefficients on this page are based on 
these reference quantities. 

The following block of information contains wing geometry for the 
wing alone. If the Mach number is nonzero, the geometry is for the wing 
in the incompressible plane. 

The last half of this page contains the predicted aerodynamic char- 
acteristics of the wing in the presence of the body and other external 
interference velocity fields. Most of these quantities are self-explanatory 
and will not be described herein. It should be noted that in the coordinate 
system for the wing alone, x w is measured from the leading edge of the 
root chord, positive forward. The same positive direction is taken for 

all x direction coefficients. A few of the more important coefficients 
w 

are defined as follows: 


OUTPUT 

ALGEBRAIC 

NOTATION 

NOTATION 

CNP 

C„ 

n w(b) ,p 

CNV 

Cn w(b) ,v le 

CNVS 

Cn W(B),V se 

CC^ 

2b 

cc 

o 

CL*C/(2*B) 

CSUC*C/(2*B) 

o 

2b 
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OUTPUT 

NOTATION 

KVLE 

KVSE 

GAM/V2PI 


ALGEBRAIC 

NOTATION 



r 

2ttV 


The last line of figure 8(d) contains the predicted leading-edge and 
side-edge vortex-lift constants, K Vle and K VgE , respectively. Since these 
values include effects of external interference on the wing, they are 
useful only in cases of wings alone with no interference effects included. 
The last two items are the lateral position and strength of the trailing 
vortex and the leading-edge separation vortex at the wing trailing edge. 

Figure 8(e) is headed by a summary of the strengths and positions of 
the vortices shed from the configuration ahead of the wing trailing edge. 
The pairs of vortices are listed in the following order. Vortex 1 is the 
right-side body vortex shed from the nose. Vortex 2 is the trailing vortex 
shed from the wing. Vortex 3 is the leading-edge separation vortex shed 
from the wing. If more than one separation vortex is requested, vortex 3 
is the vortex associated with the inboard region and vortex 4 is shed from 
the outboard region. If a vortex is missing for any reason, all following 
vortices are moved up in the table. For example, if no vortices are shed 
by the nose, vortex 1 becomes the trailing vortex shed by the wing, and so 
on. The remainder of figure 8(e) is the induced velocities at the tail 
control points. These velocities are analogous to the induced velocities 
on the wing shown in figure 8(c) . 

Figure 8(f) contains calculated results for the tail surface. All 
the quantities on this page are analogous to those described for the wing 
in figure 8(d). The last entry on this figure is a summary of the 
strengths and positions of all the vortices in the field just aft of the 
tail trailing edge. The first group of vortices are the same as described 
in connection with figure 8(e). The second group of vortices are defined 
as follows. Vortex 4 is the trailing vortex corresponding to the potential 
lift on the tail. Vortex 5 is the leading-edge separation vortex shed from 
the tail, if multiple vortices are shed from the tail leading edge, this 
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vortex would be shed from the inboard tail region and vortex 6 would be 
shed from the outboard tail region. 

The next page of output, figure 8(g), is a summary page of the force 
coefficients, pitching-moment coefficients, and centers of pressure of 
each component of the configuration and of the total configuration. The 
coefficients for the individual components are described in Table I. The 
total configuration variables are defined as follows: 

OUTPUT ALGEBRAIC 


NOTATION 

NOTATION 

: 

CN 

c « JL 
'■'N qS 


CM 

c - -82- 
Sn qSi 


XCP 

X CP ■ x m 

-SEi 

CL 

c - JL 
U L qS 


CDI 

% ’ 

tan a 

CA 

C A " qS 


CDI/CL**2 

C D/ C L 



The last page of output for this run, figure 8(h), contains a summary 
of the trajectories of the shed vortices. At the top of the page the 
vortices are identified and their strengths listed. This is followed by 
blocks of output, one block for each x-station, describing the local 
crossflow geometry of the configuration and the position of the right-side 
vortices. The x stations of each block of results are approximately 
DXOUT apart. Notice that the trajectory calculation starts at the wing 
leading edge with a pair of body vortices. As the calculation moves 
downstream, other vortices are shed and added to the calculation. The 
trajectory calculation is carried downstream to a point aft of the tail 
trailing edge. For purposes of saving space in figure 8, only selected 
portions of the trajectory calculation is presented herein. The variables 
in each block are defined as follows: 
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OUTPUT 

NOTATION 


ALGEBRAIC 

NOTATION 


X 

DX 

A 

S 

RO 


DA/DX 

SIGMA (REAL) 
SIGMA (IMAG) 


x 

Ax 

a 

s or s 
w e 

r Q , transformed circle radius 
da/dx 

y 

z 


This completes the output for one card in Item 17 of the input deck. 
Additional runs will repeat the output of figures 8(c) through (h) . The 
above set of output obtained with N PRINT * 0 is a considerable amount of 
output for production runs; therefore, an optional set of output can be 
obtained by setting NPRINT = -1 . in this case, the complete output con- 
sists of figures 8(a), (b) , and (g) , with some shed vortex positions and 
strengths added. 

Some extra output over and above that shown in figure 8 can be 
obtained when NPRINT = 10. This additional output is useful only for 
diagnostic purposes and is not described herein. This output is labeled 
and the user should have no trouble interpreting the results. 


Program Listing 

The SUBSON computer program consists of the main program and fifteen 
subroutine subprograms. Each source deck is identified in columns 73 
through 80 by a four-character identification and a three-digit number 
sequencing the cards within that deck. The program listing is given on 
the following pages. The table below will act as a table of contents 
for the listing. 

PROGRAM IDENTIFICATION PAGE NO. 

MAIN SB01 97 

LATTUS SB02 104 


95 



PROGRAM 


IDENTIFICATION 


PAGE NO 


SHAPE 

SB03 

104 

BDYVTX 

SB04 

105 

CNVNZ 

SB05 

105 

TRJTPY 

SB06 

105 

FCT 

SB07 

106 

HPCG 

SB08 

106 

OUTP 

SB09 

108 

EXTV EL 

SB10 

109 

VTXLAT 

SBll 

109 

LQADl 

SB12 

110 

INVERS 

SB13 

113 

INFWW 

SB 14 

113 

ZVTX 

SB15 

114 

CNVTX 

SB16 

114 
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Mean surface 

Mean surface tangent at P 

Local chord projected 
parallel to root chord 



Wing root chord 


Figure 3.- Mean surface detail for a wing 
with camber, twist, and dihedral. 
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Figure 5.- Canard— body-wing configuration 
sample cases 1 and 2 m 



44.92 


_ 50.37 
56.16 — 


Figure 6.- Wing-body- tail configuration for 
sample case 3 . 
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SAMPLE CASE 1 

RfcF, NASA TM X »643 
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CANARD * BOOT * *InG CONFIGURATION 
BRAOV* PAGE * AND KOENIG JANUARY 1962 
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(a) Sample case 1. 

Figure 7.- Sample input decks for SUBSON program. 
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(b) Sample case 2. 
Figure 7.— Continued. 
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3 23 *1 0 1 2 6 10 0 0 1 9 8 0 010 

SAMPLE CASE 3 WING - BUOY • TAIL CONFIGURATION 

REF, NACA Rt* A52L15 0, KUENIG APRIL 1953 

TRJR CALCULATION 
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(c) Sample case 3. 
Figure 7.- Concluded. 
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CCARDl 
NHEAD= 2 
NTBL= 35 

nprint= 0 

MPRI NT* 0 
NOSEV= l 
NTRIM= 0 
NCWW= 4 
MSWW= 10 
NCAMW= 0 
NPSI W= 0 
NSFPW= l 
NCWT = R 
MSWT= 10 
NCAMT= 0 
NPS1T= 0 
NSFPT= l 
NBODY= 0 
6 END 


SAMPLE CASE 1 

REF. NASA TM X-643 


CANARD - BODY - WING 
RRADY, PAGE, AND KOENIG 


(a) Page 1. 

Figure 8.- Sample case from SUBSON program. 
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CONFIGURATION 
JANUARY 1962 



Finw conditions 


ALPHA MACH RETA 

16.00 0.0000 1.0000 


rONF IGUPATICK GEOMETRY - 


REFERENCE AREA « 446.000 

REFERENCE LENGTH - 19.110 


WING 


TAIL 


body 


XL F 

XTF 

A/2 

XHL 

7.5000 

15.5500 

5.0300 

11.0000 



SWEEP 

ANGLE 

DIHEDRAL 

Y|RT) 

LE 

TE 

0.000 

1.030 

63.430 

0.000 


Y ( WING ) 

LE 

TF 


1.430 

63.430 

0.000 


1.830 

63.430 

0.000 


2.230 

63.430 

0.000 


2.630 

63.430 

0.000 


3.030 

63.430 

0.000 


3.430 

63.430 

0.000 


3.830 

63.430 

0.000 


4.230 

63.430 

0.000 


4.630 

63.430 

0.000 


5.030 

63.430 

0.000 

XLE 

XTE 

0/2 

XHL 

29.0000 

54.0000 

15.5500 

39.3000 



SWEEP 

ANGLE 

DIHEDRAL 

VIRT) 

LE 

TE 

0.000 

2.000 

59.000 

-10.000 


V(TAIL) 

LE 

TE 


3.355 

59.000 

-10.000 


4.710 

59.000 

-10.000 


6.065 

59.000 

-10.000 


7.420 

59.000 

-10.000 


8.775 

59.000 

-10.000 


1 0. 120 

59.000 

-10.000 


11.485 

59.000 

-10.000 


12.840 

59.000 

-10.000 


t4. 195 

59.000 

-10.000 


15.550 

59.000 

-10.000 




AVERAI 

theta 

FINENESS 

R I BASE) 

WING 

8.000 

4.464 

0.840 

1.030 


ZHL 

0.0000 


ZHL 

0.0000 


TAIL 

2.000 


CIR00T1 CITIP1 YtSEPl 

8.0500 0.0000 5.0300 


C I ROOT ) CHIP) YISEPI 

25.0000 0.0000 15.5)00 


CENTER OF NOHENTS 

X Z 0X1 

39.300 0.000 0.500 


TABLE OF BODY COORDINATES 



X 

R 

S 

DR/OX 

I 

0. 00000 

0.02000 

0.02000 

0.14000 

2 

2.00000 

0.30000 

0. 30000 

0.12750 

3 

4.00000 

0.53000 

0.53000 

0.10500 

4 

6. 00000 

0. 72 0 00 

0. 72000 

0.08750 

5 

7.50000 

0.84000 

0.84000 

0. 07906 

6 

8.29982 

0.90249 

1.43000 

0.07656 

7 

9.09965 

0.96247 

1.83000 

0.07500 

8 

9.89947 

1.02246 

2.23000 

0.07500 

9 

10. 69930 

1.08245 

2.63000 

0.07500 

10 

11.49912 

1.14243 

3.03000 

0.07500 

11 

12.29095 

1.20242 

3.43000 

0.07500 

12 

13.09877 

1.26241 

3.83000 

0.07500 

13 

13. 89060 

1.32239 

4.23000 

0.07063 

14 

14.69842 

1.37540 

4.63000 

0.06563 

15 

15.49825 

1.42738 

5.03000 

0.06500 

16 

15.55100 

1.43081 

1.43081 

0.06500 

17 

16.00000 

1.46000 

1.46000 

0.06250 

18 

18. 00000 

1.58000 

1.58000 

0.05 500 

19 

20.00000 

1.68000 

1.60000 

0.04500 

20 

22.00000 

1.76000 

1.76000 

0.03750 

21 

24.00000 

1.83000 

1.03000 

0.03000 

22 

26.00000 

1.88000 

1.08000 

0.02250 

23 

20.00000 

1.92000 

1.92000 

0.01750 

24 

29.00000 

1.93500 

1.93500 

0.01500 

25 

31.25508 

1.96883 

3.35500 

0.01333 

26 

33.51018 

1.99510 

4.71000 

0.01083 

27 

35.76527 

2.01765 

6.06500 

0.00774 

28 

38.02036 

2.03000 

7.42000 

0.00274 

29 

40.27545 

2.03000 

8.77500 

0.00000 

30 

42.51390 

2.03000 

10.12000 

-0.00173 

31 

44.78564 

2.02214 

11.48500 

-0.00708 

32 

47.04073 

t. 99439 

12.84000 

-0.01653 

33 

49.29581 

1.94760 

14.19500 

-0.02459 

34 

51.55089 

1.88347 

15.55000 

-0.03739 

35 

54.00099 

1.76993 

1.76993 

-0.05817 

36 

56.00000 

1.63000 

1.63000 

-0.07250 

37 

58.00000 

1.49000 

1.40000 

-0.09000 

38 

60.00000 

1.27000 

1.27000 

-0.11250 

39 

62.00000 

1.03000 

1.03000 

-0.12500 

40 

64.00000 

0.77000 

0.77000 

-0.15000 

41 

66.00000 

0.43000 

0.43000 

-0.29000 

42 

67.00000 

0.02000 

0.02000 

-0.41000 


(b) Page 2. 
Figure 0.- Continued. 
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alpha 

16.00 

— INCIDENCE — 
M WING TAIL 

0.00 10.000 0.000 

W 

KVL p * 
0.530 

IMG 

KVSE* 

0.000 

TAIL 

KVLF* 

0.410 

KVSF* 

0.000 

W/V LIMIT 
1.000 

Nosr vopticfs - 

GAM/2*PI*V*«* Y/PR 

0.1013 0.484 1 

Z/Pft 

.435 

XS/RB 

0.000 




VriOCITY INDUCED 

AT SPECIFIED FIELO POINTS 

ON WING 

(BESKIN) 






X 

Y 

Z 

V/V( INF) 

W/V (INF) 

1 

9.334 

1.230 

0.289 

7.3907E-02 

1.4845E-01 

2 

11.247 

1.230 

-0.043 

-1.5989E-02 

2. 2914E-01 

"X 

13. 159 

1.230 

-0.375 

- 1 • 49 27E-01 

2 • 2207E-0 1 

4 

15.07? 

1.230 

-0.707 

-2.3183E-01 

1 • 3501E-01 

5 

9. 984 

1.630 

0.176 

2. 3220F-02 

1 . 0602F-0 l 

6 

11.697 

1.630 

-0.121 

-2.0400E-02 

1 • 3666E— 01 

7 

13.409 

1.630 

-0.41ft 

-7.7487E-02 

1 .41 OOF-Ol 

ft 

15.122 

1.630 

-0.716 

-l • 2603E-0 l 

1 • l 584E-01 

0 

10.634 

2.030 

0.064 

4.8533E-03 

7.7436E-02 

10 

12.147 

2.030 

-0.199 

-1.8259E-0? 

9.2183F-02 

ll 

13.659 

2.030 

-0.462 

-4.681 0E-02 

9.7566E-02 

12 

15.17? 

2.030 

-0.724 

-7. 4273E-0? 

9.0812E-0? 

13 

11.284 

2.430 

-0.049 

— 2 .4004E-03 

5.91 41E-02 

14 

12.597 

2.430 

-0.277 

— 1 • 5569E— 02 

6 • 7343 E— 02 

1 5 

13.909 

2.430 

-0.505 

-3.12 19E-02 

7.184 IE-02 

16 

15.222 

2.430 

-0.733 

-4. 7004E-02 

7.081 0E-02 

17 

11.934 

2 • ft 30 

-0.162 

- 5 • 4094E — 03 

4. 7053E-0? 

18 

13.046 

2.830 

-0.355 

-1.3267E-02 

5. 1996E-02 

19 

14.159 

2.830 

-0.549 

-2.2243E-02 

5.5218E-02 

20 

15.272 

2.830 

-0.74? 

-3 . 1525E-02 

5. 6004E-02 

’1 

12.584 

3.230 

-0.275 

-6.6405E-03 

3.8717E-02 

22 

13.496 

3.230 

-0.433 

- 1 • 1426E-02 

4. 1803E-02 

23 

14.409 

3.230 

-0.592 

-W6674E-02 

4.3960E-02 

24 

15.32? 

3.230 

-0.750 

-2.21 57E-02 

4.51 06E-02 

25 

13.233 

3.630 

-0.388 

-7. 0755E-03 

3. 2734E-02 

26 

13.946 

3.630 

-0.51? 

—9, 9608E-03 

3.4635E-02 

27 

14.659 

3.630 

-0.635 

-1.2992E-02 

3 • 5975E— 02 

?ft 

15.37? 

3.630 

-0.759 

-1.6 176E-02 

3.6982E-02 

29 

13.883 

4.030 

-0.501 

-7 • 1 396 E— 03 

2.8290E-02 

30 

14.396 

4.030 

-0.590 

— 8.746 IE— 03 

?.9?44E-02 

31 

14.909 

4.030 

-0.679 

-1 . 0431E-02 

3.0087E-02 

32 

15.422 

4.030 

-0.768 

-1.21 84E-02 

3.0814E-02 

33 

14.533 

4.430 

-0.614 

-6.9727E-03 

2 * 4690F-02 

34 

14.846 

4.430 

-0.668 

-7.7663E-03 

2.5.1 72E-02 

35 

15. 159 

4.430 

-0.722 

-8.58 1 8E-03 

2.5622E-02 

36 

15.472 

4.430 

-0.777 

-9.4201E-03 

2. 6045E-02 

37 

15. 183 

4.830 

-0.726 

-6 . 7264E-03 

2.1S58E-02 

3 P 

15.296 

4. 830 

-0.746 

—6 • 963 4E-03 

2.2005E-0? 

39 

15.409 

4.830 

-0.766 

—7 • 2029E-03 

2.21 50F-02 

40 

15.522 

4.830 

-0.785 

-7.4449E-03 

2.2293E-02 


SUMMARY OF VORTFX STRENGTHS AND POSITION AT X - 7.500 


VPRTEX GAMMA/2*P I *V Y l 

l 0.085127 0.406 1.206 


(c) Page 3. 

Figure 8.- Continued. 
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VFLOCITY INDUCED AT SPECIFIED FIELD JOINTS ON WING 



X 

Y 

2 

V/V ( INF 1 W/VUNF) 

1 

4.334 

1.230 

0.289 

1.3685E-0? -1.A366E-02 

2 

1 1 • 747 

1.230 

-O.OA3 

1.31A6E-03 -9.5A3AE-03 

3 

13.159 

1.230 

-0.375 

-1.2206E-03 -A.A13AE-03 

A 

15.072 

1.230 

-0.707 

-I.LA52E-03 -2 ■ 10A5E-03 

5 

9.98A 

1.630 

0.176 

8.63 1 7E-03 -5. A27AE-03 

6 

11.697 

1.630 

-0.121 

3.2A71E-03 -5.1 893E-03 

7 

13.A09 

1.630 

-0.A18 

7 . 990AE-OA -3.5265E-03 

8 

15.122 

1 .630 

-0.716 

-9.2118E-06 -2. 2802E-03 

9 

10.63A 

2.030 

0.06A 

5.6603E-03 -2.2AA3E-03 

10 

12.1A7 

2.030 

-0. 199 

3. 1833E-03 -2.7330E-03 

11 

13.659 

2.030 

-0.A62 

1.5972E-03 -2 . 38 59E-03 

12 

15.172 

2.030 

-0.72A 

7.7890E-0A - 1 . 91 6AE-03 

13 

U.28A 

2.A30 

-0.0A9 

3.8802E-03 -9.A032E-0A 

t A 

12.597 

2.A30 

-0.277 

7.678AE-03 -1.A122E-03 

15 

13.909 

2.A30 

-0.505 

1.7609E-03 - 1 . A827E-03 

16 

15.222 

2.A30 

-0.733 

l. 1751E-03 -1 • 3987E-03 

17 

11.93A 

2.830 

-0.162 

2.76 75E-03 -3.6009E-0A 

18 

13.0A6 

2.830 

-0.355 

2.1667E-03 -7.0150E-0A 

19 

1 A. 159 

2.830 

-0.5A9 

1.66AAE-03 -8.61 18E-0A 

20 

15.272 

2.830 

-0.7A2 

1.2976E-03 -9. 2238E-0A 

21 

12.56A 

3.230 

-0.275 

2.0382E-03 -8.9085E-05 

22 

13.A96 

3.230 

-0.A33 

1.7328E— 03 -3.11 38E-0A 

23 

IA.A09 

3.230 

-0.592 

1.A7A2E-03 -A. 5A8AE-0A 

2A 

15.322 

3.230 

-0.750 

1.2671E-03 -5.A8 75E-0A 

23 

13.233 

3.630 

-0.388 

1.5A91E-03 A. 0093E-05 

26 

13.9A6 

3.630 

-0.512 

1.1995E-03 -9. 5A39E-05 

27 

l A. 659 

3.630 

-0.635 

1.27A3E-03 — 1.9726E-0A 

28 

15.372 

3.630 

-0.759 

1.16 50E-03 —2. 70OIE-OA 

29 

13.883 

A. 030 

-0.501 

1 . 2079E-03 1.0062E-0A 

30 

1A.396 

A. 030 

-0.590 

1.1A55E-03 2.7879E-05 

31 

1 A. 909 

A. 030 

-0.679 

1 .09 17E-03 -3.531 5E-05 

32 

15.A22 

A. 030 

-0.768 

1.03 78E-03 — 9.068AE-05 

33 

1 A. 533 

A.A30 

-0.61A 

9.7275E-0A 1.3155E-0A 

3A 

1A.6A6 

A.A30 

-0.668 

9. 5290E-0A 9.7355E-05 

35 

15.159 

A. A30 

-0.722 

9.3 188E-0A 6. 51 37E-05 

36 

15.A72 

A.A30 

-0.777 

9. 09A0E-0A 3, A76AE-05 

37 

15.183 

A. 830 

-0.726 

8.0665E-0A 1.A557E-0A 

36 

15.296 

A. 830 

-0.7A6 

8. 01A2E-0A 1.3565E-0A 

39 

15.A09 

A. 830 

-0.706 

7.960AE-OA 1.7591E-0A 

AO 

15.522 

A. 830 

-0.785 

7.9090E-0A 1 • 1633E-0A 



TOTAL INDUCED 

VELOCITY 

AT SPECIFIED FIELD POINTS ON MING 


X 

Y 

Z 

V/VI INF! 

W/V 1 INF 1 

1 

9.33A 

1.230 

0.289 

8.7592E-02 

1.3409E-01 

2 

11.247 

1.230 

-0.043 

-1.4674E-02 

2. 1960E— 01 

3 

13.159 

1.230 

-0.375 

-l • 5049E-01 

2.1766E-01 

A 

15.077 

1.230 

-0.707 

-2.3297E-01 

1.3290E-01 

5 

9.904 

1.630 

0. 176 

3. 1852E-02 

1 . 0060E-01 

6 

11.697 

1.630 

-0.121 

-1.7152E-02 

I. 31 47E-01 

7 

13.409 

1.630 

-0.418 

-7.6687E-02 

1 . 3748E-01 

8 

15.172 

1.630 

-0.716 

-1.2604E-01 

1 • 1356E-01 

9 

10.634 

2.030 

0.064 

1.0514E-02 

7.5192E-02 

10 

12.147 

2.030 

-0.199 

- 1. 5076E-02 

8. 945 OE- 02 

tl 

13.659 

2.030 

-0.462 

-4.5213E-02 

9.51 80E-02 

12 

15.172 

7.030 

-0.724 

-7. 3495E-02 

8.8895E-02 

13 

11.284 

2.430 

-0.044 

1.4798E-03 

5.820IE-02 

1A 

12.597 

7.430 

-0.277 

-1.2890E-02 

6.5930E-02 

15 

13.909 

2.430 

-0.505 

-2.9458E-02 

7.0358E-0? 

16 

15.222 

2.430 

-0.733 

-4.5829E-02 

6.9411E-02 

17 

11.934 

2.830 

-0. 162 

—2 . 6420E-03 

4.6693E-02 

18 

13.046 

2.830 

-0.355 

-1.1100E-02 

5. 1294E-0? 

19 

14.159 

2.830 

-0.549 

-2.05TQE-02 

5.4357E-02 

20 

15.272 

2.630 

-0.742 

-3. 0227E-02 

5. 5081E-02 

21 

12.504 

3.230 

-0.275 

-4.6022E-03 

3.8628E-02 

2? 

13.496 

3.230 

-0.433 

-9.6934E-03 

4.1492E-02 

23 

14.409 

3.230 

-0.592 

- 1. 5199E-02 

4. 3505E-0? 

2A 

15.322 

3.230 

-0.750 

-2.0869E-02 

4.4550E-O2 

75 

13.233 

3.630 

-0.380 

-5.5263E-03 

3.2774E-02 

76 

13.946 

3.630 

-0.512 

-8.5613E-03 

3.4540E-02 

27 

14.659 

3.630 

-0.635 

-1.1718E-02 

3. 5778E-02 

28 

15.372 

3.630 

-0.759 

- 1 . 501 IE-02 

3.6704E-02 

29 

13.883 

4.030 

-0.501 

-5.9317E-03 

2.8391E-02 

30 

14.396 

4.030 

-0.590 

-7. 60 07E-03 

2.9272E-02 

31 

14.909 

4.030 

-0.679 

-9.3397E-03 

3.0052E-02 

32 

15.422 

4.030 

-0.768 

-1.1 146E-02 

3.0723E-02 

33 

14.533 

4.430 

-0.614 

-3. 9999E— 03 

2.4822E— 02 

3A 

14.046 

4.430 

-0.668 

-6. 8134E-03 

2.5269E-02 

35 

15.159 

4.430 

-0.722 

-7.6499E-03 

2.5687E-02 

36 

15.472 

4.430 

-0.777 

-8.5107E-01 

2. 6080E-02 

37 

15. 183 

4.630 

-0.726 

-5.9198E-01 

2.2003E-02 

38 

15.296 

4.830 

-0 746 

-6. 1620E-03 

2.2141E-02 

39 

15.409 

4.830 

-0.766 

-6.4069E-03 

2. 2276E-02 

AO 

15.527 

4.630 

-0.785 

-6.6540E-03 

2.2409E-02 


(c) Concluded. 
Figure 8.- Continued 



W|Nc 

LIFTING SURFACE C ALCUL AT IHN 

FLOW CONO I T IONS - ALPHA BETA 

26 • OOO 1.0000 

LATTICE ARRANGEMENT - A CMORDWISE ANO 10 SPANUISE ELEMENTS 

REFERENCE QUANTITIES - PLANFORM AREA AVERAGE CHORD 
S - 32.386 CAVG - 4.05 


C I ROOT | CITIPI B/2 DIHEDRAL XCM ZC N YtSEPI ETA1SEPI 

fl .050 0.000 4.000 0.000 0.000 0.000 4.000 1.000 

STATION Y SWEEP LE SWEEP TE 

1 0.400 63.430 0.000 

? 0.800 63.430 0.000 

3 1.200 63.430 0.000 

4 1.600 63.430 0.000 

5 2.000 63.430 0.000 

6 2.400 63.430 0.000 

7 2.800 63.430 0.000 

8 3.200 63.430 0.000 

9 3.600 63.430 0.000 

10 4.000 63.430 0.000 


CNP r.NV CNVS CLP CMP CDI CA CDI/CL2 CXV CVV CRP 

1.00133 0.79972 0.00000 1.07112 -1.11890 0.08809 -0.39037 0.07678 0.19519 -0.34899 1.07474 


X COORDINATES OF CENTERS OF PRESSURE DUE TO 
POTENTIAL LIFT LEADING EDGE SUCTION SIDE EDGE SUCTION OVERALL 
-4.5236 -4.3412 -0.0000 -4.4426 


SPANWISE LOAD DISTRIBUTION 


station 

Y/18/21 

LOCAL CHORD , C CL*C/CL*CAVG 

CL*C/I2*R> 

CX*C/(2*BI 

CY*C/(2*B) 

CSUC*C/C2*8I 

ETA, DEG. 

1 

0.0500 

7.6501 

1.1919 

0.3230 

0.0225 

-0.0443 

0.0497 

90.3978 

? 

0.1500 

6.8503 

1.2551 

0.3402 

0.083? 

-0.1434 

0.1658 

93.5590 

3 

0.2500 

6.0504 

1.2371 

0.3353 

0.0973 

-0.1709 

0. 1966 

93.0956 

4 

0.3500 

5.2506 

1.2043 

0.3264 

0.1106 

-0.1966 

0.2256 

92.7830 

5 

0.4500 

4.4508 

1.1505 

0.3118 

0.1195 

-0.2139 

0.2450 

92.6166 

6 

0.5500 

3.6510 

1.0746 

0.2912 

0.1241 

-0.2230 

0.2552 

92.5148 

7 

0.6500 

2.8511 

0.9743 

0.2641 

0.1246 

-0.2245 

0.2568 

92.4564 

8 

0.7500 

2.0513 

0.8438 

0.2287 

0.1205 

-0.2172 

0.2484 

92.4529 

9 

0.8500 

1.2515 

0.6676 

0.1809 

0.1085 

-0.1952 

0.2233 

92.5091 

10 

0.9500 

0.4258 

0.4008 

0.1086 

0.0768 

-0.1369 

0.1570 

92.7293 

PREDICTED — 

TRAILING 

vorticity separation VORTICITY 





KVLF 

KVSE 

Y 

GAM/V2PI Y 

GAM/V2PI 





4.16155 

- 0.00000 

3.3559 

0.34558 2.1711 

0.22611 






(d) Page 4. 
Figure 8.- Continued. 
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SUMMARY OF VORTEX STRENGTHS AND POSITION AT X - 16.001 


VORTEX GAMNA/2*PI*V V 

1 0.085127 0.397 

7 0.3*5578 A. 386 

3 0.226110 3.201 


2 

1.788 

-0.790 

0.162 



VELOCITY 1 

INDUCED AT 

SPECIFIED 

FIELD POINTS ON 

TAIL (BESKIN) 


X 

Y 

Z 

V/V I INF 1 

W/V (INF ) 

1 

32. 35* 

2.677 

0.000 

0.0000 

1.5098E-01 

2 

35. 32 A 

2.677 

0.000 

0.0000 

1 • 558AE-01 

3 

38.293 

2.677 

0.000 

0.0000 

1.58AAE-01 

A 

A 1 .262 

2.677 

0.000 

0.0000 

l . 58AAE-01 

5 

AA. 231 

2.677 

0.000 

0.0000 

1.5752E-01 

6 

*7. 200 

2.677 

0.000 

0.0000 

1. 52A3E-01 

7 

50.169 

2.677 

0.000 

0.0000 

1.A215E-01 

8 

53. 138 

2.677 

0.000 

0.0000 

1.2595E-01 

9 

3A.376 

*.032 

0.000 

0.0000 

6.8058E-0? 

10 

37.033 

*.032 

0.000 

0.0000 

6.9A81F-0? 

11 

3 9.690 

*.032 

0. 000 

0.0000 

6.9852E-02 

12 

*2.368 

*.032 

0.000 

0.0000 

6.9852F-02 

13 

A5.005 

*.032 

0.000 

0.0000 

6.91 28E-02 

1A 

*7.662 

*.032 

0.000 

0.0000 

6.655AE-0? 

15 

50.320 

*.032 

0.000 

0.0000 

6. 2388E-02 

16 

52.977 

*.032 

0.000 

0* 0300 

5.5986E-02 

17 

36.397 

5.387 

0.000 

0.0000 

3.0792E-O2 

18 

38.7A3 

5.387 

0.000 

0.0000 

3.913AE-02 

19 

A 1.088 

5.387 

0.000 

0.0000 

3.913AE-02 

20 

A3.A3A 

5.307 

0.000 

0.0000 

3 • 90 1 IE-02 

21 

*5. 779 

5.387 

0.000 

0.0000 

3.8363E-02 

22 

*8.125 

5.387 

0.000 

0.0000 

3. 6926E-02 

23 

50. *71 

5.387 

0.000 

0.0000 

3.A796E-02 

2 A 

52.816 

5.387 

0.000 

0.0000 

3.1623E-02 

25 

38. *18 

6.7*2 

0.000 

0.0000 

2. A985F-02 

26 

A0.A52 

6.7*2 

0.000 

0.0000 

2. A985E-02 

27 

A2.A86 

6. 7*2 

0.000 

0.0000 

2 • A985E-02 

28 

**.520 

6.7*2 

0.000 

0.0000 

2. A815E-Q2 

29 

A6.55A 

6. 7*2 

0.000 

0.0000 

2.A262E— 02 

30 

*8. 588 

6.7*2 

0.000 

0.0000 

2.33A7E-02 

31 

50.621 

6.7*2 

0.000 

0.0000 

2.21 17E-0? 

32 

52.655 

6.7*2 

0.000 

0.0000 

2.0356E— 02 

33 

A0.A39 

8.097 

0.000 

0.0000 

1.7323E-0? 

3* 

*2. 162 

8.097 

0.000 

0.0000 

1.7323E-0? 

35 

A3.88A 

8.097 

0.000 

0.0000 

1.72A2E-02 

36 

*5.606 

8.097 

0.000 

0.0000 

I.T01SE-02 

37 

*7.328 

8.097 

0.000 

0.0000 

1 .6621 E-02 

3« 

A9.050 

8.097 

0.000 

0.0000 

1.6029E-02 

39 

50.772 

8.097 

0.000 

0.0000 

1.5265E-02 

AO 

52.A9A 

8.097 

0.000 

0.0000 

1.A228E-02 

A 1 

*2. *53 

9. A* 7 

0.000 

0.0000 

1.2726F-0? 

A2 

A3. 865 

9.AA7 

0.000 

0.0000 

1 • 2668E-02 

A3 

A5.276 

9.AA7 

0.000 

0.0000 

1.2552E-02 

AA 

A6.688 

9.AA7 

0.000 

0.0000 

1.2337F-02 

A 5 

A8.099 

9.AA7 

0.000 

0.0000 

1* 201AE-02 

*6 

*9.511 

9.AA7 

0.000 

0.0000 

1. 16A1E— 02 

A7 

50.922 

9.AA7 

0.000 

0.0000 

1.11 6AE— 02 

A8 

52.33* 

9.AA7 

0.000 

0.0000 

W0537E-02 

A9 

**.*75 

10.802 

0.000 

0.0000 

9.6688E-03 

50 

*5.57* 

10.802 

0.000 

0.0000 

9 . 566 l F-03 

51 

*6.67* 

10.802 

0.000 

0.0000 

9.A378E-03 

52 

*7.77* 

10.002 

0.000 

0.0000 

9.2525E-03 

53 

*8.87* 

10.802 

0.000 

0.0000 

9.0A04E-03 

5A 

*9.973 

10.802 

0.000 

0.0000 

8.7832E-03 

55 

51.073 

10. 00? 

0.000 

0.0000 

8 . 5006 E -03 

56 

52.173 

10.802 

0.000 

0.0000 

8.12A8E-03 

57 

*6.503 

12.162 

0.000 

0.0000 

7.A608E-03 

58 

*7.290 

12.162 

0.000 

0.0000 

7.3732E-03 

59 

A8.077 

12. 162 

0.000 

0.0000 

7.2527E-03 

60 

*8.86* 

12.16? 

0.000 

0.0000 

7.1331F-03 

61 

*9.651 

12.162 

0.000 

0.0000 

6 .99A9E—03 

62 

50. *38 

12.162 

0.000 

0.0000 

6. 83*26-03 

63 

51.22* 

12.162 

0.000 

0.0000 

6.6755E-03 

6A 

52.011 

12.162 

0.000 

0.0000 

6.A612E— 03 

65 

*8.525 

13.517 

0.000 

0.0000 

5.816AE-03 

66 

*9.000 

13.517 

0.000 

0.0000 

5. 7581E— 03 

67 

*9. *75 

13.517 

0.000 

0.0000 

5.6921E— 03 

68 

*9.950 

13.517 

0.000 

0.0000 

5.61 32E-03 

69 

50. *25 

13.517 

0.000 

0.0000 

5.53A8E-03 

TO 

50.900 

13.517 

0.000 

0.0000 

5.A570E-03 

71 

51.375 

13.517 

0.000 

0.0000 

5. 3798E-03 

72 

51.850 

13.517 

0.000 

0.0000 

5.2728E-03 

73 

50.5*6 

1A.872 

0.000 

0.0000 

*• 5559E-03 

7A 

50.709 

1*. 872 

0.000 

0.0000 

*. 5337E-03 

75 

50.073 

l*. 872 

0.000 

0.0000 

A.5117E-03 

76 

51.036 

IA.072 

0.000 

0.0000 

*. A897F-03 

77 

51. 199 

1*. 872 

0.000 

0.0000 

A.A677E-03 

78 

51.363 

IA.872 

0.000 

0.0000 

A.AA58E-03 

79 

51.526 

1*. 872 

0.000 

0.0000 

A. A2A0F-03 

80 

51 .689 

IA.872 

0.000 

0.0000 

*. 3906E-03 


(•) Page 5. 
Figure 8.- Continued. 


131 


VELOCITY INDUCED AT SPECIFIED FIELD POINTS ON TAIL 



X 

Y 

Z 

V/VIINF) 

W/VtINF) 

1 

32,354 

2.677 

0.000 

1.7874E-01 

-2.24026-01 

2 

35.324 

2.677 

0.000 

3- 1029F-01 

-3. 2266E-01 

3 

38.293 

2.677 

0.000 

1.07576-01 

-3.3855E-01 

4 

41.262 

2.677 

0.000 

7.03806—02 

-2.68076-01 

5 

44.231 

2.677 

0.000 

5.8058E-02 

-2.28446-01 

6 

47.200 

2.677 

0.000 

5.3R41E-02 

-2.01386-01 

7 

50.169 

2.677 

0.000 

5.72586-02 

-1.84606-01 

8 

5 3.138 

2.677 

0.000 

6.5031 E-02 

-1.73656-01 

9 

34.376 

4.032 

0.000 

1.4847E-01 

-6.55286-02 

10 

37.033 

4.032 

0.000 

1.49036-01 

-2.93896-02 

11 

39.690 

4.032 

0.000 

2.4147E-01 

3.15396-02 

12 

42.368 

4.032 

0.000 

2.58 1 56-01 

-2.19876-01 

13 

45. 005 

4.032 

0.000 

1.3125E— 01 

-1.82536-01 

14 

47.662 

4.032 

0.000 

l . 0230E— 01 

-1.48196-01 

15 

50.320 

4.032 

0.000 

9.3787E-02 

-1.26 72F-01 

16 

52.977 

4.032 

0.000 

9,501 IE— 02 

-1.13636-01 

17 

36.397 

5.387 

0.000 

1. 10786-01 

3.45 82E-03 

18 

38.743 

5.387 

0.000 

1.0657E-01 

1.17296-02 

19 

41.088 

5.387 

0.000 

1.065 16-01 

2.54126-02 

20 

43.434 

5.387 

0.000 

1. 1753E-01 

4.76286-02 

21 

45.779 

5.387 

0.000 

1.61 106-01 

7.2286F-02 

22 

48.125 

5.387 

0.000 

2. 4523E-01 

3. 5050E— 02 

23 

50.471 

5.387 

0.000 

2.3979E-01 

-6.78616-02 

24 

52.816 

5.387 

0.000 

1.9439E-01 

-8.6268E— 02 

25 

38.418 

6.742 

0.000 

7.29756—02 

2.64556-02 

26 

40.452 

6.742 

0.000 

7.14816-02 

2.85256-02 

27 

42.486 

6.742 

0.000 

7. 03896-02 

3* 1153E-02 

28 

44.520 

6.742 

0.000 

6.9659E— 02 

3,47146-02 

29 

46.554 

6.742 

0.000 

6. 96466—02 

3.96546-02 

30 

48.588 

6.742 

o. ooo 

7.12546-02 

4.6268E— 02 

31 

50.621 

6.742 

0.000 

7.59346-02 

5.43846-02 

32 

57.655 

6.742 

0.000 

8. 18556—02 

5.94786-02 

33 

40.439 

8.097 

0.000 

4.79036-02 

3.0158E-02 

34 

42.162 

8.097 

0.000 

4.74876-02 

3.06216-02 

35 

43.884 

8.097 

0-000 

4.70446-02 

3.10466-02 

36 

45 .606 

8.097 

0.000 

4.64216-02 

3.15246-02 

37 

47.328 

8.097 

0.000 

4.5594E— 02 

3.23006-02 

38 

49.050 

8.097 

0.000 

4.47446-02 

3,36956-02 

39 

50.772 

8.097 

0.000 

4.42456-02 

3.5946E-02 

40 

52.494 

8.097 

0.000 

4.41116-02 

3.80536-02 

41 

42.453 

9.447 

0.000 

3.25066-02 

2. 7844E— 02 

42 

43.865 

9.447 

0.000 

3.23976-02 

2. 78 5 86 - 02 

43 

45.276 

9.447 

0.000 

3.22036-02 

2.77846-02 

44 

46.688 

9.447 

0.000 

3.18676-02 

2.77126-02 

45 

48.099 

9.447 

0.000 

3.13806-02 

2.7769E-02 

46 

49.511 

9.447 

0.000 

3.07896-02 

2.B095E— 02 

47 

50.922 

9.447 

0.000 

3.02316-02 

2.88186-02 

48 

52.334 

9.447 

0.000 

2.97636-02 

2.96216-02 

49 

44.475 

10.802 

0.000 

2.28346-02 

2.41446-02 

50 

45.574 

10*802 

0.000 

2.27826-02 

2.40416-02 

51 

46.674 

10.802 

0.000 

2.26666-02 

2. 3903 E-02 

52 

47.774 

10.802 

0.000 

2.24746-02 

2.37906-02 

53 

48.874 

10.802 

0.000 

2.2204E— 02 

2.37436-02 

54 

49*973 

10*802 

0. 000 

2.18816-02 

2.38326-02 

55 

51.073 

10.802 

0.000 

2.15406-02 

2.40886-02 

56 

52.173 

10.802 

0.000 

2.12166-02 

2.4417E-02 

57 

46.503 

12.162 

0.000 

1.65266-02 

2.05136-02 

58 

47.290 

12.162 

0.000 

1.6472E-02 

2.0416F-02 

59 

48.077 

12.162 

0.000 

1.6389E-02 

2.03356-02 

60 

48.864 

12.162 

0.000 

1.62696-02 

2.02 706 - 02 

61 

49.651 

12.162 

0.000 

1.61226-02 

2.02516-02 

62 

50.438 

12.162 

0.000 

1.5955E-02 

2.02926-02 

63 

51.224 

12.162 

0.000 

1.57756-02 

2.03 806-02 

64 

52.011 

12.162 

0.000 

I.S597E— 02 

2.05 34E-02 

65 

48.525 

13.517 

0.000 

1.22506-02 

1.74056-02 

66 

49.000 

13.517 

0.000 

1.22046-02 

1.73666-02 

67 

49.475 

13.517 

0.000 

1.21456-02 

l. 7343E-02 

68 

49.950 

13.517 

0.000 

1 . 2082E— 02 

1.73346-02 

69 

50.425 

13.517 

0.000 

1.20096-02 

1.73426-02 

70 

5 0.900 

13.517 

0.000 

1.19336-02 

1.73586-02 

71 

51.375 

13.517 

0.000 

1.18516-02 

1.73876-02 

72 

51.850 

13.517 

0.000 

I. 17706-02 

1.74426-02 

73 

50.546 

14.872 

0.000 

9.21796-03 

1.4923E— 02 

74 

50.709 

14.872 

0.000 

9.19906-03 

1.49246-02 

75 

50.873 

14.872 

0.000 

9.17986-03 

1.49266-02 

76 

51.036 

14.872 

0.000 

9. 16006-07 

1.49286-02 

77 

51.199 

14.872 

0.000 

9.13886-03 

1.49336-02 

78 

51.363 

14.872 

0.000 

9.1174E-03 

1.49386-02 

79 

51.526 

14.872 

0.000 

9.09596-03 

1.4943E-02 

80 

51.689 

14.872 

0.000 

9.07506-03 

1.49566-02 


(a) Continued. 
Figure 8.- Continued 



TOTAL INDUCED VELOCITY AT SPECIFIED FIELD POINTS ON TAIL 



X 

Y 

l 

V/V (INF) 

W/VI INFI 

1 

32.354 

2.677 

0.000 

1.78746—01 

-7.30406-02 

2 

35.374 

2.677 

0.000 

3. 10296-01 

-1. 66836-01 

3 

38.293 

2.677 

0.000 

1.07576-01 

-1.00116-01 

4 

41.262 

2.677 

0.000 

7.03806-02 

-1.0962E-01 

5 

44.231 

2.677 

0.000 

5.00586-02 

-7.092 36 -02 

6 

47.200 

2.677 

0.000 

5 . 3841E-02 

-4.89576-02 

7 

50.164 

2.677 

0.000 

5. 7258E—0? 

-4. 2450E-02 

8 

53.138 

2.677 

0.000 

6.50316-02 

-4.77026-02 

9 

34.376 

4.032 

0.000 

1.4847E-01 

2 • 5303E-03 

10 

37.033 

4.03? 

0.000 

1.4903E-01 

4.0092E-02 

11 

39.690 

4.032 

0.000 

2.4147E— 01 

1.0139E-01 

12 

42.340 

4.032 

0.000 

2.58156-01 

-1.5002E-01 

13 

45.005 

4.032 

0.000 

1.3125E-01 

-1.13406-01 

14 

47.662 

4.032 

0.000 

1 • 0230E-01 

-8.16416-02 

15 

50.320 

4.032 

0.000 

9.37876-02 

-6.4330E-02 

16 

52.977 

4.032 

0.000 

9.501 16-02 

-5.76456-02 

17 

36.397 

5.387 

0. 000 

l • 10786-01 

4.22506-02 

18 

38*743 

5.387 

0.000 

1.06576-01 

5.08636-02 

19 

41.088 

5.387 

0. 000 

1.06516-01 

6.45466-02 

20 

43.434 

5.387 

0.000 

1. 1753E-01 

8.66396-02 

21 

45.779 

5.387 

0.000 

1.61106-01 

1 ■ 1065E-01 

22 

48.125 

5.387 

0.000 

2.45236-01 

7.19766-02 

73 

50.471 

5.387 

0.000 

2.3979E-01 

-3.3064E-02 

24 

52.816 

5.387 

0.000 

1.9439E-01 

-5.46446-02 

25 

38.410 

6.742 

0.000 

7.29756-02 

5.14416-02 

26 

40.452 

6.74 2 

0.000 

7 .140 IE-02 

5.35116-02 

27 

42.486 

6.742 

0.000 

7 .03896— 02 

5.61 39E-02 

28 

44.520 

6.742 

0.000 

6. 96596-02 

5.95296-02 

29 

46.554 

6.742 

0.000 

6.9646E-02 

6. 3915E-02 

30 

48.588 

6.742 

0.000 

7. 1254E-02 

6.96 15E-02 

31 

50.621 

6.742 

0.000 

7. 59346-02 

7.65016-02 

32 

52.655 

6.742 

0.000 

8.18556-02 

7.98336-02 

33 

40.439 

8.097 

0.000 

4.79036-02 

4.74816-02 

34 

42.162 

0.097 

0.000 

4.7487E-02 

4. 7944E-02 

35 

43.884 

0.097 

0.000 

4.7044E— 02 

4. 8288F— 02 

36 

45.606 

8.097 

0.000 

4.6421E— 02 

4. 8542E-02 

37 

47.328 

8.097 

0.000 

4.5594E-02 

4.89216-02 

38 

49.050 

8.097 

0.000 

4.4744E-02 

4.9724E— 02 

39 

50.772 

8.097 

0.000 

4.42456-02 

5. 12126-02 

40 

52.494 

8.097 

0.000 

4.41116-02 

5.22826-0 2 

41 

42.453 

9.447 

0.000 

3.25066-02 

4. 05 TOE— 02 

42 

43.865 

9.447 

0.000 

3.2397E-02 

4.0526E-02 

43 

45.276 

9.44 7 

0.000 

3.2203E— 02 

4.03376-02 

44 

46.688 

9.447 

0.000 

3.1867E-02 

4. 0049E— 02 

45 

48.099 

9.447 

0.000 

3 . 1380E-02 

3.97 84E— 02 

46 

49.511 

9.447 

0.000 

3.0789E-02 

3.97366-02 

47 

50.922 

9.447 

0.000 

3.0231E— 02 

3.99826-0? 

48 

52.334 

9.447 

0.000 

2.9763E-02 

4.01 58E-02 

49 

44.475 

10.802 

0.000 

2.28346-02 

3.3813E-02 

50 

45.574 

10.80? 

0.000 

2.27826-02 

3.36076-02 

51 

46.674 

10.802 

0.000 

2.26666-02 

3.33416-02 

52 

47.774 

10.802 

0.000 

2.24746-02 

3.30426-0? 

53 

48.874 

10.802 

0.000 

2. 22046-02 

3.27836-02 

54 

49.973 

10.802 

0.000 

2.18816-02 

3.26166-02 

55 

51.073 

10.802 

0.000 

2. 15406-02 

3.25896-02 

56 

52.173 

10.802 

0.000 

2. 1 2 16E-02 

3.25426-02 

57 

46.503 

12.162 

0. 000 

1. 65266-02 

2.79746-02 

58 

47.290 

12.16? 

0.000 

1.6472E-02 

2.77896-02 

59 

48.077 

12.162 

0.000 

1.63896-02 

2.7588E-02 

60 

48.064 

12.162 

0.000 

1.62696-02 

2.74036-02 

61 

49 . 6*1 

12.162 

0.000 

1.6122E-02 

2. 72*66-02 

62 

50.438 

12.162 

0.000 

1. 5955E-02 

2.71 26E-02 

63 

51.224 

12.162 

0.000 

1.5775E-0? 

2.70556-0? 

64 

52.011 

12.162 

0.000 

1.55976-02 

2.69966-0? 

65 

48.525 

13.517 

0.000 

1.22506-02 

2. 32226-02 

66 

49.000 

13.517 

0.000 

1.2204E-02 

2.3124E-02 

67 

49.475 

13.517 

0.000 

1.21456-02 

2. 3036F— 02 

68 

49.950 

13.517 

0.000 

1. 20826-02 

2.29476-02 

69 

50.425 

13.517 

0.000 

1 , 2009E-02 

2.28766-0 2 

70 

50.900 

13.517 

0.000 

1.19336-02 

2.281 5E-02 

71 

51.375 

13.517 

0.000 

1.18516-02 

2.27676-02 

72 

51.850 

13.517 

0.000 

l • 1770E-02 

2.27146-02 

73 

50.546 

14.872 

0.000 

9.21796-03 

1.9479E-02 

74 

50.709 

14.872 

0.000 

9. 1990E-03 

1.9458E-02 

75 

50.873 

14.872 

0.000 

9.1 798E— 03 

1.9437E-02 

76 

51.036 

14.872 

0.000 

9. 1600E-03 

1 . 94 1 BE— 02 

77 

51.199 

14.872 

0.000 

9. 1388E-03 

1.94006-02 

78 

51.363 

14.872 

0.000 

9. 11746-03 

1.93836-02 

79 

51.526 

14.872 

0.000 

9.0959E-03 

1.9367E-02 

80 

51.689 

14.872 

0.000 

9.07506-03 

1.93466-02 


(e) Concluded. 
Figure 8.- Continued. 
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TAIL 

LIFTING SURFACE CALCULATION 

FLOW CONDITIONS - ALPHA BETA 

16.000 1.0000 

LATTICE ARRANGEMENT - 8 CHORDWISE AND 10 SPANWISE ELEMENTS 


REFERENCE QUANTITIES - PLANFORM AREA AVERAGE CHORD 




S 

- 339.479 

CAVG - 

12.53 

r I ROOT ) 

CITIP) 

B/2 

DIHEDRAL 

xr.M 

ZCN 

25.000 

0.000 

13.550 

0.000 

0.000 

0.000 

STATION 

Y 

SWEEP LF 

SWEEP TE 



l 

1.355 

59.000 

-10.000 



2 

2.710 

59.000 

-10.000 



3 

4.065 

59.000 

-10.000 



4 

5.420 

59.000 

-10.000 



5 

6.775 

59.000 

-10.000 



6 

8.120 

59.000 

-10.000 



7 

9.485 

59.000 

-10.000 



B 

10.840 

59.000 

-10.000 



9 

12.195 

59.000 

-10.000 



10 

13.550 

59.000 

-10.000 



CNP 

CNV 

CNV S 

CLP 

CMP 

COI 

0.60736 

0.24320 

0.00000 

0.70679 -0.75117 

0.06511 


YISEP) ETA ( SEP) 
13.550 1.000 


CA C0I/CL2 CXV 

-0.13223 0.13035 0.06611 


CVV CRP 

-0.10206 0.70978 


X COORDINATES OF CENTERS OF PRESSURE DUE TO 
POTFNTI At. LIFT LFADING EDGE SUCTION SIDE EDGE SUCTION 
-13.4937 -12.7614 -0.0000 


OVERALL 

-13.3043 


SPANWISE LOAD DISTRIBUTION 


STATION 

Y/ 1 9/21 

LOCAL CHORD, C CL*C/CL*CAVG 

CL*C/1 2*61 

CX*C/1 2*B) 

CV*C/ 1 2*81 

CSUC*C/(2*BI 

ET A, DEG. 

1 

0.0500 

23.7530 

1.1359 

0.1855 

0.0076 

-0.0101 

0.0126 

96.2013 

2 

0.1500 

21.2590 

1.1961 

0.1954 

0.0175 

-0.0258 

0.0312 

93.1809 

3 

0.2500 

18.7650 

1.2785 

0.2089 

0.0257 

-0.0397 

0.0473 

91.8755 

4 

0.3500 

16.2709 

1.2946 

0.2115 

0.0331 

-0.0513 

0.0611 

91.8737 

5 

0.4500 

13.7769 

1.2120 

0.1980 

0. 0379 

-0.0589 

0.0700 

91.7272 

6 

0.5496 

11.2921 

1.1069 

0.1808 

0.0406 

-0.0633 

0.0753 

91.6849 

7 

0.6496 

8.7981 

0.9797 

0.1600 

0.0416 

-0.0648 

0.0^70 

91.6808 

8 

0.7500 

6.2949 

0.8243 

0.1347 

0.0404 

-0.0629 

0.0748 

91.6962 

9 

0.8500 

3.8009 

0.6268 

0.1024 

0.0361 

-0.0562 

0.0668 

91.7396 

10 

0.9500 

1.2769 

0.3463 

0.0566 

0.0251 

-0.0387 

0.0461 

91.8927 

PREDICTED— 

TRAILING 

VORTICITY SEPARATION VORTICITY 





KVLF 

KVSE 

Y 

GAM/V2PI V 

GAM/V2PI 





3.20097 

-0.00000 

11.9301 

0.75897 7.6678 

0.16884 






SUMMARY OF VORTEX STRENGTHS ANO POSITION AT X - 55.001 


VORTFX 

GAHMA/2*PI*V 

Y 

Z 

1 

0.085127 

5.660 

1.135 

2 

0.345578 

2.941 

5.499 

3 

0.226110 

3.062 

3.738 

4 

0.758969 

13.930 

0.000 

5 

0.168842 

9.668 

1.918 


(f) Page 6. 

Figure 8.- Continued. 
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ALPHA 

16.00 


— INCIDENCE — 
M WING TAIL 

0.00 10.000 0.000 


WING 


TAIL 


KVLF* KVSF* KVLE* 

0.530 0.000 0.410 


KVSE* W/V LIMIT 

0.000 1.000 


SUMMARY OF FORCE AND PITCHING MOMENT COEFFICIENTS 




CN 

CM 

XCP 


CL 

CD I 

CA 

NOSE... 


POTENTIAL 

VORTEX 

2.850E-03 

1.005E-03 

5.191E-03 

1.855E-03 

4.497F 

4.026F 

00 

00 

2.740F-03 

9.662E-04 

7.857E-04 

2.771E-04 


WING... 

W(ft) 

POTENTIAL 

7.161E-02 

1.022E-01 

1.202E 

01 

6.436E-02 

1 • 845E-02 

1.263E-02 

W ( B ) 

VORTEXf LE 

3.031E-02 

4.355E-02 

1.184E 

01 

2.724E-02 

7.812E-03 

5.345F-03 

WIB) 

VORTEX, SE 

0.000 

0.000 

7. 5 OOF 

00 

0.000 

0.000 

0.000 

BODY... 

BCW 1 

POTENTIAL 

1.682E-02 

2.400E-02 

1.202E 

01 

1.616E-02 

4.635E-03 


8 ( W) 

VORTEX, LE 

9.765E-03 

1.403E-02 

1.184E 

01 

9.386E-03 

2.692E-03 


B(W) 

VORTEX, SE 

0.000 

0.000 

7.500E 

00 

0.000 

0.000 


AFTERBODY... 

1. 755E-03 

4. 069E-04 

3.487E 

01 

1.687E-03 

4. 836E-04 


TAIL... 

TIB) 

POTENTIAL 

5.308E-01 

-8.871E-02 

4.249E 

01 

5. 102E-01 

1 • 463E— 01 

0.000 

TIBI 

VORTFX,LE 

7.590F-02 

-9.776E-03 

4.176E 

01 

7.296E-02 

2.092E-02 

0.000 

TIBI 

VORTEX, SE 

0.000 

0.000 

2.900E 

01 

0.000 

0.000 

0.000 

BODY... 

BIT 1 

POTENTIAL 

7.621E-02 

- 1. 274E-02 

4.249E 

01 

7.326E-02 

2. 101E-02 


BITl 

VORTEX, LE 

1.586E-02 

-2.042E-03 

4.176E 

01 

1.524E-02 

4.370E-0** 


BIT) 

VORTEX, SE 

0.000 

0.000 

2.900E 

01 

0.000 

0.000 


TOTAL CONFIGURATION... 

8. 329E-01 

7.798E-02 

3.751E 

01 

7.942E-01 

2.277E-01 

1 • 797F-02 


CDI/CL*P2 * 3.610 E- 01 


Cg) Page 7. 
Figure 8.- Continued. 
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SUMMARY OF FREE VORTFX TR A JECTGRI F S 


VPPTEX l POOY VORTFX FROM NOSE 

VORTFX 2 W TNG TRAILING VORTFX 

VORTFX 3 WING SEPARATION VORTFX 

VPPTFX GAMMA/ ?*P I *V 

1 0.085127 

2 0.345578 

3 0.226110 


X DX A S RO OA/DX 

7.5000 0.5000000 0.8400 0.8400 0.3400 0.0791 

VORTFX SIGMA(REAL) (IMAG) 
l 4*06 5 F— 0 1 1.206F 00 


X OX AS RO OA/OX 

8.5000 0.2500000 0.9175 1.5301 1.0401 0.0762 

VORTFX SIGMA(RFAL) (IMAG) 

1 3.P93F-01 1.338F 00 


X DX A S RO PA/OX 

9.5000 0.5000000 0.9925 2.0302 1.2577 0.0750 

VPPTEX SIGMA(RFAL) (IMAG) 

1 3. 735F-01 1 • 44 IF 00 

• • • 

X nx A S RO OA/DX 

16.0000 0.5000000 1.4600 1.4600 1.4600 0.0625 

VORTFX S I GMA (PEAL) (IMAG) 

1 3.973F-01 1.7B8F 00 


X OX A S RO DA/DX 

16.0010 0.5000030 1.4601 1.4601 1.4601 0.0625 

VORTEX SIGMA(REAL) (IMAG) 

1 3.973E-01 1.788E 00 

2 4.386E 00 -7.901E-01 

3 3. 201 E OC 1.622F-01 


X DX A S RO DA/DX 

17.0010 0.2500000 1.5231 1.5201 1.5201 0.0587 

VPPTFX SIGMA (PEAL) (IMAG) 

1 4.430E-01 1.810E 00 

2 4.4RIF 00 — 4.6 10E-0 l 

3 3.117F 00 l • 570F-01 

• • • 


(h) Page 8. 

Figure 8.- Continued. 
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29.0010 Q.S0C0300 


1 .9350 


0.0150 


1.9350 


1.9356 


VnPTTX SIGMA(PEAL) 


1.909F 00 
3.55?F CO 
4. 8 8 OF 00 


( IMAG) 
1.010E 00 
2.539F 00 
1.847E 00 


30.0010 


OX 

0.5000000 


A 

1.9500 


2.5653 


RO 

2.0233 


VHPTEX SIGMA(RFAL) (IMAG) 

1 2.051E 00 8.237E-01 

2 3 . 489 F 00 2.6C0E 00 

3 4.950F OC 2.284E 00 


31.0010 


nx 

0.5000000 


A 

1.9650 


s 

3. I960 


RO 

2 .201 R 


VORTFX SIGMA(REAL) (IMAG) 

1 2.165E 00 6.690E-01 

2 3.457E 00 2.648E 00 

3 4.957F 00 2.733E 00 


61.0010 


DX 

0.5000000 


A 

1.8991 


15.2196 


RO 

7.7283 


VORTEX SIGMA(PFAL) (IMAG) 

1 5.563F 00 4.832E-01 

2 3.557F 00 4.993E 00 

3 2.528F 00 3.710F 00 


62.0010 


DX 

C. 5000000 


A 

l. 8626 


13. 0185 


RO 

6.6425 


VOPTFX SIGMA(PFAL) (IMAG) 

1 5.660E 00 5.058E-01 

2 3.410F 00 5.087E 00 

3 2.668F 00 3.586E 00 


53.0010 


DX 

0.5000000 


A 

1.8163 


5 

7.3942 


RO 

3.9202 


VnPTFX SIGMA(PEAL) I I MAG ) 

1 5.655F 00 5.482E-01 

2 *.?43F 00 5.197E 00 

3 2. POOF 00 3.536E 00 


64.0010 


DX 

0.1260000 


A 

1.7699 


S 

1.7699 


RO 

1. 7699 


VORTEX SIGMAIPEAL) ( IMAG) 

1 5.607F 00 8.390E-01 

2 3.085E 00 5.350E 00 

3 2.923C 00 3.603E 00 


55.0010 


DX 

0.2500000 


A 

1.6999 


S 

1.6999 


RO 

1.6999 


VORTEX 

SlGMA(REAL) 

(IMAG) 

1 

5.660F 

00 

I.135E 

00 

2 

2 • 941 E 

00 

5.499E 

00 

3 

3.062E 

00 

3.738E 

00 


DA/DX 

0.3143 


DA/DX 

0.0135 


OA/DX 

-0.0343 


DA/DX 

-0.0412 


DA/DX 

- 0.0497 


DA/DX 

-0.0582 


DA/OX 

-0.0653 


(h) Concluded. 
Figure 8.- Concluded 



PART IV - SUPSON COMPUTER PROGRAM 


Introduction 

This computer program predicts the static longitudinal aerodynamic 
characteristics of wing-body-tail combinations at supersonic speeds. It is 
an extension of the method of reference 2 to angles of attack for which 
symmetrical body vortices are shed from the nose of the configuration and 
leading-edge and side-edge separation vortices are shed from the wing and 
tail. A lifting-surface method (refs. 12 and 13), modified to include 
interference velocity fields in the form of induced camber and modified to 
compute leading-edge suction distributions, represents the wing and tail 
surfaces. 

The program is written in FORTRAN IV for the IBM 360 series machines. 
No tapes, drums, or disks other than the standard input/ output units are 
required. Minor changes are required to run the program on other machines 
such as the CDC 6600. Typical running time on the IBM 360/67 for a wing- 
body-tail configuration is approximately 1 minute per angle of attack. 
Actual running time depends on the number of panels representing the 
lifting- surf aces. Some specific running times are noted in the discussion 
of the sample cases. 

The following sections present descriptions of the method, program, 
input, and output. A program listing and sample cases are included. The 
algebraic notation used in this section is the same as that used in 
reference 1. A list of symbols from reference 1 is included. 

List of Symbols 

jr aspect ratio 

a local body radius 

b semispan 

C D induced drag coefficient 

C lift coefficient, ~ 

L qS 

C„ axial-force coefficient 

A 

C pitching-moment coefficient, -—r 
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c 


c 

n 


c 

s 


c 

X 


c 

y 


normal-force coefficient, 
local chord 

TE 

section normal-force coefficient, — / dx 

c J q 

LE 

section leading-edge suction coefficient 
x-direction section suction coefficient 
y-direction section suction coefficient 


K* vortex-lift ratio 

v 

L lift force 

£ reference length 

M pitching moment about center of moments, or free-stream Mach 

number 

N normal force 

Ap static pressure difference between lower and upper surfaces of 

lifting surface 

q free-stream dynamic pressure 

r body radius 

r radius of base of nose 

N 

S reference area 

s semispan of lifting surface 

u,v,w perturbation velocities along x,y,z directions, respectively 

V free-stream velocity 

x,y,z configuration coordinates with origin at body nose 


x ac 

X 

location 

of 

aerodynamic center 

X HL 

X 

location 

of 

lifting-surface hinge line 

X m 

X 

location 

of 

center of moments 
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X 

z 

a 

3 


HL 


B 


n 


6 

a 

A 

P 


N 


A 

avg 
B (T) 
B(W) 
e 

HL 

LE 

N 

P 

root 

SE 
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x position for onset of separation from body nose 

center-of-pressure location 
height of lifting-surface hinge line 

body angle of attack 

Vm 2 - 1 

right body-vortex strength, positive counterclockwise when 
viewed from rear of configurations 

n'th separation vortex strength on right wing panel 
trailing- vortex strength on right wing panel 

lifting-surface deflection angle, positive trailing edge down 
nose angle, degrees 

sweep angle 
density 

complex vortex position, y + iz 

Subscripts 

afterbody 

average 

body in presence of tail 
body in presence of wing 
tail or empennage 
hinge line 
leading edge 
nose 

potential 
root chord 
side edge 



T(B) 

tail in presence 

of body 

TE 

trailing edge 


tip 

tip chord 


V 

vortex 


W(B) 

wing in presence 

of body 

w 

wing 



Description of Method 

A brief description of the method is presented herein* The user 
should consult reference 1 for a complete description and details of the 
theoretical approach. 

An axisymmetric nose at some moderate angle of attack sheds a symmetric 
pair of body vortices. These shed body vortices, whose strength and initial 
position are determined from data correlations, are tracked downstream past 
the wing using slender-body techniques in the crossflow plane. One excep- 
tion to this trajectory calculation is that the vortices move parallel to 
the wing if the leading edge is supersonic. The vortex -rinduced velocities 
are computed at the wing control points and combined with the Beskin upwash 
induced by the body to obtain the total upwash induced on each wing panel. 
This, added to the free-stream contribution, results in a total local 
incidence angle distribution over the wing. 

The wing is modeled by a constant pressure panel scheme in the form 
of a prediction program obtained from R. Carmichael of the NASA, Ames 
Research Center. The result from this method is a pressure distribution 
on the wing. An equivalent circulation distribution on the wing is ob- 
tained from the pressure loading and this is used to compute a distribution 
of leading-edge suction and side-edge suction (if present) and their associ- 
ated vortex positions and strengths. The vortex lift on the leading edge 
and side edge is obtained from the suction distribution through the 
Polhamus vortex-lift analogy (ref. 9) and correlation curves. 

The trajectories of the body vortices, the wing trailing vortex, and 
the wing leading-edge separation vortices are computed downstream past 
the afterbody and horizontal tail. These trajectories are computed in the 
crossflow plane considering mutual interference between the vortices and 
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interference from their images in the body. If the tail leading edge is 
supersonic, the vortices move parallel to the tail surface. The induced 
velocity field on the tail panels is computed, and the tail loading is 
obtained in a manner similar to that just described for the wing. The 
forces on the body in the presence of a wing and tail are computed by the 
method of reference 2. The free vortex- induced forces on the nose and 
afterbody are computed using the method of Sacks (ref. 10) . 

The forces and moments on the entire configuration are obtained by 
summing the contributions of the various components. The forces are 
resolved into normal force and axial force (excluding frictional drag and 
wave drag) , and lift and induced drag. 

Description of Program 

The SUPSON computer program consists of a main program, nine function 
subprograms, and twenty-five subroutine subprograms. The main program 
(SP01) accepts most of the input, prints a portion of the output, and 
generally directs the flow of the calculation. The subprograms provide 
specific services to the main program during the calculation procedure. 

Since the lifting-surface portion of the program was a separate wing-alone 
prediction program which was incorporated into the prediction scheme, there 
are several options present in this part of the program which are not used. 
These were left in the program so that the original lifting-surface proced- 
ure could be left intact as much as possible; however, this does result in 
portions of subprograms being carried along but never used. These surplus 
calculations are short and do not cause a large core storage penalty by 
their presence. The following is a list of subprograms and their general 
purpose . 

Subroutine SSWING is the former main program of the Carmichael- 
Woodward lifting-surface method. It accepts description of the wing (or 
tail) geometry through namelist input and then sets up locations of 
control points and divides the wing up into the specified number of constant 
pressure panels. It then calculates an array of influence functions which 
are used later in computing wing-loading distributions. 

Subroutine TABLE sets up a table of coordinates describing the config- 
uration at specific axial stations required in the vortex trajectory 
calculations. The body radius at the required stations are obtained by 
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linear interpolation in the input table of body coordinates. Body slopes 
are computed by a simple differencing scheme. 

Subroutine FORCE computes the loading on the lifting surface consi- 
dering both the geometric camber and induced camber. The bound circulation 
distribution corresponding to the potential loading is computed as is the 
strength and position of the equivalent trailing vortex. 

Subroutine EDGFRC computes the distribution of leading-edge and side- 
edge suction. The leading-edge and side-edge vortex lift is computed along 
with the strength and lateral position of the associated separation vortex. 

Subroutine MATRIX computes the aerodynamic influence matrix for each 
lifting surface. 

Subroutine CNVTX computes the vortex- induced force and center of 
pressure on the afterbody using the method of reference 10. 

Subroutine BDYVTX uses tables derived from data correlations to look 
up the strength and position of the pair of symmetric vortices shed from 
the body nose. 

Subroutine CNVNZ computes the nose vortex- induced normal force and 
center of pressure on the body nose using the method of reference 10. 

Subroutine SHAPE does a table look-up for the body radius and slope 
and local lifting-surface semispan at any prescribed axial station. 

Subroutine FILL is used to fill in intermediate locations in an array 
using a linear interpolation procedure. 

Subroutine ZSECT is a general airfoil section calculation routine 
which is used only when the lifting surfaces are assumed to have finite 
thickness. The available section shapes are specified by the variable 
SECT as described in the discussion of input preparation. 

Subroutine TCOMP computes the pressure function for a wedge-shaped 
airfoil of triangular planform. This routine is used only when the 
lifting surface has a finite thickness. 

Subroutine TAINT is used to interpolate in tables of airfoil section 
thicknesses and surface slopes. 

Subroutines LINEQS and SOLVE are used together to solve the set of 
linear simultaneous equations for the loading on the lifting surfaces. 
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Subroutine COMP computes the downwash function fora uniformly loaded 
triangular shaped panel. 

Subroutine EXTVEL computes the free vortex- induced velocity at wing 
or tail control points. 

Subroutine TRJTBY computes the trajectories of the free vortices past 
the configuration using the subroutines FCT, OUTP, HPCG, and SHAPE. 

Subroutine FCT computes the derivatives for the equations of motion 
for each free vortex. 

Subroutine HPCG is a predictor-corrector integration package which 
uses a Runge-Kutta starting procedure. 

Subroutine OUTPUT stores the vortex positions in a table at specified 
intervals in x. When necessary, some diagnostic information on the vortex 
trajectories is available as optional output. 

Subroutine ZVTX determines the vertical position of the leading-edge 
separation vortex using a table look— up of correlated data for delta wings. 

Subroutine WPANL calculates the upwash induced by the constant 
pressure panels comprising the lifting surface. 

Subprograms KFACT, CH1416, EQ14, CHRT8, EQ21, EQ24, EQ24L, EQ26, EQ30, 
EQ30L, and EQ31 are used to compute the lift and center of pressure on the 
body in the presence of a lifting surface. These subprograms perform the 
same functions in SUPSON as they do in R1307 and the individual descriptions 
are not repeated here. 


Description of Input 

Variable definitions.- The format of the input cards for the SUBSON 
program is shown in figure 1. In this figure the program variable name 
is shown as well as the card columns in which the value is punched and the 
format in which it is punched. The following is a table of the input 
variables along with the algebraic synfool where applicable. The variable 
is defined and its limits shown where necessary. The algebraic notation 
used in defining the configuration is shown in figure 2. A discussion 
of the preparation of the input is presented in the section following 
the table. 
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PROGRAM 

NOTATION 

Item 1 

NHEAD 

NTBL 

NPRINT 


MPRINT 

NOSEV 

NREGNW 

NALFW 

NSEPW 

NREGNT 


ALGEBRAIC 

NOTATION DEFINITION 


Number of heading cards. 

Number of entries in table of body- 
coordinates . 5 < NTBL < (96 less the 

total number of columns on both wing 
and tail) . 

Index controlling optional output: 

NPRINT = -2 Minimum output. Final aero- 
dynamic characteristics only. 

= -1 Abbreviated output. Final 

aerodynamic characteristics 
plus details of lifting- 
surface results. 

= 0 No optional output. 

= 1 Some optional output. 

= 4 Large amount of optional 
output . 

> 4 All optional output including 

diagnostic information (not 
recommended for general use) . 

Index controlling diagnostic output: 

MPRINT - 0 No additional output. 

= 1 Output vortex trajectories 

during calculations. (This 
option should be used only 
if program has terminated 
execution during a previous 
trajectory calculation.) 

NOSEV = 0 No nose separation vortex 
pair. 

= 1 Nose separation vortex pair 
included . 

Number of regions into which wing is 
divided. 1 < NREGNW < 20 

NALPW = 0 Uncambered and untwisted wing. 

> 0 Cambered and/or twisted wing, 

local angles must be input. 

Number of leading-edge separation vortices 
shed from wing. 1 <C NSEPW <_ 2 

Number of regions into which tail is 
divided. 1 £ NREGNT < 20 
NREGNT — 0 No tail. 
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PROGRAM 

NOTATION 

NALPT 

NSEPT 

NBODY 

NAFT 

Item 2 
TITLE 

Item 3 

EM 

REFS 

REFL 

XM 

THETAN 

DXOUT 

DXI 

Item 4 
XBDY 

RBDY 


ALGEBRAIC 

NOTATION 


DEFINITION 


NALPT = 0 Uncambered and untwisted tail. 

> 0 Cambered and/or twisted tail, 
local angles must be input. 

Number of leading-edge separation vortices 
shed from tail. 1 <_ NSEPT <_ 2 

Index controlling body upwash: 

NBODY = 0 Body upwash included in wing 
and tail interference 
calculation. 

< 0 Body upwash not included. 

Index controlling presence of afterbody 
behind rear lifting surface. 

NAFT = 0 No afterbody. 

= 1 Afterbody included. 


Any alphabetical or numerical identifi- 
cation information. 


M 

S 

l 


x 


r 


Mach number, M > 1.0. 

Reference area. 

Reference length. 

x-location of center of moments 

Nose semiapex angle, degrees. 

x— increment in table of free vortex 
trajectories. DXOUT < RAVGW , 
typically. 

Maximum integration interval for vortex 
trajectory calculation. DXI ^ DXOUT, 
typically. 


x-stations at which body coordinates are 
defined. 

Body radius at above x-station. 
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PROGRAM 

ALGEBRAIC 


NOTATION 

NOTATION 

DEFINITION 

Item 5 



RAVGW 

r 

avg 

^w 

Average body radius at wing. 


YSEPW 

y 

se Pw 

Spanwise location of second leading-edge 
separation vortex. If NSEPW = 1, 

YSEPW = b/2 . 

XHLW 

Xrl 

W 

x-coordinate of wing hinge line at wing- 
body juncture. 

ZHLW 

Z HL 

w 

z-coordinate of wing hinge line at wing- 
body juncture. 

XWCP 

X CP 

w 

x-coordinate of alternate wing center-of- 
pressure location. 

Item 6 


Item 6 is made up of NREGNW namelist decks. 

$ INPUT 


Namelist identification. 

PER 


Location of wing control point in fraction 
of panel chord. Typically, PER = 0.95. 

ROWS 


Number of chordwise divisions into which 
wing region is divided. 1 < ROWS <C 10 

COLS 


Number of spanwise divisions into which 
wing region is divided. 1 COLS 20 


NOTE : Sum of- ROWS X COLS over all wing regions must not exceed 100 and 

total number of columns on wing semispan must not exceed 20. 


ROOTLE 


x-coordinate of leading-edge of root chord 
of wing region now under consideration. 

ROOTTE 

x mn 

TE w 

x-coordinate of trailing edge of root 
chord. x TE > x^ 

ROOTY 


y-coordinate of root chord. Typically, 
ROOTY = RAVGW in first wing region. 

TIPLE 

"“tip 

x-coordinate of leading edge of tip chord. 

TIPTE 

^tip 

x-coordinate of trailing edge of tip 
chord. x TE > x^ E 

TIPY 


y-coordinate of tip chord of this wing 
region. If NREGNW = 1, TIPY = b/2. 
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PROGRAM 

NOTATION 

ALGEBRAIC 

jioraTioN. 

DEFINITION 

TCROOT 

t//c I root 

Root chord thickness ratio. May be equal 
to zero. 

TCTIP 

t/C |tip 

Tip chord thickness ratio. May be equal 
to zero. 

SECT 


Specification of airfoil section. 
SECT = 1 Parabolic arc. 

= 2 Double wedge. 

= 3 30-70 hexagon. 

= 4 Wedge. 

= 5 NACA 000X 
= 6 NACA 6400X 
- 7 NACA 6500X 
= 8 RAE 101 

$END 


End of namelist. 

Item 7 



ALPHAW 

a 

w 

Local angle of wing mean surface due to 
camber and twist, radians. 

Item 8 



RAVGT 


Average body radius at tail . 

YSEPT 

y sep e 

Spanwise location of second leading-edge 
separation vortex. If NSEPT = 1, 

YSEPT = b/2 . 

XHLT 


x-coordinate of tail hinge line at wing- 
tail juncture. 

ZHLT 

*». 

z— coordinate of tail hinge line at wing- 
tail juncture. 

XTCP 


x— coordinate of alternate tail center— of— 
pressure location. 

Item 9 


Item 9 is made up of NREGNT namelist decks. 

$ INPUT 


Namelist identification. 

PER 


Location of tail control point in fraction 
of panel chord. Typically, PER = 0.95. 
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PROGRAM 

NOTATION 

ROWS ^ 

COLS 

ROOTLE 

ROOTTE 

ROOTY 

TIPLE \ 

TIPTE ( 

TIPY 

TCROOT 

TCTIP 

SECT 

$END ^ 

Item 10 


ALPHAT 


Item 11 
NDEX 


ALPHAD 

ALPIW 

AKVLW1 

AKVLW2 

AKVSW 

ALPIT 

AKVLTl 


ALGEBRAIC 

NOTATION DEFINITION 


Corresponding tail variables having same 
limits and restrictions as those on wing. 


a 


e 


Local angle of tail mean surface due to 
camber and twist, radians. 


Index controlling next case of input. 

NDEX = 1 Execute program using variables 
on this card. 

= 0 Ignore this card and return to 
beginning for new case. 

a Angle of attack of configuration, degrees. 


w 


Incidence angle of wing relative to body 
axis, degrees. 


K 



Fraction of leading-edge suction converted 
to lift in inboard wing region. 
(0 <K* le <1.0) 



Fraction of leading-edge suction converted 
to lift in outboard wing region. 

(0 < Kv le < 1.0) 



Fraction of side-edge suction converted 
to lift on wing. (0 < K* SE < 1.0) 


6 

e 


Incidence angle of tail relative to body 
axis, degrees. 



Fraction of leading-edge suction converted 
to lift in inboard tail region. 

(0<K* le <1.0) 
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PROGRAM 

NOTATION 

ALGEBRAIC 

notation 

DEFINITION 

AKVLT2 


Fraction of leading-edge suction converted 
to lift in outboard tail region. 

(o < k* le < 1.0) 

AKVST 

*V. 

Fraction of side-edge suction converted to 
lift on tail. (0 < K* SE < 1.0) 

WLIMIT 

(w/V) „ 
v ' max 

Limit on vortex- induced velocities at wing 
and tail control points. 

(0 < WLIMIT < 1.0) 


input preparation.- A discussion of the input variables is presented 
in this section as an aid in the preparation of the input data deck. If 
a configuration has one set of lifting surfaces, this is denoted the wing 
regardless of its axial location on the body and data describing this 
lifting surface are input as wing data. If there are two lifting surfaces, 
the set nearest the nose is the wing and the aft set is the tail. For 
example, in a canard-body-wing configuration, the canard data are input 
as wing data and the wing data are input as tail data. In the following 
discussion, necessary geometric relations are illustrated in figure 2. 

Item number 1 of the input data deck (fig. 1) is a card containing 
indices specifying particular program options. NHEAD indicates the number 
of identification cards following in Item 2. NTBL is the number of entries 
in the table describing the body shape. N PRINT is an index which deter- 
mines the quantity of output obtained from the program. Typically, this 
number is zero, but varying degrees of additional output can be obtained 
by increasing this index from 1 to 4. An abbreviated output summary is 
obtained by setting N PRINT = -1 or -2. MPRINT is a special index con- 
trolling the quantity of output obtained during vortex trajectory calcu- 
lations. Typically, this index is zero, but if it is greater than zero, 
the trajectory coordinates are printed as they are computed. This addi- 
tional trajectory information is useful only if the program fails to 
compute a trajectory for some particular case, and the additional output 
may give some clue as to why the calculation failed. This option should 
only be used after a computational problem has been discovered. 

The index NOSEV specifies whether there is (NOSEV = 1) or is not 
(NOSEV = 0) a symmetrical pair of vortices shed from the nose of the 
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configuration. If the nose angle (THETAN) is less than or equal to four 
degrees, NOSEV is automatically set equal to zero. 

The following three indices are associated with the wing. NREGNW is 
the number of spanwise regions describing the wing. This number is 
typically equal to zero unless there are breaks in sweep or it is desired 
that the panel spacing be different in different wing regions. NALPW 
specifies if wing camberline slopes are nonzero and must be input (NALPW = 1) 
or are zero and need not be input (NALPW = 0) . NSEPW specifies the number 
of leading-edge separation vortices shed from the wing. This index must 
be either 1 or 2, and if NSEPW =2, a special value of YSEPW is required 
in Item 5. 

The provision for multiple separation vortices is included to handle 
wings with breaks in leading— edge sweep. It has been observed that leading- 
edge separation vortices are shed from the wing regions inboard and out- 
board of the break. A maximum of two leading-edge separation vortices are 
allowed. 

The next three indices in Item 1, NREGNT , NALPT, and NSEPT are tail 
indices analogous to the preceeding wing indices and subject to exactly 
the same restrictions. If no tail is present, all three must be zero. 

The next index in Item 1, NBODY, determines whether the upwash field 
around the body is included (NBODY = 0) in the wing-and tail-interference 
calculations or not included (NBODY < 0) . This index is used to determine 
the magnitude of the body-interference effect and generally should be 
equal to zero. 

The last index, NAFT, is used to specify the presence of an afterbody 
behind the aft lifting surface. If there is an afterbody, NAFT = 1, and 
if not, NAFT = 0. 

Item 2 is a group of NHEAD cards containing identification information 
which is printed on the first page of output. The information on the last 

of this group is retained and printed as a heading card at various 
points in the output. 

Item 3 is a single card containing EM, the free— stream Mach number 
(greater than 1) ; REFS and REFL, the reference area and reference length, 
respectively; XM, the x-coordinate of the center of moments; and THETAN, 
the nose semiapex angle in degrees. The final two variables are associated 
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with the free vortex trajectory calculations. DXOUT is the approximate 
increment in x at which trajectory coordinates are stored for use in 
induced velocity calculations. A lower limit for this variable is about 
0.5 percent of the overall length of the body because of storage limitations 
Typically, a reasonable value for DXOUT is about one half the maximum 
radius of the body. DXI is the initial integration interval for the tra- 
jectory calculations. The integration package will cut the interval m 
half if necessary for reasonable accuracy, and this halving process can 
occur ten times before the program automatically terminates execution with 
an appropriate message. If the input value of DXI is too large, the 
program will stop because of unacceptable accuracy, and if DXI is too 
small, the running time will become large. Experience has indicated a 
value of DXI between 2 and 5 percent of the body length will work for 
most cases. Under rare circumstances when two vortices get very close 
together or when a vortex gets very near the wing or body, a smaller value 
of DXI may be required. DXI should not be larger than DXOUT. 

Item 4 is a group of NTBL cards describing the body shape. Each card 
contains an x-station, XBDY , and the corresponding body radius, RBDY. The 
cards should be in ascending order in x and there should be less than 
75 cards in this item in a typical run. The program internally sets up 
its own table of coordinates which is stored in the XBDY and RBDY arrays 
and is limited to 100 entries. A good rule of thumb to follow in inputting 
Item 4 is the following. 

NTBL < (96 — Total number of spanwise columns on wing and tail) 

Some care is required when describing the body shape via XBDY and RBDY. 
Linear interpolation is used throughout; therefore, if the body shape is 
changing rapidly, more points are required. There should be a minimum 
of five entries in the nose region ahead of the wing and there must be 
entries at x-stations identically equal to XWLE and XWTE, and XTLE and 
XTTE if a tail is present. The last entry in the table must be greater 
than XWTE or XTTE, whichever is greater, by an amount not less than DXI. 

If the body is made up of a nose section followed by a cylindrical after- 
body, there should be two points on the cylinder very close together near 
the beginning of the cylinder. Points on a cylinder can be spaced large 
distances apart, but if the cylinder is followed by a section with changing 
radius, the last two points on the cylinder should be close together. 
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Item 5 contains geometric information for the wing. RAVGW is the 
average body radius in the vicinity of the wing. YSEPW is the y-station 
at which the wing is assumed broken for purposes of having two leading- 
edge separation vortices. If NSEFW = 1, then YSEPW must be equal to the 
wing semispan. If NSEPW = 2, YSEPW must have some value greater than 
RAVGW and less than the wing semispan. The chosen value should coincide 
with one of the breaks in the panel layout. It is advised that there be 
at least three columns of panels on either side of YSEPW to achieve some 
reasonable accuracy in the separation vortex strength and position calcu- 
lation. The next two variables, XHLW and ZHLW, are the coordinates of 
the wing hinge line at the wing-body juncture. The last variable, XWCP, 
is an experimental center-of-pressure location for the wing potential lift 
which may be used if desired. If this is not to be used, XWCP must be 
identically zero. 

Item 6 is a namelist (INPUT) describing the geometry of the wing and 
the panel arrangement. There must be one complete namelist deck for each 
region of the wing (NREGNW) and they must be in order from wing root to 
wing tip. The first variable, PER, is the chordwise location of the 
control point of each constant pressure panel presented as a fraction of 
the local panel chord. This number is typically 0.95 for supersonic flow. 
The next two variables, ROWS and COLS, are the number of chordwise and 
spanwise divisions into which the wing region is divided. Thus, ROWS x COLS 
is the number of panels in the wing region under consideration. There are 
certain limitations on the size of these numbers, as the total number of 
panels on the entire wing semispan cannot exceed 100 and the total number 
of spanwise divisions on the entire wing semispan cannot exceed 20. 

The following six variables specify the planform geometry of the wing 
region. ROOTLE is the x-station of the leading edge of the root chord of 
the region, ROOTTE is the x-station of the trailing edge of the root chord, 
and ROOTY is the y-station of the root chord. If this is the first wing 
region, these values correspond to the actual wing root chord where ROOTLE 
and ROOTTE are the intersections of the leading and trailing edges with 
the body, and ROOTY has the same value as RAVGW in Item 5. If NREGNW > 1, 
these values specify the chord of the inboard side of the wing region, in 
a similar manner, TIPLE and TIPTE are the x-stations of the tip leading 
and trailing edges, respectively. TIPY is the y-station of the outboard 
side of the region. If the region under consideration is the last (or only) 
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region of the wing, TIPY must be equal to the semispan of the wing. Note 
that on wings with multiple regions and continuous leading and trailing 
edges, the tip chord of one region is identical to the root chord of the 

adjacent region. 

The remaining variables in the namelist are optional. TCROOT and 
TCTIP are the thickness-to-chord ratios of the root and tip of the wing 
region. These values may be input as zero or omitted entirely from the 
deck. Since wing thickness is not a prime consideration in the calculation 
procedure, it is suggested that these values be omitted. If wing thickness 
is to be input, a value of SECT must be input to specify the wing section 
shape. The options available for airfoil sections are noted in the tables 
of variable definitions. This ends the namelist description of the wing 
region. All other regions should follow in the input deck immediately; 
therefore, there will be NREGNW sets of the namelist INPUT decks. 

Item 7 is optional input and is included only if the wing is cambered 
or twisted (NALPW > 0) . If such is the case, ALFHAW, the local panel angle, 
a of each element of the wing must be input. These angles, in radians, 
are relative to the wing root chord and thus do not include any wing 
incidence. There are eight values per card and the angles are input from 
leading edge to trailing edge, from wing root to tip. There are a total 
of NPANLW values, where NPANLW is the total number of panels on the wing 
and is computed internally in the program. 

Values of ALPHAW are obtained as follows. Consider the sketch in 
figure 2 (b) which shows the cambered and twisted section of the lifting 
surface at some spanwise station. At point P, corresponding to a control 
point on the wing mean surface, a tangent to the wing mean surface is 
constructed, which makes an angle a_g with the wing root chord. The 
positive sense of is shown in this figure. The input value required 

is a in radians. Near the leading edge of the section shown in fig- 
ure 2(b), a,£ is negative. Item 10 completes the input description of the 

wing. 

If NREGNT = 0; that is, no tail is present, go directly to Item 11. 

If a tail is present. Item 8 is the next card in the input deck. The 
variables in Item 8, RAVGT, YSEPT, XHLT, ZHLT, and XTCP are analogous to 
the wing parameters in Item 5 and subject to the same restrictions and 
limitations . 
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Item 9 is a namelist INPUT describing the tail and there must be 
NREGNT decks of this namelist. The tail variables in this namelist are 
in direct correspondence with the wing variables described in Item 6. 

If NALPT is not zero, the local panel incidence angles in radians of 
the tail are input in Item 10 in an analogous manner to item 7 for the 
wing. 

Item 11 is a group of cards, one card for each run, which specifies 
the variables which are considered changeable for a given geometric con- 
figuration. The first entry on the card is the index NDEX which is simply 
used to control the stacking of additional cases. NDEX =1 on each card 
represents a new flow condition. If NDEX = 0, the card is ignored and 
the program returns to read in a completely new case beginning with item 1. 
Thus, a blank card is used to separate calculations with different configu- 
rations. When NDEX ^ 0, the next value on the card is the configuration 
angle of attack in degrees, ALPHAD, taken as the angle between the axis of 
the body and the free-stream velocity. The second quantity is the incidence 
angle of the wing root chord in degrees, ALPIW. Its sense is such that a 
positive incidence is a leading-edge up condition. The next three vari- 
ables are the K* factors which relate the actual realized vortex lift 
from the leading and side edges to that which is theoretically available. 
AKVLWl is the fraction of leading-edge separation vortex lift which is 
obtained on the inboard portion of the wing if NSEPW =2 or on the entire 
wing if NSEPW =1. If NSEPW = 2, AKVLW2 is the leading-edge suction 
factor for the outboard portion of the wing. The values of AKVLWl and 
AKVLW2 are obtained from the correlation curve in figure 3. This figure 
is reproduced from figure 9 of reference 1. The data on which the curve 
is based are for sharp— edged delta wings. To use this curve for another 
wing shape, an effective aspect ratio computed as 

zo = 4 

tan a L e 

should be used. As described in reference 1, care must be taken when 
using this correlation curve for anything other than sharp-edged wings. 

A rounded leading edge can cause the suction factor, K*_, to go to zero. 

V LE 

AKVSW is the side-edge suction factor for wings with tips. The side- 
edge suction lift is usually small compared to the potential lift except 
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for very low-aspect-ratio wings. Comparisons with subsonic rectangular 
wing data indicate that this factor should be unity, but the magnitude 
for supersonic flow is not known because of a lack of appropriate data 
for correlation purposes. Unless other information is available to the 
user, it is suggested that a value of 1.0 be used for all wings with non- 
zero* tip chords. When there is no wing tip chord, such as on a delta 
wing, AKVSW should be identically zero. 

The next four variables in Item 11, ALPIT, AKVLTl, AKVLT2, and AKVST 
are the corresponding tail parameters. They fall under the same rules 
and guidelines set up for the respective wing parameters. If no tail is 
present, all four values should be set equal to zero. 

The final quantity on this card is WLIMIT, the maximum allowable 
vortex- induced velocity nondimens ionalized by the free-stream velocity. 

The purpose of this variable is to limit the magnitude of the vortex- 
induced velocities on the wing or tail. In the course of program develop- 
ment, a canard-body-wing configuration exhibited a discontinuous lift curve 
at some angle of attack around 8°. Close investigation showed that this 
occurred when the canard trailing vortex passed very near to some wing 
control points, inducing inordinately large velocities, which resulted 
in large local angles of attack. The lifting-surface method reacted 
accordingly and large loading gradients were predicted. This, of course, 
is an unrealistic situation because a true viscous vortex does not behave 
as a potential vortex and induce infinite velocities at its center. For 
this reason a limit was introduced which arbitrarily sets any vortex- 
induced velocity greater than WLIMIT, equal to WLIMIT. 

Generally, WLIMIT should be set equal to 1.0. If , in the process 
of running the program, unusual variations in the lift or pitching moment 
with angle of attack occur which can be attributed to unrealistic vortex- 
induced interference, WLIMIT can be used to limit the magnitude of the 
large induced velocities causing the problem. A value of WLIMIT =0.1 has 
been used in some specific examples to reduce the apparent discontinuity 
in the predicted lift and moment curves and resulted in good agreement 
with experiment. This discussion is not meant to suggest that an arbi- 
trary velocity limit will cure the problems with the near flow fields of 
potential vortices. It is simply included to note that a simple, approxi- 
mate fix is available. If WLIMIT = 0.0, the effect of the free vortices 
on the lifting surfaces is completely eliminated. 
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The above discussion includes all the input required for a typical 
run. The sample cases in the following section cover most of the options 
available in the program. 

Sample cases.- Listings of the input data decks for three sample cases 
are presented in figure 4 and sketches of the configurations chosen for 
these sample cases are shown in figures 5 and 6. These configurations were 
used in the tests of references 14 and 15, respectively. 

Sample case 1 is the canard-body-wing combination shown in figure 5 
at a Mach number of 2.01. No optional output is requested and the canard 
and wing are each defined by a single region. A complete range of angles 
of attack for one canard deflection angle is specified. This input deck, 
shown in figure 4(a), requires approximately 300 seconds running time on 
the IBM 360/67. 

Sample case 2, shown in figure 4(b), is the same configuration with 
the canard removed. No optional output is requested. For demonstration 
purposes, the wing is specified by two regions. The inboard and outboard 
regions are given the same panel spacing as that used in sample case 1. 

This is not normally the situation in laying out this type of lifting 
surface in multiple regions. Generally, some change in spacing between 
the regions is necessary. Of course, a break in the sweep angle of the 
leading or trailing edge dictates the use of multiple regions and the 
panel spacing may or may not be changed. 

Sample case 3 is the wing-body-tail combination shown in figure 6. 
Minimum output is requested and the wing and tail surfaces are each 
defined by a single region. A complete range in angle of attack for three 
tail deflection angles is specified. This input deck, shown in figure 4(c), 
requires approximately 540 seconds on the IBM 360/67. 

Description of Output 

The output produced by the SUPSON computer program for sample case 1 
is shown in figure 7. The first page of output from the program, fig- 
ure 7(a), is a tabulation of most of the input data in Items 1, 2, and 3 
of figure 1. The next page of output, figure 7(b), is a listing of the 
namelist INPUT which describes the geometry of the wing. If the wing is 
divided into more than one region, the variables in namelist INPUT des- 
cribing each region are listed here. The next page of output, figure 7(c), 


157 



is a similar listing of the namelist INPUT describing the geometry of the 
tail. The first entry in the namelist distinguishes between the wing 
(NSF = 1) or the tail (NSF = 2) . 

The following page of output, figure 7(d), summarizes the geometry 
of the configuration by component. The first line at the top of the page 
is the last identification card from the heading in Item 2 of figure 1. 
The wing quantities which are tabulated are: 

OUTPUT ALGEBRAIC 


NOTATION 

NOTATION 

XLE 

X LE „ 

root 

XTE 

XTE root 

XLE (TIP) 

^tip 

XTE (TIP) 


CROOT 

c root 

CTIP 

c tip 

B/2 

b/2 

RAVG 

r 

avg 

3 w 

YSEP 

YSEPW (Item 5) 

XHL 

X, 

ZHL 

z “„ 

XCP 

x cp or XWCP ( 


The same quantities are tabulated for the tail surface if one is 
present. The following quantities are listed for the body. 
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OUTPUT 

ALGEBRAIC 

NOTATION 

NOTATION 

THETA 

0 N (Item 3) 

FINENESS 

X L E / r N 

w 

R(BASE) 

r N 

AVERAGE RADIUS WING 

r avg 

^w 

AVERAGE RADIUS TAIL 

r avg e 

CENTER OF MOMENTS X 

x m 

CENTER OF MOMENTS Z 

Z 

m 

DXI 

Ax ^ 


> Item 3 

DXOUT 

- J 

X 

X 

R 

r 

S 

s or s 
w e 

DR/DX 

dr/dx 


This concludes the general geometric description of the configuration. 
This information is output once at the beginning of each case. The follow- 
ing output is dependent on the information input on Item 11; that is, the 
angle of attack, incidence angles, K* factors, and induced velocity limit. 

Figure 7 (e) is the first page of output for each run within the series 
of runs making up sample case 1. The first line summarizes the information 
input in Item 11 as follows: 

OUTPUT ALGEBRAIC 

NOTATION NOTATION 


ALPHA a 

M M 


INCIDENCE WING 


& 


w 
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OUTPUT 

NOTATION 

ALGEBRAIC 

NOTATION 

INCIDENCE TAIL 

6 e 

WING KVLE* 

K v 

LE w 

WING KVSE* 

K* 

VsE w 

TAIL KVLE* 

K 

LE e 

TAIL KVSE* 

e 

W/V LIMIT 

(w/V) max 


The next block of output on this page is a summary of the strength 
and position of the right-hand vortex (if present) of the symmetrical pair 
of vortices shed from the nose of the body. 


OUTPUT 

NOTATION 

ALGEBRAIC 

JJOTATION, 

GAM/2*PI*V*RB 

r B 

2lrVr N 

Y/RB 

y B /r N 

Z/RB 

Z B //r N 

XS/RB 

V r N 

GAM/2*PI*V 

r B 

2ttV 

Y 

y b 

Z 

Z B 


The following block of output is the induced Beskin upwash at the 
wing control points (x,y,z) due to the presence of the body. The induced 
velocities, V/V(INF) and W/V(INF), expressed as a fraction of the free- 
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stream velocity, are positive in the positive y and z directions, 
respectively. The next block of data are the induced velocities induced 
at the same wing control points by the vortex pair shed from the nose of 
the body. These velocities have the same positive sense as the body- 
induced upwash above. The final block of data in figure 7 (e) is the total 
induced velocity at each control point. 

The next page of output, figure 7(f), contains the results from the 
lifting-surface calculations for the wing. The first several lines are 
reiteration of some input quantities. The remainder of the results are 
calculated quantities associated with the potential forces and moments on 
the wing. The potential span-loading distribution is followed by the span- 
wise position and strength of the trailing vortex. The above output in 
figure 7(f) is printed in subroutine FORCE. The definitions of this output 
are as follows. 


OUTPUT 

NOTATION 

ALGEBRAIC 

NOTATION 

CNP 

C M 

N W(B) , p 

DCN/DALPHA 

dc N /da 

CMP 

c 

^f(B) ,p 

DCM/DCN 

dC n/ dC N 

XAC 

x ac 

XCP 

X CP 

YCP 

y C P 

B 

b 

CCN/2B 

cc 

n 

2b 

GAM/V2PI 

P t 

2irV 
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The next block of output on this page is the spanwise distribution of 
the leading-edge suction. If the wing has a nonzero tip chord, a second 
block of output containing the distribution of side-edge suction is printed. 
This latter output is printed in subroutine EDGFRC and is defined as 
follows. 


OUTPUT 

NOTATION 

ALGEBRAIC 

NOTATION 


Y/(B/2) 

JL. 

b/2 


CCS/2B 

cc s 

2b 


EPS (DEG) 

€ 


CCX/2B 



CCY/2B 



CYC/2B 


(This variable does not 
appear in the sample case 
because c^ = 0.) 


Figure 7 (g) is headed by a summary of the strengths and positions of 
the vortices shed from the configuration ahead of the wing trailing edge. 
The pairs of vortices are listed in the following order. Vortex 1 is the 
right-side body vortex shed from the nose. Vortex 2 is the trailing vortex 
shed from the wing. Vortex 3 is the leading-edge separation vortex shed 
from the wing. If more than one separation vortex is requested, vortex 3 
is the vortex associated with the inboard region and vortex 4 is shed from 
the outboard region. If a vortex is missing for any reason, all following 
vortices are moved up in the table. For example, if no vortices are shed 
by the nose, vortex 1 becomes the trailing vortex shed by the wing, and so 
on. The remainder of figure 7 (g) indicates the induced velocities at the 
tail control points. These velocities are analogous to the induced 
velocities on the wing shown in figure 7(e). 

Figure 7(h) contains calculated results for the tail surface. All 
the quantities on this page are analogous to those described for the wing 
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in figure 7(f). The last entry on this figure is a summary of the strengths 
and positions of all the vortices in the field just aft of the tail trailing 
edge. The first group of vortices are the same as described in connection 
with figure 7(g). The second group of vortices are defined as follows. 
Vortex 4 is the trailing vortex corresponding to the potential lift on 
the tail. Vortex 5 is the leading-edge separation vortex shed from the 
tail. If multiple vortices are shed from the tail leading edge, this 
vortex would be shed from the inboard tail region and vortex 6 would be 
shed from the outboard tail region. 

The next page of output, figure 7(i), is a summary page of the force 
coefficients, pitching-moment coefficients, and centers of pressure on each 
component of the configuration and of the total configuration. The coeffi- 
cients for the individual components are described in Table I. The total 
configuration variables are defined as follows. 


OUTPUT 

ALGEBRAIC 

NOTATION 

_NOTATION_ 

CN 

c = JL 
N qS 

CM 

c = JL- 

m qSi 

XCP 

_ C 

x__ = x - •— £ 
CP m C. T 

N 

CL 

C L = qS 

CDI 

C_ = C_ tan a 
D i L 

CA 

*1$ 

II 

a* 

CDI/CL**2 

c /c 2 
d/^l 


The last page of output for this run, figure 7(j), contains a summary 
of the trajectories of the shed vortices. At the top of the page the 
vortices are identified and their strengths listed. This is followed by 
blocks of output, one block for each x station, describing the local 
crossflow geometry of the configuration and the position of the right-side 
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vortices. Each block of results is separated by approximately DXOUT. 

Notice that the trajectory calculation starts at the wing leading edge 
with a pair of body vortices. As the calculation moves downstream, other 
vortices are shed and added to the calculation. The trajectory calculation 
is carried downstream to a point aft of the tail trailing edge. The vari- 
ables in each block are defined as follows. 


OUTPUT 

ALGEBRAIC 

NOTATION 

NOTATION 

X 

X 

DX 

Ax 

A 

a 

s 

s or s 
w e 

RO 

r o 

DA/DX 

da/d x 

SIGMA (REAL) 

y 

SIGMA (IMAG) 

z 


This completes the output for one card in Item 11 of the input deck. 
Additional runs will repeat the output of figures 7(e) through ( j) . The 
above set of output obtained with N PRINT =0 is a considerable amount 
of output for production runs; therefore, an optional set of output can 
be obtained by setting NPRINT = -2. In this case, the complete output 
consists of figures 7(a), (b) , (c) , (d) , and (i) with some shed vortex 

positions and strengths added. 

Some extra output over and above that shown in figure 7 can be 
obtained when NPRINT > 0. This additional output is useful only for 
diagnostic purposes and is not described herein. This output is labeled 
and the user should have no trouble interpreting the results. 

Program Listing 

The SUPSON computer program consists of the main program, nine func- 
tion subprograms, and twenty-five subroutine subprograms. Each source 
deck is identified in columns 73 through 80 by a four-character identifi- 
cation and a three-digit number sequencing the cards within that deck. 
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The program listing is given on the following pages. The table below will 
act as a table of contents for the listing. 


PROGRAM 

IDENTIFICATION 

PAGE NO. 

MAIN 

SP01 

167 

Subroutines: 



S SWING 

SP02 

172 

TABLE 

SP03 

175 

FORCE 

SP04 

176 

EDGFRC 

SP05 

177 

MATRIX 

SP06 

179 

CNVTX 

SP07 

181 

BDYVTX 

SP08 

181 

CNVNZ 

SP09 

182 

SHAPE 

SP10 

182 

FILL 

SP11 

182 

ZSECT 

SP12 

182 

TCOMP 

SP13 

183 

TAINT 

SP14 

183 

LINEQS 

SP15 

184 

SOLVE 

SP16 

184 

COMP 

SP17 

184 

EXTV EL 

SP18 

185 

TRJTKY 

SP19 

185 

FCT 

SP20 

185 

HPCG 

SP21 

186 

OUTP 

SP22 

188 

ZVTX 

SP23 

188 

WPANL 

SP24 

188 
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PROGRAM 


IDENTIFICATION 


PAGE NO 


KFACT 

CH1416 

Functions : 

EQ14 

EQ21 

EQ24 

EQ24L 

EQ26 

EQ30 

EQ30L 

EQ31 

CHRT8 


SP25 

SP26 

SP27 

SP28 

SP29 

SP30 

SP31 

SP32 

SP33 

SP34 

SP35 


190 

190 

192 

192 

192 

192 

192 

192 

193 
193 
193 
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O330l**«« 

mauMiiii 


k. z a x > 


OOOOOO 3 

Okkkkk O 


►k * 

C On 

» O K» 


^ounloH * 

K «KK «K 3 
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c - a* o-^«»i»ar-a* 




^ r .KNr>.NKC«ac«««aaa»»a»»»»»»»oooooooooet«*x>«<*'« 


eaoco«a«««««a««ea«oa«««««a«««oaa*«a«a«««a««««« 4 «« 4 *««««««««« 

\ma'\.M'M«»nifliniivniNivni(V<M««niM«ni(«MiwKi«NivnjM!vivmivniniN'vniniMiM)wniivniwNAiMNfliMM«ii«iiiMM 

nwnnnnanMannaanMKMaaNananaannNnAnanaaaanaaaaaaaaaaaaaaaaaa** 



o<«r««t*ii> 4 K«»o^niia*iA«Ka» o^mh** 4 n*#o- 35 *| 5 «n«*;h 55 *£«n 

K«o««a*«a««*»»»»»»»»»oeooeoooooM<«« *••*»» 

a x> « « a « « ««««*« «« « • a a « « aaaaa a aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

<VrVftiAiAi'VMnnVN)«MNaNNNN(«NNaNNNNNNNNNNNWaN«Naa«NWM««NaNM««iaN*NNMI««NnNa«aNNNMNMMNNNN 

an nnan «nnaa««n«aaaaa«annaa«aa«aaa««aaaaa«aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa»ana« 
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o-oirt » l(» 

3 0-50 300000000^^00 
O 0900900009000 ooo 

Lai&xa.ai&ifta&&icka 

MMMOnwnnNMMMWftAMi* 






OOOOOOOOO OO OO ooooooooo oooo 




* «c 

Ui 
9 *- 

o« 

« 

0 o 

*o 

K 

aw 

w*~ 

w 

55 

K 9 
9 9 • 
9 >- 

1 oo 


o o 

►- t- 

ss 


saw 
B« N 

w • UI 


u > w 

K 

w »<• 


r r 1 


» O w at 

> L>(IU 

a k coo 

«. l» 9 J® 

CO * < « 


vw« > w« aa^^KHMOooooo 

k.k.ocaot<auoikiuaijaw9iiiz 

owswKwatw 


OOOOOOOOOww — — 

ooooooo oooooooooo 
W^WIWWf^W^IWI^MWWIVIMWWIIJI 

aaiassasassssssss 

nanM«<i««0s«<> M<, * AN 






toooaoooo^xw^w^i 


► ooo oo o « 
INNNMMMI 


►ooooooooo 

►000000990000 


4 a 

X 2 
*- 9 
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o .? a » ao aoaoo’aooooooDoaoooooooooooooDooooooooooooo 






uuu u 








ooooooo ooooooo oooo 
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Figure 1.- Input format for SUPSON program, 







Item 7 if NALPW 



Figure 1.- Concluded- 












0 


e so j <■ \ i u \ \ * i »* 

SAMP^t CASf- 1 kfcF, UASA 111 X-32B StPT # I960 C, DKlVt« 
CANAPU • ^InG - HHUt LU fc( 6 1 A T I UN 


2,01 

o. 

192,0 

l . 

15,55 

25.0 

10, 

,297 

• »* i 6 




*956 

.25* 





I. 9ab ,^u5 

2.60b .b/5 

3.267 ,f>H2 

3,929 , /O 

9,27 5 ,82b 

b,?SS ,9a 

b.92 .90* 

6,b85 l.oa* 

II, 375 1,51 

1 7 , 7a9 1,667 

17. 7S 1.667 

17,6 1,667 

16,61 1,667 

2b. 1,667 

50, 1,667 

57, 1,667 

56, 1,667 

1,39 3,37 9,12b 0. 0, 

♦INPUT 
PEn» C , 9b * 
run S s a, 

C'JlS* a# 

RUUTLfcs a, 27b, 

RUL) T T E * 1 1 ,37b, 

rout v * 1 , 5 a, 

TiPLfr. *11,375# 

TIPUM1.5/S, 

TIPV* 3 , 37, 

TC.RUU1* 0,0# 

TC T IPs u,0, 

♦ fcND 

1,667 b , 0 1 55 25, 0, 0, 

♦Input 
PLRs 0,95, 

RUMS* 6, 

CUL8»10, 

MUUT L 6*16,61, 

MUUTTt s37,0, 

RUUTYS 1,667, 

TIPUK«3/,0, 

TlPU*i7,0, 

TlPYs 6,36, 

TUROUTe 0,0, 

TCTIP* 0,0, 

♦ t NO 

1 00,0 10, ,5 u« 0, 0, • * 0, 0, 

1 02,0 10, ,5 0, 0, 0, 0, 0. 

i oa,o to, ,5 0 . 0 , 0, ,5 0, 0, 

1 06,0 10, ,5 Of 0, 0, ,5 0. 0, 

1 06,0 10, ,5 0, 0. 0, #5 0, 0, 

1 10,0 10, ,5 0, 0, 0, .5 0, 0, 

112,010, ,5 Of 0, 0, ,5 0, 0, 

1 16,0 10, ,5 0, 0, 0, 0, 0, 

1 20,0 10, ,5 0, 0, 0, ,5 0, 0, 

1 20.0 10, ,G 0, 0, 0, ,0 0, 0, 


(a) Sample case 1. 

Figure 4.- Sample input decks for SUPSON program. 
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* t 'v 

i) 

(i 1 f* ^ 

10 0 0 

0 0 

SAMPLE 

tA$t * 

pwi • (,* A SuPSt'* 



* I Ml 

- HOo > 

Li M 1 Gt'* a t i L ik 

F , '5ASA T m *•!*& 

C, OWlVtR 

*,01 

19?. r 

lS.35 *5,0 

10. 1, 

1, 

9. 

ft. 




.*97 

• o 7 to 




.9bto 

.23* 




l,9ub 

, P«.b 




*,6Ub 

,b73 




i,*67 

, toH? 




3.9*9 

• 7 8 




9. *75 

,b*b 




*.*55 

.9« 




*.9* 

,996 




6,5M 

1 , Uto* 




U.37S 

t .31 




17, 7«9 

1 ,66 i 




17.7b 

\ .667 




1 7 , to 

1 , dto / 




1 to , 6 1 

1 , * 6 7 




**• 

1 , Ob 7 




30. 

1 .667 




37. 

1,667 




5tt. 

1,667 




1 ,6b7 

ft, 36 

25. 

0, 



♦ J NPl'T 
KU^SSfc# 

C^LSbS, 

St L I s i , 

KttUTL t S- i 6,6 1 0 
NUUT Ifc •.W, (U 
RUUTtsi ,6b/, 

TlPUs<>7.fto5, 

TiHU«i7 f n, 

Tl^V«S # CliS, 

TLhOtitsO,'/. 

TCII^bO.O, 

ttKiO 

♦INPUT 

WUwSs6, 

tOLSES, 

St l 1 x ) , 

«UUTLts?7.ft0b, 

WUUtUsiJ.O, 

«UU1 Yxb.Ol ib, 

UPlfcs57 # P # 

TXPTK* 47*0# 

TlPy«ft,3h, 

TtMulaO .0# 

TU JHbO.O# 

♦ t NO 

1 1 1> , 0, ,5 0.0 0,0 0,0 0,0 0,0 0,0 


(b) Sample case 2. 


Figure 4.- Continued 



? t 4 -2 

SA-PlF CASE J 
kjNf, • hi‘Oy - 


1,30 

34Q, 

o. 

.0 


,984 

10, 

1,54 

lb. 

1.^3 

20,29 

2.2 

25.16 

2.34 

29,75 

2.36 

31.11 

2.376 

34.26 

2.34 

38, 

2.24 

40,65 

2.13 

«2,93 

2.05 

9b, 11 

1.62 

97. 

1,75 

2,36 

16.42 

♦Input 


PfcKxG,95, 


R0*S=6 , 


CUIS=1 0, 


ROUTLF S20 

.29, 

R0UTTE834 

,26. 

RUOT Y *2 , 38 , 

TIPLES25, 

It*# 

TIPTU31 , 

11# 


n^yath.a?, 

TCNDUTbU.O* 

TCl^sC.n, 

♦ tNU 

2, 6,34 

♦Input 
PtR«0,95# 
kU*S«4 , 

COiSzS, 

RUCjTie 840,85, 
RUUTTEsab.il, 
ROUT Y*2.0, 
TlPlfesai ,8, 
TlPU«43,89, 
TlPy 86,34, 
TCRUUT *0 1 0 , 
TCTIPeO.O, 

♦ Enq 

1 2. 0. 

l a. o. 

1 6, 0, 

1 &. 0, 

1 10, 0. 

1 12, 0. 

t IP. 0, 

1 2, 0. 

1 a, o, 

1 p, o, 

1 «. 0, 

l 10, 0 , 

1 12. 0, 

1 IP, 0, 

l 2, o. 

1 4, 0, 

1 p , 0 1 

i a. 0, 

i io. f), 

1 12, 0, 

1 IP. 0, 


'j l 1 0 1 1 0 

REF, NASA f-fMl' 6-11-59A PAY 1959 
TAIL Cl’* F lOUkAT TUN 

11,32 23.59 12. 


23,59 0, 0, 


42,4* 0, 0, 


0. 0, 0. 0, 

0 . 0, 0, 0 , 

0, 0, 0, 0. 

0, 0, 0, 0, 

0. 0, 0. 0, 

0, 0, 0. 0. 

0, 0. 0, 0, 

0. 0, 0, *10, 

0, 0, 0, *10. 

0, 0, 0. *10. 

0, 0, 0, *10. 

0, 0, 0, -10, 

0, 0, 0. -10, 

0, 0, 0, -10, 

0. 0, 0, -20, 

0, 0. 0, *20, 

0, 0, 0 , -20, 

0, 0, 0. *20, 

0. 0, 0. *20, 

o, 0 , o, -20, 

U, 0, 0, -20, 


(c) Sample case 3. 
Figure 4.- Concluded. 


1 0 0 
PETERSON AND HfcNfcES 

1. 1. 


0, 0, 0, 

0. 0. 0, 

0, 0 , 0 , 

0, 0. 0, 

0. 0. 0, 

o. 0, 0, 

0 , 0 , 0 , 

0, 0, 0, 

0. 0, 0, 

0 , 0 , 0 , 

0. 0, 0, 

0, 0, 0. 

0, 0. 0, 

0, 0, 0, 

0, 0, 0, 

0, 0. 0, 

o, 0, 0, 

0, 0, 0, 

0, 0, 0, 

0, 0. 0, 

0, 0, 0, 
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46.90 


Figure 6.- Wing-body-tail configuration 
for sample case 3. 




EINPUTM 
NHE A0= 2 
NTDL= 20 
NPP I NT= 0 
MPR INT= 0 
NOSE V= t 
NP FGNW= 1 
NALPW= 0 
NSFPW= 1 
NP FGNT = 1 
NALPT = 0 
NSFPT= l 
NRODY= 0 
NAFT = 0 
CFND 

SAMPLE CASE 1 REF. NASA TM X-328 SEPT. I960 C. DRIVER 
CANARD - MING - BODY CONFIGURATION 

M REFS PFFL XCM OXOUT DXI 

2.01000 192.00000 15.33000 25.00000 l. 00000 1.00000 


(a) Page 1. 

Figure 7.- Output from SUPSON program 
for sample case 1. 
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£!NPUT 
NSF= 1 
ROWS= 4 
rOLS= 4 
P EP = 0.950 
ROOTLE= 4.2750 
ROOTTE= 11.3750 
ROOTY= 1.339999 
T I PLE* 11.3750 
TI PTF= 11.3750 
TIPY= 3.370 
TCROOT = 0.0 
TCTIP= 0.0 
ASYM= F 
SFCT= 1 
ALT= -1.0 
SC AL E= 1.0 
TYPE= 4 
SEND 


(b) Page 2. 
Figure 7.- Continued 



EINPUT 
NS F= ? 

ROWS= 6 
cnLS= 10 
PFR= 0.950 
R00TLE= 18.60999 
ROOTTF= 37.0 
R0PTY= 1.6670 
TIPLE= 37.0 
TTPTE= 37.0 
T I PY= 8.360 
TCPnOT= 0.0 
TC T I P= 0.0 
ASYM= F 
S ECT = 1 
ALT= -1.0 
SC AL E= l.o 
TYPE= 4 
F.EN0 


(c) Page 3. 
Figure 7.- Continued. 



CANAPD - WING - RODY CONFIGURATION 


XLF 

4.27500 

XTE 

11.17500 

XLEITIP) 

11.17500 

XTEITIP) 

11.37500 

CROOT 

7.10000 

CT IP 
0.00000 

B/2 

3.37000 

PAVG 

1.34000 

YSEP 

3.37000 

XHL 

9.12500 

ZHL 

0.00000 

XCP 

0.00000 



XLE 

18* 6099*5 

XTE 

37.00000 

XLEITIP) 

37.00000 

XTEITIP) 

37.00000 

CROOT 

18.39001 

CTIP 

0.00000 

B/2 

8.36000 

PAVG 

1.66700 

YSEP 

8.36000 

XHL 

25.00000 

ZHL 

0.00000 

XCP 

0.00000 




BODY 


nose 

THETA FINENESS R C BASE I 
10.000 2.508 0.826 


AVERAGE RADIUS CENTER OF MOMENTS 
WING TAIL X Z OX I 

1.340 1.667 25.000 0.000 1.000 


TABLF OF BODY COORDINATES 



X 

R 

S 

OR/OX 

l 

0.00000 

0.00000 

0.00000 

0.25589 

2 

0.29700 

0.07600 

0.07600 

0.24707 

3 

0.95600 

0.23300 

0.23300 

0.22630 

4 

1.94500 

0.44500 

0.44500 

0.20415 

5 

2.60500 

0.57300 

0.57300 

0.17930 

6 

3.26700 

0.68200 

0.68200 

0.15634 

7 

3.92900 

0.70000 

0.78000 

0.14049 

8 

4.27500 

0.82600 

0.82600 

0.11690 

9 

6.05000 

1.00502 

1.84750 

0.08041 

10 

7.82500 

1.11146 

2.35500 

0.05795 

11 

9.60000 

1.21073 

2.86250 

0.05593 

12 

11.37500 

1.31000 

3.37000 

0.05561 

13 

11.37600 

1.31005 

1.31005 

0.05565 

14 

17.74899 

1.66700 

1.66700 

0.02800 

15 

17.75000 

1.66700 

1.66700 

0.00000 

16 

17.79999 

1.66700 

1.66700 

0.00000 

17 

18.60999 

1.66700 

1.66700 

0.00000 

18 

20.44897 

1.66700 

2. *3630 

0.00000 

19 

22.28798 

1.66700 

3.00560 

0.00000 

20 

24.12698 

1.66700 

3.67490 

0.00000 

21 

25.96597 

1.66700 

4.34420 

0.00000 

22 

27.80498 

1.66700 

5.01350 

0.00000 

23 

29.64397 

1.66700 

5.68280 

0.00000 

24 

31.48297 

1.66700 

6.35209 

0.00000 

25 

33.3219R 

1.66700 

7.02139 

0.00000 

26 

35.16096 

1.66700 

7.69069 

0.00000 

27 

36.99998 

1.66700 

8.36000 

0.00000 

28 

17.00099 

1.66700 

1.66700 

0.00000 

29 

38.00000 

1.66700 

1.66700 

0.00000 


(d) Page 4. 
Figure 7.- Continued. 


DXOUT 

1.000 


! 
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ALPHA 

16.00 


-- INCIDENCE — 
M WING TAIL 

2.01 10.000 0.000 


WING TAIL 

KVLE* KVSE* KVLE* KVSE* W/V 

0.500 0.000 0.500 0.000 l 


NOSE VORTICES - GAM/?*Pl *V*RB Y/RR 
0.0906 0.646 


Z/RB XS/RB 

1.311 0.000 


SUMMARY OF VOPTFX STRENGTHS AND POSITION AT X * 4.275 


VORTFX GAMMA /? *0 1 *V Y 

l 0.074799 0.533 


Z 

1.083 


VELOCITY INDUCED AT SPECIFIED FIELD POINTS ON WING {BESKIN) 



X 

Y 

Z 

V/VUNF) 

W/VCINE) 

1 

6.606 

1.582 

0.437 

5.7449F-02 

9.5909E-0? 

2 

8. 169 

1.582 

0.166 

2.9297E-02 

1 • 3809E-01 

3 

9.733 

1.582 

-0.106 

-2. 1923E-02 

1.6345E-01 

4 

l 1.297 

1.58? 

-0.377 

-8. 1 243E-02 

l • 6068F-01 

5 

7.946 

2.084 

0.205 

1 • 5486E-02 

7. 8089E-02 

6 

9.070 

2.084 

0.009 

8. 1490E-04 

8. 9664E-02 

7 

10.195 

2.084 

-0.186 

- t • 7450E-02 

9.7132E-02 

8 

11.319 

2.084 

-0.381 

-3 • 7576E-02 

9.93636-02 

9 

9.270 

2.581 

-0.025 

-1. 1 602E-03 

5.9592E-02 

1 0 

9.960 

2.581 

-0.145 

-7.0939E-03 

6.2933E-02 

1 1 

10.650 

2.581 

-0.265 

-1.3584E-02 

6.5480E-02 

12 

1 1 .340 

2.58 l 

-0.385 

-2.04 74E-02 

6.7140E-02 

13 

10.473 

3.032 

-0.234 

-7 .3537E-03 

4. 7348E-02 

14 

10.769 

3.032 

-0.285 

-9. 1522E-03 

4.81 80E-02 

15 

11.064 

3.032 

-0.337 

-1 . 1005E-02 

4. 8926E-02 

16 

11.360 

3.03? 

-0.388 

-1.2907E-02 

4.9582E-02 



VELOCITY 

INDUCED AT 

SPECIFIED 

FIELD POINTS ON WING (VORTEX) 


X 

Y 

Z 

V/VUNF) W/V ( INF ) 

1 

6.606 

1.582 

0.437 

1 . 6056E-02 -5.3901E-03 

? 

8. 169 

1.58? 

0.166 

7.68 2 3E— 03 -7.9791E-03 

3 

9.733 

1.58? 

-0.106 

2. 7704F-03 -6.3888E-03 

4 

11.297 

1.58? 

-0.377 

6. 1 840E-04 -4. 2710E-03 

5 

7.946 

2.004 

0.205 

7. 2 1 36E-03 -2. 1 890E-03 

6 

9.070 

2.084 

0.009 

5. 052 IE-03 -3. 0338E-03 

7 

10.195 

2.084 

-0.186 

3. 3 1 62E-03 -3.091 IE-03 

0 

11.319 

2.084 

-0.381 

2.0716 E-03 -2. 7641E-03 

9 

9.270 

2.581 

-0.025 

4.01 01E-03 -9.7394E-04 

10 

9.960 

2.581 

-0.145 

3 . 3966E-03 -1.2613E-03 

11 

10.650 

2.581 

-0.265 

?• 8322E-03 -1.4172E-03 

1 2 

n.340 

2.581 

-0.385 

2. 3305E-03 -1.4747E-03 

13 

10.473 

3.032 

-0.234 

2.5851E-03 -4. 9036E-04 

14 

10.769 

3.03? 

-0.285 

2 • 4 36 7E— 03 -5.6838E-04 

15 

11.064 

3.032 

-0.337 

2 • 2920E-03 -6.3389E-04 

16 

11.360 

3.03? 

-0.388 

2.151 1 E-03 -6. 88 15E-04 


TOTAL INDUCED VELOCITY AT SPECIFIED FIELD POINTS ON WING 



X 

Y 

Z 

V/VUNF) 

W/V( INF) 

1 

6.606 

1.582 

0.437 

1 • 6056E-02 

9.051 9E-02 

2 

8.169 

1.582 

0. 166 

7.6823E-03 

1.301 1 E-Ol 

3 

9.733 

1.58? 

-0.106 

2.7704E-03 

1 • 5706E— 01 

4 

11.297 

1.582 

-0.377 

6. 1840E-04 

1.5641E-01 

5 

7.946 

2.084 

0.205 

7.2 136E-03 

7.5900E-02 

6 

9.070 

2.084 

0.009 

5.0521E-03 

8. 6630E-02 

7 

10.195 

2.004 

-0.186 

3.3162E-03 

9. 4041E-02 

8 

1 1.319 

2.084 

-0.381 

2 ♦ 0716E-03 

9. 6599E-02 

9 

9.270 

2.581 

-0.025 

4.0 I 01E-03 

5.86 18E-02 

10 

9.960 

2.581 

-0.145 

3 . 3966E-03 

6. 1672E-02 

l l 

10.650 

2.581 

-0.265 

2.8322E-03 

6.4063E-02 

12 

11.340 

2.581 

-0.385 

2. 3305E-03 

6. 5665E-02 

13 

10.473 

3.03? 

-0.234 

2.585 IE-03 

4. 6858E-02 

14 

10.769 

3.032 

-0.285 

2.4367E-03 

4.7611E-02 

15 

11.064 

3.032 

-0.337 

2 • 2920E-03 

4.8292E-02 

16 

11.360 

3.032 

-0.388 

2.1511 E-03 

4.8894E-02 





(e) Page 

5. 


Figure 7.- Continued. 


LIMIT 

000 
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NASA/AMFS wing program 


♦****WING***** 


CANARD - WING 

flow conditions - 


- BODY CONF I GURA T ION 

MACH BETA 

2.01000 1.74359 


PFFFRFNCF QUANTITIES - AREA LFNGTH XMOM 

192.000 15.330 25.000 


B/2 PANELS 

3.370 16 


nCN/D ALPHA 

CNP (RAD**-1I ( DEG**— 1 ) * BETA CMP 

0.06369 0.11375 0.00199 0.19B33 0.06657 


DCM/DCN 

1.04514 


AN L P AO DISTPIBUTION 
Y CCN/7R 

1 . 6fl 167 0.30599 

2 • 0843 1 0.26667 

?. 58055 0.21207 

3.03166 0.10R9R 


?B * 13.4R00 


XAC 

8.978 


TRAILING VORTEX .... Y = 2.B223 GAM/V2PI = 3.2823E-01 


XCP 

8.977 


SP 

REGION COL 

l l 

1 2 

l 3 

l 4 


ANWISE DISTRIBUTION 


Y/ (B/2) 
0.47292 
0.62352 
0.77411 
0.92470 


CCS/2B 

0.07600 

0.14431 

0.22702 

0.27376 


OF L.E. SUCTION 


FPS(DEG) 

90.460 

89.745 

38.989 

90.778 


CCX/2B 

-0.02170 

-0.03905 

-0.05855 

-0.08027 


CCY/2B 
0.07366 
0. 1389? 
0.21934 
0.26696 


(f) Page 6. 
Figure 7.- Continued. 


YCP 

2.145 
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SUMMARY HF VOPTEX STRENGTHS AND POSITION AT X - 11.776 


VORTEX GAMM A/2*P I *V Y 

1 0.074799 0.373 

2 0.328235 2.82? 

3 0. 1*50326 2.595 


Z 

1.968 

-0.391 

0.360 


VELOCITY INDUCED AT SPECIFIED FIELD POINTS ON 



X 

Y 

Z 

V/V (INF 

l 

22.282 

1.996 

0.000 

0.0000 

2 

25.196 

1.996 

0.000 

0.0000 

3 

28.111 

1.996 

0.000 

0.0000 

4 

31.025 

1 .996 

0.000 

0.0000 

5 

33.940 

1 .996 

0.000 

0.0000 

6 

36.854 

1.996 

0.000 

0.0000 

7 

23.828 

2.664 

0.000 

0.0000 

8 

26.437 

2.664 

0.000 

0.0000 

9 

29.045 

2.664 

0.000 

0.0000 

10 

31.653 

2.664 

0.000 

0.0000 

11 

34.261 

2.664 

0.000 

0.0000 

12 

36.870 

2.664 

0.000 

0.0000 

13 

25.374 

3.333 

0.000 

0.0000 

14 

27.676 

3.333 

0.000 

0.0000 

15 

29.978 

3.333 

0.000 

0.0000 

16 

32.281 

3.333 

0.000 

0.0000 

17 

34.583 

3.333 

0.000 

0.0000 

18 

36.885 

3.333 

0.000 

0.0000 

19 

26.919 

4.001 

0.000 

0.0000 

20 

28.915 

4.001 

0.000 

0.0000 

21 

30.912 

4.001 

0.000 

0.0000 

2? 

32.908 

4.001 

0.000 

0.0000 

23 

34.904 

4.001 

0.000 

0.0000 

24 

36.900 

4.001 

0.000 

0.0000 

25 

28.463 

4.669 

0.000 

0.0000 

26 

30. t54 

4.669 

0.000 

0.0000 

27 

31.844 

4.669 

0.000 

0.0000 

28 

33.535 

4.669 

0.000 

0.0000 

29 

35.225 

4.669 

0.000 

0,0000 

30 

36.915 

4.669 

0.000 

0.0000 

31 

30.006 

5.336 

0.000 

0.0000 

32 

31.391 

5.336 

0.000 

0.0000 

33 

32.776 

5.336 

0.000 

0.0000 

34 

34.161 

5.336 

0.000 

0.0000 

35 

35.546 

5.336 

0.000 

0.0000 

36 

36.931 

5.336 

0.000 

0.0000 

37 

31.546 

6.002 

0.000 

0.0000 

38 

32.626 

6.002 

0. 000 

0.0000 

39 

33.706 

6.002 

0.000 

0.0000 

40 

34.786 

6.002 

0.000 

0.0000 

41 

35.866 

6.002 

0.000 

0.0000 

42 

36.946 

6.002 

0.000 

0.0000 

43 

33.079 

6.664 

0.000 

0.0000 

44 

33.855 

6.664 

0.000 

0.0000 

45 

34.632 

6.664 

0.000 

0.0000 

46 

35.408 

6.664 

0.000 

0.0000 

47 

36.185 

6.664 

0.000 

0.0000 

48 

36.961 

6.664 

0.000 

0.0000 

49 

34.592 

7.319 

0.000 

0.0000 

50 

35.069 

7.319 

0.000 

0.0000 

51 

35.546 

7.319 

0.000 

0.0000 

52 

36.023 

7.319 

0.000 

0.0000 

53 

36.499 

7.319 

0.000 

0.0000 

54 

36.976 

7.319 

0.000 

0.0000 

55 

35.968 

7.914 

0.000 

0.0000 

56 

36. 172 

7.914 

0.000 

0.0000 

57 

36.377 

7.914 

0.000 

0.0000 

58 

36.581 

7.914 

0.000 

0.0000 

59 

36.785 

7.914 

0.000 

0.0000 

60 

36.990 

7.914 

0.000 

0.0000 


(g) Page 7. 

Figure 7.- Continued. 


TAIL ( BESKIN) 

W/ V ( INF ) 
l • 9482E-01 
1. 9482E-01 
l • 9482E-0 1 
l • 9482F— 01 
1.9482E-01 
1.9482E-01 
l • 093 IE— 01 
1.0931E-01 
1 • 0931E-01 
1 • 093 IE-01 
1.0931E-01 
1. 0931E-01 
6.9863E-02 
6.9863 E-02 
6. 9863E-0? 
6.9863E— 02 
6.9863E-02 
6.9863E-02 
4.84 77E-02 
4. 8477E-02 
4. 8477E-02 
4.84 77E— 02 
4.8477E-02 
4.8477E-02 
3.5602E-02 
3.5602E-02 
3.5602E-02 
3.5602E-02 
3. 5602E-02 
3 • 5602E-02 
2.7257E-02 
2.72 57E-02 
2.72 57E-02 
2.7257E-02 
2 • 7257E-02 
2.7257E-02 
2 • l 545E-02 
2.15456-0 2 
2*1 545E-0? 
2* 1545E-02 
2 • 1545E-02 
2. 1545E-02 
l • 7472E-02 
1* 74 72E-02 
l • 7472E-02 
1 • 7472E-02 
1*74 72E-02 
1*74 7ZE-02 
l • 4487E-02 
1*44 87E-02 
1.4487E-02 
1.4487E-02 
1 • 4487E-02 
1.4487E-02 
1 • 2391E-02 
1.2391E-02 
1 • 2391E-02 
1.2391E-02 
1*239 IE-02 
1 • 2391E-02 



VELOCITY INDUCEO AT SPECIFIED FIELD POINTS ON TAIL (VORTEX) 


X 


1 

22.282 

2 

25.196 

3 

28*111 

4 

31.025 

5 

33.940 

6 

36.854 

7 

23.828 

8 

26.437 

9 

29.045 

10 

31.653 

11 

34.261 

12 

36.870 

13 

25.374 

14 

27.676 

15 

29.978 

16 

32.281 

17 

34.583 

18 

36.885 

19 

26.919 

20 

28.915 

21 

30.9X2 

22 

32.908 

23 

34.904 

24 

36.900 

25 

28.463 

26 

30.154 

27 

31.844 

28 

33.535 

29 

35.225 

30 

36.915 

31 

30.006 

32 

31.391 

33 

32.776 

34 

34.161 

35 

35.546 

36 

36.931 

37 

31.546 

38 

32.626 

39 

33.706 

40 

34.786 

41 

35.866 

42 

36.946 

43 

33.079 

44 

33.855 

45 

34.632 

46 

35.408 

47 

36.185 

48 

36.961 

49 

34.592 

50 

35.069 

51 

35.546 

5? 

36.023 

53 

36.499 

54 

36.976 

55 

35.968 

56 

36.172 

57 

36.377 

58 

36.581 

59 

36.785 

60 

36.990 


Y 

1 

1.996 

0.000 

1.996 

0.000 

1.996 

0.000 

1.996 

0.000 

1.996 

0.000 

1.996 

0.000 

2.664 

0.000 

2.664 

0.000 

2.664 

0.000 

2.664 

0.000 

2.664 

0.000 

2.664 

0.000 

3.333 

0.000 

3.333 

0.000 

3.333 

0.000 

3.333 

0.000 

3.333 

0.000 

3.333 

0.000 

4.001 

0.000 

4.001 

0.000 

4.001 

0.000 

4.001 

0.000 

4.001 

0.000 

4.001 

0.000 

4.669 

0.000 

4.669 

0.000 

4.669 

0.000 

4.669 

0.000 

4.669 

0.000 

4.669 

0.000 

5.336 

0.000 

5.336 

0.000 

5.336 

0.000 

5.336 

0.000 

5.336 

0.000 

5.336 

0.000 

6.002 

0.000 

6.002 

0.000 

6.002 

0.000 

6.002 

0.000 

6.002 

0.000 

6.002 

0.000 

6.664 

0.000 

6.664 

0.000 

6.664 

0.000 

6.664 

0.000 

6.664 

0.000 

6.664 

0.000 

7.319 

0.000 

7.319 

0.000 

7.319 

0.000 

7.319 

0.000 

7.319 

0.000 

7.319 

0.000 

7.914 

0.000 

7.914 

0.000 

7.914 

0.000 

7.914 

0.000 

7.914 

0.000 

7.914 

0.000 


V/V ( INF ) W/V < INF ) 
1.175TE-01 -3.8098E-01 
1.3732E-01 -3.1656E-01 
1.3330E-0I — 4.6734E-01 
6.8128E-02 -3.2909E-01 
5.6476E-02 -2.7537E-01 
5. 1282E-02 -2.5169E-01 
1.9245E-01 -1.9361E-01 
1.7936E-01 -1.1332E-01 
4. 7276E-0I -1.7480E-01 
1.4355E-01 -2.5320E-01 
1.1303E-01 -2.0003E-01 
1.0179E-01 -1.7988E-01 
1.6672E-01 -5.9287E-02 
1.5478E-01 -2. 8686E-02 
1.9363E-01 5.7063E-02 

3.2934E-01 -1.9808E-01 
1.6661E-01 -1.7577E-01 
I.3429E-01 -1.4616E-01 


1.2679E-01 

3. 8934E-03 

1.1837E-01 

1.4526E-02 

1.2399E-01 

4.8676E-02 

1.8053E-01 

8. 6594E-02 

2.8443E-01 

2.4177E-02 

2.4060E-01 -8.0302E-02 

9.2526E-02 

2.6850F-02 

8.6790E-02 

3.35 16E-02 

9.0593E-02 

4.7019E-02 

9.5842E-02 

5. 8034E-02 

1.1234E-01 

7.1089E-02 

1.3756E-01 

6. 9675E-02 

6.565 IE-02 

3.4335E-02 

6.5420E-02 

3.9962E— 02 

6.6716E-02 

4. 3205E-02 

6.6593E-02 

4. 6975E-02 

7.0247E-02 

5.2231E-02 

7. 5330E-02 

5.4015E-02 

4.9274E-02 

3.72 14E-02 

4.9898E-02 

3 . 8327E-02 

4.9317E-02 

3.8843E-02 

4.9281E-02 

4.0998E— 02 

5. 0686E-02 

4. 3004E-02 

5.2233E-02 

4. 3356E-02 

3.8066E-02 

3.4316E— 02 

3. 7679E-0? 

3.4392E-02 

3.7439E-02 

3.5232E-02 

3.7757E-02 

3.6245E-02 

3.8319E-02 

3.6853E-02 

3. 8964E-02 

3.681 IE— 02 

2.9216E-02 

3.0988E-02 

2.9235E-02 

3. 1407E-02 

2.93 78E-02 

3.18 I7E-02 

2.9593E-02 

3.2080E-02 

2.98 13E-02 

3.21 63E-02 

3.0I25E-02 

3.2031E-02 

2.3844E-02 

2. 8540E-02 

2.39 16E-02 

2 . 8576E-02 

2. 3987E-02 

2. 8600E-02 

2.4055E-02 

2. 8612E-02 

2.4127E-02 

2.8607E-02 

2. 4282E-02 

2. 8490E-02 


(g) continued. 
Figure 7. — Continued 



TOTAL INDUCED VELOCITY AT SPECIFIED FIELD POINTS IN TAIL 



X 

Y 

Z 

V/VUNFI W/VUNF) 

l 

22.28? 

1.996 

0.000 

1 » 1757E— 01 -1. 86 15E-01 

2 

25. 196 

1.996 

0.000 

1.3732E-01 -1.21 74E-01 

3 

28.111 

1.996 

0.000 

1.3330E-01 -2.7251E-01 

4 

31.025 

1.996 

0.000 

6.8128E-02 -1.3427E-01 

5 

33.940 

1.996 

0.000 

5.6476E-02 -8.0545E-02 

6 

36.854 

1.996 

0.000 

5. 1282E-02 -5. 6861E-02 

7 

23.828 

2.664 

0.000 

1.9245E-01 -8.4293E-02 

B 

26.437 

2.664 

0.000 

1.7936E-01 -4. 0053E-03 

9 

29.045 

2.664 

0.000 

4.7276E-01 -6.5490E-02 

10 

31.653 

2.664 

0.000 

1.4355E-01 -1.4389E-01 

II 

34.261 

2.664 

0.000 

1.1 303E-01 -9.0713E-02 

12 

36.870 

2.664 

0.000 

1.0179E-01 -7.0564E-02 

13 

25.374 

3.333 

0.000 

1.6672E-01 1.0576E-02 

14 

27.676 

3.333 

0.000 

1.5478E-01 4.11 77E-02 

15 

29.978 

3.333 

0.000 

1.9363E-01 1.2693E-01 

16 

32.281 

3.333 

0.000 

3.2934E-01 -1.2822E-01 

17 

34.583 

3.333 

0.000 

1.6661E-01 -1.0590E-01 

18 

36.885 

3.333 

0.000 

1.3429E-01 -7.6296E-02 

19 

26.919 

4.001 

0.000 

1.2679E-01 5.2371E-02 

20 

28.915 

4.001 

0.000 

1.1 837E-01 6. 3003E-02 

21 

30.912 

4.001 

0.000 

1.2399E-01 9.7154E— 02 

22 

32.908 

4.001 

0.000 

1. 8053E-01 l . 3507E-01 

23 

34.904 

4.001 

0.000 

2 . 8443E— 01 7 . 2655E-02 

24 

36.900 

4.001 

0.000 

2.4060E— 01 -3. 1824E-02 

25 

28.463 

4.669 

0.000 

9.2526E-02 6.2452E-02 

26 

30. 154 

4.669 

0.000 

8.6790E— 02 6.9118E-02 

27 

31.844 

4.669 

0.000 

9.0593E-02 8.2621E-02 

28 

33.535 

4.669 

0.000 

9.5842E-02 9.3636E-02 

29 

35.225 

4.669 

0.000 

1.1234E-01 1.0669E-01 

30 

36.915 

4.669 

0.000 

1. 3756E-01 1.0528E-01 

31 

30.006 

5.336 

0.000 

6.5651E— 02 6. 1 592E-02 

32 

31.391 

5.336 

0.000 

6. 5420E-02 6.7219E-02 

33 

32.776 

5.336 

0.000 

6.6716E-02 7.0462E-02 

^4 

34.161 

5.336 

0.000 

6.6593E-02 7.4232E-02 

35 

35.546 

5.336 

0.000 

7.0247E-O2 7.9488E-0? 

36 

36.931 

5.336 

0.000 

7.5330E-02 8. 1272E-02 

37 

31.546 

6.002 

0.000 

4.9274E-02 5.8759E-02 

38 

32.626 

6.002 

0.000 

4.9898E-02 5.9872E-02 

39 

33.706 

6.002 

0.000 

4.9317E— 02 6.0388E-02 

40 

34.786 

6.002 

0.000 

4.9281E— 02 6.2543E-02 

41 

35.866 

6.002 

0.000 

5.0686E— 02 6.4549E-02 

42 

36.946 

6.002 

0.000 

5.2233E— 02 6.4901E-02 

43 

33.079 

6.664 

0.000 

3.8066E-0? 5. 1788E-02 

44 

33.855 

6.664 

0.000 

3 . 7679E-02 5. 1 864E-02 

45 

34.632 

6.664 

0.000 

3. 7439E-02 5.2704E-02 

46 

35.408 

6.664 

0.000 

3.7757E-02 5.3717E-02 

47 

36.185 

6.664 

0.000 

3.8319E-02 5. 4325E-02 

48 

36.961 

6.664 

0.000 

3.8964E-02 5.4283E-02 

49 

34.592 

7.319 

o.coo 

2.9216E-02 4. 5475E-02 

50 

35.069 

7.319 

0.000 

2.9235E-02 4.5894E-02 

51 

35.546 

7.319 

0.000 

2.9378E— 02 4.6304E-02 

52 

36.023 

7.319 

0.000 

2.9593E— 02 4.6567E-02 

53 

36.499 

7.319 

0.000 

2.9813E-02 4.6650E-02 

54 

36.976 

7.319 

0.000 

3.0125E-02 4.6519E-02 

55 

35.968 

7.914 

0.000 

2.3844E-0 2 4.0931E-02 

56 

36.172 

7.914 

0.000 

2. 3916E-02 4. 0967E-02 

57 

36.377 

7.914 

0.000 

2.3987E-02 4.0991E-02 

58 

36.581 

7.914 

0.000 

2.4055E-02 4.1003E-02 

59 

36.785 

7.914 

0.000 

2.4127E-02 4. 0998E-02 

60 

36.990 

7.914 

0.000 

2.4282E-02 4.0881E-02 


(g) Concluded. 
Figure 7.- Continued. 
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NASA/AMES WING PROGRAM 


*****TAIL ***** 


CANARD - WING - BODY CONFIGURATION 


FLOW CONDITIONS - MACH 

2.01000 


BETA 

1.74359 


REFERENCE QUANTITIES - AREA 

1 92.000 


LENGTH 
15. 330 


XMOM 

25.000 


B/2 
8 • 360 


PANELS 

60 


CNP 

0.30934 

SPAN LOAD 
Y 

1.99578 
2.66439 
3.33281 
4.00097 
4.66870 
5. 33575 
6.00151 
6.66443 
7.31886 
7.91379 


DCN/DALPHA 

t R AD** — 1 ) ( DEG**- 1 ) * BETA CMP 

1.12533 0.01964 1.96210 -0.12993 


DCM/DCN 

-0.39058 


XAC 

30.988 


XCP 

31.439 


DISTRIBUTION 

CCN/2B 

0.15634 

0.15320 

0.15181 

0.15348 

0.15809 

0.15104 

0.13822 

0.11990 

0.09521 

0.04957 


2B * 33.4400 


TRAILING VORTEX .... 


Y = 7.3476 


GAM/V2PI 


4. 1602E-01 



SPANWISE DISTRIBUTION 

OF L.E. SUCTION 

REGION 

COL 

Y/ ( B/2 1 

CCS/2B 

EPS ( DEG ) 

CCX/2B 

1 

1 

0.23943 

0.00666 

86.565 

-0.00190 

l 

2 

0.31949 

0.00425 

92.983 

-0.00166 

1 

3 

0.39955 

0.00854 

104.707 

-0.00486 

1 

4 

0.47961 

0.01663 

93.895 

-0.00674 

l 

5 

0.55967 

0.02867 

92.705 

-0.01107 

1 

6 

0.63973 

0.04133 

90.049 

-0.01417 

t 

7 

0.71979 

0.05274 

89.082 

-0.01724 

1 

8 

0.79985 

0.06464 

88.450 

-0.02045 

1 

9 

0.87991 

0.07851 

87.837 

-0.02405 

l 

10 

0.95997 

0.08993 

90.052 

-0.03083 

SUMMARY OF VORTEX STRENGTHS AND POSITION AT I 


VORTEX GAMMA/2*PI*V 

Y 

l 



1 

0.074799 

4.16 2 

0.546 



2 

0.328235 

2.55? 

3.083 



3 

0.150326 

1.794 

2.933 



4 

0.416016 

7.348 

0.000 



5 

0.069799 

6.440 

1.411 


CCY/2B 

0.00638 

0.00391 

0.00702 

0.01521 

0.02645 

0.03882 

0.04984 

0.06132 

0.07474 

0.08448 

- 37.276 


(h) Page 8. 

Figure 7.- Continued. 


YCP 

4.599 
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ALPHA 

16.00 

M 

2.01 

— INCIDENCE — 
WING TAIL 

10.000 0.000 

W 

KVLE* 

0.500 

I NG 

KVSE* 

0.000 

TAIL 

KVLE* 

0.500 

KVSE* W/V 

0.000 1. 

LIMIT 

000 

SUMMARY DF FORCE AND PITCHING MOMENT COEFFICIENTS 







CN 

CM 

XCP 

CL 

CD I 

CA 

NOSE... 

POTENTIAL 

VORTEX 

6.402E-03 

2.777E-03 

9. 357E-03 
4.141E-03 

2.596E 00 
2.138E 00 

6.1546-03 

2.6696-03 

1 .7656-03 
7.6536-04 


WING... 

WtB) 

W(B> 

WtB) 

POTENTIAL 
VORTEX, LE 
VORTEX, SE 

6.2726-02 

2.531E-02 

0.000 

6. 555E-02 
2.723E-02 
0.000 

8.977E 00 
8.505E 00 
0.000 

5.637E-02 
2.275E— 02 
0.000 

1 .6166-02 
6.523E-03 
0.000 

1.106E-02 
4. 462E-03 
0.000 

BODY... 
B(W) 
B(W) 
B C W) 

POTENTIAL 
VORTEX, LE 
VORTEX, SE 

?• 3196-02 
9.3586—03 
0.000 

2.223E-02 

8.971E-03 

0.000 

1.030E 01 
1.030E 01 
1.030E 01 

2. 229E-02 
8.995E-03 
0.000 

6.392E-03 

2.579E-03 

0.000 


AFTERBODY... 

-6.246E— 03 

-4. 3846-03 

1.424E 01 

-6 .0046— 03 

- 1.722E-03 


TAIL... 
T(B) 
T ( 8) 
T ( B 1 

POTENTIAL 
VORTEX, LE 
VORTEX, SE 

3.0936-01 

4.537E-02 

0.000 

-1.2996-01 
-*1 .9566-02 
0.000 

3.1446 01 
3.1616 01 
0.000 

2.9746-01 

4.361E-02 

0.000 

8.527E-02 
1. 2516-02 
0.000 

0.000 

0.000 

0.000 

BODY... 

B (T ) 
8 ( T) 
BtTI 

POTENTIAL 
VORTEX, LE 
VORTEX, SE 

/ 

5.034E-02 

7.383E-03 

0.000 

-1 .6296-02 
—2 .3906-03 
0.000 

2.996E 01 
2.996E 01 
2.996E 01 

4.8396-02 

7.097E-03 

0.000 

1.388E-02 

2.0356-03 

0.000 


TOTAL CONFIGURATION 

... 5.3596-01 

— 3.508E-02 

2.600E 01 

5.0976-01 

1.461E-01 

l* 552E-02 



CDI/CL**2 - 

5.626E- 

01 





(i) Page 9. 
Figure 7.- Continued. 
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SUMMARY OF FREE VORTEX TRAJECTORIES 


VORTEX 1 BODY VORTEX FROM NOSE 

VORTEX 2 WING TRAILING VORTEX 

VORTEX 3 WING SEPARATION VORTEX 

VORTEX GAMMA/2*P1*V 

X 0.074799 

2 0.328235 

3 0.150326 


X DX 

A 

S 

RO 

DA/DX 

4.2750 1.0000000 

0.8260 

0.8260 

0.8260 

0.1169 

VORTEX SIGMA! REAL ) 

( IMAG) 




X 5.332E-01 1. 

083E 00 




X OX 

A 

S 

RO 

OA/DX 

5.2750 0.5000000 

0.9269 

1.4015 

1.0072 

0.0963 

VORTEX SIGMAIREAU 

I IMAG) 




1 5.I71E-01 1. 

276E 00 




X DX 

A 

S 

RO 

OA/OX 

6.2750 0.5000000 

1.0185 

1.9118 

1.2272 

0.0776 

VORTEX SIGMAIREAU 

I IMAG) 




I 4.963E-01 1. 

*4326 00 




X DX 

A 

S 

RO 

DA/DX 

7.2750 0.5000000 

1.0785 

2.1977 

1.3635 

0.0649 

VORTEX SIGMAIREAU 

I IMAG) 




1 4.731E-01 1. 

p562E 00 





• • 

• 



X DX 

A 

S 

RO 

DA/DX 

11.7750 0.5000000 

1.3324 

1.3324 

1.3324 

0.0539 

VORTEX SIGMA(REAL) 

I IMAG) 




1 3.732E-01 1. 

»968E 00 




X DX 

A 

S 

RO 

OA/DX 

11.7760 0.5000000 

1.3325 

1.3325 

1.3325 

0.0539 

VORTEX SIGMAIREAU 

{ IMAG) 




1 3.732E-0 1 1. 

. 968E 00 




2 2.822E 00 -3 

.907E-01 




3 2.595E 00 3, 

.603E-01 




X DX 

A 

S 

RO 

DA/DX 

12.7760 0.2500000 

1.3885 

1.3885 

1.3885 

0.0496 

VORTEX SIGMAIREAU 

1 IMAG) 




1 3.622E-01 1 

•984E 00 




2 2.9S8E 00 -9 

•758E-02 




3 2.321E 00 2 

• 655E-01 





• II 

(j) page 10. 

Figure 7.- Continued. 
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x nx a 

24.7760 0. 5000000 1.6670 

s 

3*911 l 

RO 

2.3108 

DA/OX 

0.0000 

VORTEX SIGMA(PEAL) ( I MAG ) 

1 1.717F 00 6.512E-01 

2 2.551E 00 7.158F 00 

3 3.282E 00 1.953E 00 





X OX A 

25.7760 0.5000000 1.6670 

s 

4*2751 

RO 

2.4625 

DA/DX 

0.0000 

VORTEX SIGMA(PEAL) (IMAG) 

1 1.875E 00 3.794E-01 

2 2.544F 00 2.087F 00 

3 3. 1 55F 00 2.528E 00 





X OX 

26.7760 0.5000000 

A 

1.6670 

s 

4.6390 

RO 

2.6190 

DA/DX 

0.0000 

VORTEX SIGMA(REAL) (IMAG) 

1 2. 05 IE 00 2.585E-01 

2 7.659F 00 2.096E 00 

3 2.761F 00 2.879E 00 





• • 

• 



X OX 

36.7760 0.5000000 

A 

1.6670 

S 

8*2785 

RO 

4.3071 

DA/DX 

0.0000 

VOPTF X SIGMA(REAL) 

1 4.116E 00 

2 2.565F 00 

3 1.804E 00 

( IMAG) 
4.262E-01 
2.927F 00 
3.094E 00 





X OX A 

37.0260 0.2500000 1.6670 

s 

1.6670 

RO 

1.6670 

DA/DX 

0.0000 

VDPTEX SI GMAI REAL ) (IMAG) 

1 4 . 1 4 OF 00 4.597E-01 
7 7.562F 00 2.992E 00 
3 1.789E 00 3.004E 00 




X DX A 

37.1510 0.1250000 1.6670 

S 

1.6670 

RO 

1*6670 

DA/DX 

0.0000 

VORTEX SIGMAIREAL) (IMAG) 

1 4.151E 00 5.068E-01 
7 7.558F 00 3.039E 00 
3 1.789F 00 7.970E 00 





(j) Concluded. 
Figure 7.- Concluded 
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